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Abstract 


VNI is a general-purpose Monte-Carlo event-generator, which includes the simulation of lepton-lepton, 
lepton-hadron, lepton-nucleus, hadron-hadron, hadron-nucleus, and nucleus-nucleus collisions. On the basis 
of renormalization-group improved parton description and quantum-kinetic theory, it uses the real-time 
evolution of parton cascades in conjunction with a self-consistent hadronization scheme that is governed by 
the dynamics itself. The causal evolution from a specific initial state (determined by the colliding beam 
particles) is followed by the time-development of the phase-space densities of partons, pre-hadronic parton 
clusters, and final-state hadrons, in position-space, momentum-space and color-space. The parton-evolution 
is described in terms of a space-time generalization of the familiar momentum-space description of multipl 
(semi) hard interactions in QCD, involving 2 —> 2 parton collisions, 2^1 parton fusion processes, and 1 —> 2 
radiation processes. The formation of color-singlet pre-hadronic clusters and their decays into hadrons, on 
the other hand, is treated by using a spatial criterion motivated by confinement and a non-perturbative 
model for hadronization. This article gives a brief review of the physics underlying VNI, which is followed 
by a detailed description of the program itself. The latter program description emphasizes easy-to-use 
pragmatism and explains how to use the program (including a simple example), annotates input and control 
parameters, and discusses output data provided by it. 


Program obtainable from: http://rhie.phys.Columbia.edu/rhie/vni 
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1 PROGRAM SUMMARY 


Title of program: VNI-3.1 

Computer fo which the program has been designed and others on which it has been tested: 

IBM RS-6000, Sun Sparc, Hewlett Packard UX A-9000 

Operating systems: IBM-AIX, Sun-OS, and any other UNIX operating systems, as well as LINUX. 
Programming language used: Fortran 77 

Memory required to execute with typical data: 2000 kwords 
No. of bits in a word: 32 

No. of lines in distributed program: 25760 lines of main program, plus 49 and 244 lines for two example 
programs. 

Keywords: Monte Carlo simulation, event generator, QCD kinetic theory, parton cascades, parton coales¬ 
cence, hadronic final states. 

Nature of physical problem: 

In high-energy particle collisions certain phase-space regions can be populated by a large number of quanta, 
such that statistical correlations among them (e.g., in space, momentum, or color) become of essential im¬ 
portance. Examples are deep-inelastic lepton-hadron scattering and hadron-hadron collisions in the region 
of very small Bjorken-x, or, collisions involving heavy nuclei in the central rapidity region. In these cases the 
produced particles evolve in a complicated non-equilibrium environment created by the presence of neigh¬ 
boring ones. The ‘deterministic’ quantum evolution of particle states due to self-interactions (depending 
only on the particle itself), receives a new ‘statistical’ kinetic contribution due to mutual interactions (de¬ 
pending crucially on the local density). The theoretical basis for addressing the solution for the dynamics 
of such particle systems is a quantum-kinetic formulation of the QCD equations of motion, an approxima¬ 
tion that combines field-theoretical aspects associated with the renormalization group (including well-known 
resummation techniques) with aspects of transport theory associated with non-equilibrium multi-particle 
dynamics (including important quantum effects beyond the classical level). 

Method of solution: 

The solution of the underlying quantum-kinetic equations of motion for non-equilibrium multi-particle QCD 
by Monte-Carlo simulation of collisions allowing for a variety of combinations of beam and target particles. 
To simulate the real-time evolution of the collision system in position space and momentum space on the 
basis of the equations of motions, the procedure is three-fold: i) the construction of the initial state including 
the decomposition of the beam and target particles into their partonic substructure, (ii) the evolution of 
parton cascades including multiple scatterings, emission- and fusion-processes, and (iii) the self-generated 
conversion of partons into hadrons using a phenomenological model for parton-coalescence into pre-hadronic 
clusters and subsequent decay into final-state particles. 

Restriction on the complexity of the problem: 

For very high collision energy (y^ ^ 10 TeV in hadronic collisions, and ^/s ^ 5 TeV/nucleon in nuclear 
collisions) numerical inaccuracies due to repeated Lorentz boosts, etc., may accumulate to cause problems. 
Although the most concerned parts of the program use double precision, for extreme energies the code would 
require a conversion in full to double precision format (which is planned in the near future). 

Typical running time: 

The CPU time for a typical simulation is strongly dependent on the type of beam and target, the magnitude 
of collision energy, as well as on the time interval At chosen to follow an event in its real-time evolution. 
Examples are (for At = 35 fm): a) -|- e~ at y/s = 100 GeV: 10000 events/hour; h) p + p at y/s = 200 
GeV: 5000 events/hour; c) p + Au at y/s = 200 GeV/nucleon; 100 events/hour; d) Au + Au at y/s = 200 
GeV/nucleon: 1 event/hour; All of the above quotes are approximate, and refer to a typical 133 Mhz or 166 
Mhz processor on a modern Power-Workstation or Power-PC. 
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LONG WRITE-UP 


2.1 Introduction 

VNI [| is the Monte Carlo implementation of a relativistic quantum-kinetic approach ii to the dynamics 
of high-energy particle collisions, inspired by the QCD parton picture of hadronic interactions i, ii- It 
is a product of several years of development in both the improving physics understanding of high-energy 
multiparticle dynamics in QCD, as well as the technical implementation in the form of a computer simulation 
program. The most relevant references for the following are Refs. g i I, § 0, 0] , where details of the 
main issues, discussed below, can be found. The puropse of VNI is to provide a comprehensive description 
of particle collisions involving beams of leptons, hadrons, or nuclei, in terms of the space-time evolution 
of parton cascades and parton-hadron conversion. The program VNI is concepted as a useful tool (and 
nothing more) to study the causal development of the collision dynamics in real time from a specified initial 
state of beam and target particles, all the way to the final yield of hadrons and other observable particles. 
The collision dynamics is traced in detail on the microscopic level of quark and gluon interactions in the 
framework of perturbative QCD, supplemented by a phenomenological treatment of the non-perturbative 
dynamics, including soft parton-collisions and parton-hadronization. The generic structure of the simulation 
concept is illustrated in Fig. 1. 

The main strength of VNI lies in addressing the physics of high-density QCD, which becomes an in¬ 
creasingly popular object of research, both from the experimental, phenomenological interest, and from the 
theoretical, fundamental point of view. Presently, and in the near future, the collider facilities HERA (ep, 
eA?), Tevatron (pp, pA), RHIC and LHC (pp, AA) are able to probe new regimes of dense quark-gluon matter 
at very small Bjorken-x or/and at large A, with rather different dynamical properties. The common feature 
of high-density QCD matter that can be produced in these experiments, is an expected novel exhibition 
of the interplay between the high-momentum (short-distance) perturbative regime and the low-momentum 
(long-wavelength) non-perturbative physics. For example, with HERA and Tevatron experiments, one hopes 
to gain insight into problems concerning the saturation of the strong rise of the proton strncture functions at 
small Bjorken-x, possibly due to color-screening effects that are associated with the overlappping of a large 
number of small-x partons. Another example is the anticipated formation of a quark-gluon plasma in RHIC 
and LHC heavy ion collisions, where multiple parton rescattering and cascading generates a high-density 
environment, in which the collective motion of the quanta imust be taken into consideration. In this context 
the most advantageous and novel feature of VNI is the space-time cascade picture that provides a potentially 
powerful tool to study high-density particle systems in QCD, where accounting for the dynamically evolving 
space-time structure of the collisions is most important. The necessity of including space-time variables in 
addition to energy-momentum variables for high-density systems may be common knowledge in the field 
of heavy-ion physics, where the space-time aspect of many-body transport theory is essential ingredient to 
describe nucleus-nucleus collisions, but it is new to many high-energy particle physicists which only recently 
began to acknowleged the necessity to include time and space on top of the commonly used, pure momentum 
space description of, e.g., lepton and hadron collisions. 


Where does VNI fit into the diverse family of modern event generators for physics simulation of particle 
collisions (where ‘particle’ stands for any beam/target particle from an electron to a heavy nucleus)? 

The Monte Carlo models for high-energy particle collisions that are on the market, may be crudely divided 
into two classes: 


(i) The first class embodies event generators that restrict to ’’clean” reactions involving lepton or proton 
beams only, and which aim at a high-precision description of experimental tests of QCD’s first prin¬ 
ciples. Popular examples are JETSET/PYTHIA HERWIG g, ARIADNE @, LEPTO ||, 
and IS A JET [^|. The common feature of these event generators is a combination of well understood 
perturbative QCD parton-shower description and a non-perturbative hadronization prescription to 
convert the partonic final state into hadrons. 


^ The three letters VNI do not mean anything profound. VNI is pronounced ’’Vinnie”, short for ’’Vincent Le CuCurullo Con 
GiGinello”, the little guy who likes to hang out with his pals, the quarks (crtcrtrullos) and gluons (jiginellos). That is QGD in 
Wonderland, and that is the whole, true story. 


3 







VNI 



DATA 


Figure 1: Generic flow-chart of the simulation concept of VNI: It starts with the initial beam target 
particles A and B, decomposing them (except for leptonic A and/or B) in their hadronic constituents with 
partonic substructure, then proceeds through the parallely evolving stages of parton-cascade, parton-hadron 
conversion, and fragmentation of beam/target remnants, and finishes up with the final particle yield that 
reflects observables in a detector. 
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(ii) The second class are event generators that aim to describe ’’dirty” reactions involving nuclei, much less 
based on first-principle knowledge, but instead rely on phenomenological models to mimic the unknown 
details of the underlying physics. Here the widely used concept is to visualize a nuclear collision in terms 
of nucleon-nucleon collisions on the basis of a constituent valence-quark picture plus string-excitation 
and -fragmentation. Examples for these models are FRITIOF [|^], DPM [^], VENUS |Q, RQMD 
|49|, and HIJET Distinct from these is HIJING |51], which also incorporates a perturbative QCD 
approach to multiple minijet production, however it does not incorporate a space-time description. 


With the exception of HERWIG, all of the above Monte Carlo generators utilize some form of string frag¬ 
mentation phenomenology to model the non-perturbative hadronization and final-state particle production. 
Most commonly used is the Lund string model HERWIG on the other hand is built on a very dif¬ 
ferent parton-cluster formation/fragmentation approach, which forms the basis of the hadronization scheme 
developed in VNI. 

What are the shortcomings of the above-mentioned Monte Carlo models with respect to the particle dynamics 
in finite-density regions created by high-energy collisions? 

The high-energy particle physics generators of the first class lack the inclusion space-time variables in 
the dynamics description in both the perturbative QCD parton evolution and the non-perturbative process 
of parton-hadron transition. These models therefore cannot account for statistical interactions due to the 
presence of a finite density of particles close-by in space, such as rescattering, absorption of recombination 
processes. Hence, although these models to large extent use QCD’s fundamental quark-gluon degrees of 
freedom, important aspects of parton dynamics and interactions at finite density are left out, because the 
particles are assumed to propagate unscathed in free space. 

The event generators of the second class, on the other hand, mostly do utilize a space-time description, 
however on the level of hadronic degrees of freedom (strings, baryons, mesons and resonances) rather than 
partonic degrees of freedom. For ultra-relativistic nucleus-nucleus collisions the parton approach appears 
to be more realistic than the hadronic or string picture, as it has been realized that short-range parton 
interactions play a major role for heavy ion collisions at collider energies of ^/s ^ 100 GeV, at least during 
the early and most dissipative stage of the first few fm. Here copiously produced quark-gluon mini-jets 
cannot be considered as isolated rare events, but must be embedded in complicated multiple cascade-type 
processes. Thus, the short range character of these interactions implies that perturbative QCD can and must 
be used, and that the picture of comparably large distance excitations of strings or hadronic resonances does 
not apply in this kinematic regime. 

In view of this discussion, the program VNI can be viewed in between the above two classes of Monte 
Carlo models: It provides a kinetic space-time description of parton evolution by utilizing well-developed 
techniques for perturbative QCD simulations at zero density or free space, as the event generators of the 
hrst class. On the other hand, it applies this concept also to the physics of finite-density particle systems, 
e.g., in collisions involving nuclei, as the event generators of the second class. Comparing VNI with the 
above-mentioned Monte Carlo models, the essential differences and partly new aspects are the following: 


a) the aspects of the space-time evolution of the particle distributions in addition to the evolution in 

momentum space Q , and the concepts of quantum kinetic theory and statistical physics |]|, . 

b) the self-consistent interplay of coherent (angular ordered) perturbative parton evolution according to the 

DGLAP 1^] equations, with the fully dynamical cluster-hadronization according to the phenomeno¬ 
logical model of Ref. P] . 


c) the microscopic tracing of color degrees of freedom and the effect of color-correlations by using explicit 

color-labels for each parton, which allows to investigate final state interactions in the process of hadron 
formation 0. 

d) the diverse advantages of a stochastic simulation technique with which the various particle interaction 

processes are determined by the dynamics itself, through the local density of particles as they evolve 
causally in time. 


the statistical many-particle description for general non-equilibrium systems, which allows to study 
thermodynamic behaviour of the bulk matter |52|, such as the evolution of macroscopic energy density. 
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pressure, etc., or the dynamical development of the sytem to thermal/chemical equilibrium in heavy 
ion reactions. 

In summary, the improvement to be expected from VNI for the physics simulation of high-energy par¬ 
ticle collisions lies clearly in the ‘dirty’ high-density parton regime, where the space-time aspects are most 
important, and which currently and in the future is of central interest in experiments at, e.g., HERA, RHIC 
and LHC. On the other hand, VNI may also be perceived as a valuable alternative to high-energy event 
generators for the study of ‘clean’, zero-density collisions as e'’'e“ annihilation or pp collisions, where the 
space-time aspects can provide useful additional insight in the collision dynamics for experiments at, e.g., 
LEP or the Tevatron. 


2.2 General concept 

The central element in the physics description implemented in VNI is the use of QCD transport theory ||8| 
and quantum field kinetics |l| to follow the evolution of a generally mixed multi-particle system of partons 
and hadrons in 7-dimensional phase-space d^rd?kdE. Included are both the parton-cascade development 
i, & 0 0] which embodies the renormalization-group improved evolution of multiple parton collisions 
including inelastic (radiative) processes, and the phenomenological parton-hadron conversion model of Refs, 
iiilll , in which the hadronization mechanism is described in terms of dynamical parton-cluster formation 
as a local, statistical process that depends on the spatial separation and color of nearest-neighbor partons, 
followed by the decay of clusters into hadrons. 

In contrast to the commonly-used momentum-space description, the microscopic history of the 
dynamically-evolving particle system is traced in space-time and momentum space, so that the correla¬ 
tions of partons in space, time, and color can be taken into account for both the cascade evolution and 
the hadronization mechanism. It is to be emphasized, that the interplay between perturbative and non- 
perturbative regimes is controlled locally by the space-time evolution of the mixed parton-hadron system 
itself (i.e., the time-dependent local parton density), rather than by an arbitrary global division between 
parton and hadron degrees of freedom (i.e., a parametric energy/momentum cut-off). In particular the 
parallel evolution of the mixed system of partons, pre-hadronic clusters, and hadrons, with the relative 
proportions determined by the dynamics itself, is a novel feature that is only possible by keeping track of 
both space-time and energy-momentum variables. 

Probably the greatest strength of this approach lies in its application to the collision dynamics of com¬ 
plicated multi-particle systems, as for example in collisions involving nuclei {eA, pA and AB), for which a 
causal time evolution in position space and momentum space is essential: Here statistical, non-deterministic 
particle interactions are most important, which can only be accounted for by following the time-evolution 
of the particle densities in space and energy-momentum. This approach allows to study the time evolution 
of an initially prepared beam/target collision system in complete phase-space from the instant of collisional 
contact, through the QCD-evolution of parton distributions, up to the formation of final hadronic states. 
It provides a self-consistent scheme to solve the underlying equations of motion for the particle densities as 
determined by the microscopic dynamics. 

The model as a whole consists of three major building-blocks, each of which is illustrated schematically 
in Eig. 2: 

1. The initial state associated with the incoming collision partners (leptons, hadrons, or nuclei). Except 

for lepton beams, this involves the phenomenological construction of hadrons or nuclei in terms of 
the partons’ phase-space distributions on the basis of the experimentally measured hadron (nucleon) 
structure functions and elastic form-factors. 

2. The parton cascade development with mutual- and self-interactions of the system of quarks and gluons. 

This includes multiple inelastic processes, described as sequences of elementary 2^2 scatterings, 
1 —> 2 emissions and 2^1 fusions. Moreover, correlations are accounted for between primary 
virtual partons, emerging as unscathed remainders from the initial state, and secondary real partons, 
materialized or produced through the partonic interactions. 
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3. The hadronization dynamics of the evolving system in terms of parton-coalescence to color-neutral clusters 
as a local, statistical process that depends on the spatial separation and color of nearest-neighbor 
partons. Each pre-hadronic parton-cluster fragments through isotropic two-body decay into primary 
hadrons, according to the density of states, followed by the decay of the latter into final stable hadrons. 


Such a pragmatical division, which assumes complex interference between the different physics regimes to 
be negligible, is possible if the respective dynamical scales are such that the short-range (semi)hard parton 
interactions (scattering, radiation, fusion) of perturbative nature, and the non-perturbative mechanism 
of hadron formation (parton-coalescence and cluster-decay), occur on well-separated space-time scales (or 
momentum scales, by virtue of the uncertainty principle). Loosely speaking, the typical momentum scale 
associated with parton collisions, radiative emissions, or parton fusion, has to be larger than the inverse 
‘conhnement length scale’ ~ 1 fm which separates perturbative and non-perturbative domains. Further 
discussion of this condition of validity can be found in, e.g., |^, ]d|. 


2.3 Equations of motion from qnantum kinetics for multi-particle dynamics 

A space-time description of multiparticle systems in high-energy QCD processes can be derived systemati¬ 
cally from quantum-kinetic theory on the basis of QCD’s first principles in a stepwise approximation scheme 
(see e.g., Refs. ©I and references therein). Applied to the concept sketched in the preceeding subsec¬ 
tion, this framework allows to cast the time evolution of the mixed system of individual partons, composite 
parton-clusters, and physical hadrons in terms of a closed set of integro-differential equations for the phase- 
space densities of the different particle excitations. The dehnition of these phase-space densities, denoted 
by Fa, where a = p,c,h labels the species of partons, pre-hadronic clusters, or hadrons, respectively, is: 

. F^(t,nE,k) = , (1) 

where = E'^ — k'^ can be offshell or on-shell, as will be discussed in the following subsections. The 
densities (||) measure the number of particles of type a at time t with position in r-|- dr, momentum in 
k + dk, and energy in E + dE (or equivalently invariant mass in k‘^ + dk‘^). The Fa are the quantum analogues 
of the classical phase-space distributions, including both off-shell and on-shell particles, and hence contain 
the essential microscopic information required for a statistical description of the time evolution of a many- 
particle system in complete 7-dimensional phase-space d^rd^kdE, thereby providing the basis for calculating 
macroscopic observables. 

The phase-space densities (||) are determined by the self-consistent solutions of a set of transport equa¬ 
tions (in space-time) coupled with renormalization-group-type evolution equations (in momentum space). 
Referring to Refs. [^, ^ for details, these equations can be generically expressed as convolutions of the den¬ 
sities Fa of particle species a, interacting with specific cross sections Ij for the processes j. The resulting 
coupled equations for the space-time development of the densities of partons Fp, clusters Fc and hadrons 
Ffi is a self-consistent set in which the change of the densities Fa is governed by the balance of the various 
possible interaction processes among the particles. Fig. 3 represents these equations pictorially. For the 
densities of partons the transport equation (governing the space-time change with r^) and the evolution 
equation (controlling the change with momentum scale k^), read, respectively. 
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For the densities of pre-hadronic clusters and hadrons, the evolution equations are homogeneous to good 
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a) initial state 



b) parton cascade development 
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partons 

c) parton-clnster formation & liadronization 

central 



beam cluster B 

parton clusters 


Figure 2; The three components of the model: a) the initial state constructed in terms of the parton 
distribution of the incoming nuclei; b) the time-evolution of parton cascades in 7-dimensional phase-space 
c) the formation of color neutral clusters from secondary partons emerging from cascading, as well of the 
remnant primary partons from the initial state, followed by the fragmentation of the clusters into final 
hadrons. 8 



approximation, so that one is left with non-trivial transport equations only, ^ 
Fc{r,k) = FpFpf o I{pp ^ c) - Fco i{c ^ h) , 
kfj,-^ Fh{r,k) = Fcoi{c^h) + Fh'o i{h' ^ h) - Fuo i{h ^ h') 




0 (4) 

0 (5) 


In (§)-(P), each convolution F o / of the density of particles F entering a particular vertex I includes a 
sum over contributing subprocesses, and a phase-space integration weighted with the associated subprocess 
probability distribution of the squared amplitude. Explicit expressions are given in Refs. |^, |9|. 

The terms on the right-hand side of the transport- and evolution-equations @-(||) corresponds to one 
of the following categories: 


(i) parton scattering and parton fnsion through 2-body collisions. 


(ii) parton mnltiplication through radiative emission processes on the perturbative level. 


(iii) colorless cluster formation through parton recombination depending on the local color and spatial 
confignration. 


(iv) hadron formation throngh decays of the clnster excitations into final-state hadrons. 

As mentioned before, the eqnations ®-(i reflect a probabilistic interpretation of the multi-particle 
evolution in space-time and momentum space in terms of sequentially-ordered interaction processes j, in 
which the rate of change of the particle distribntions Fa {a = p, c, h) in a phase-space element d?rd^k is 
governed by the balance of gain (-I-) and loss (—) terms. The left-hand side describes free propagation 
of a quantnm of species a, whereas on the right-hand side the interaction kernels I are integral operators 
that incorporate the effects of the particles’ self- and mutual interactions. This probabilistic character is 
essentially an effect of time dilation, because in any frame where the particles move close to the speed of 
light, the associated wave-packets are highly localized to short space-time extent, so that comparatively 
long-distance qnantum interference effects are generally very small. 


2.4 Scheme of solution in global Lorentz-frame of reference 

In the above kinetic approach to the mnlti-particle dynamics, the probabilistic character of the transport- 
and evolution equations (§)-(^) allows one to solve for the phase-space densities Fa{r,k) by simulating the 
dynamical development as a Markovian process causally in time. Because it is an ‘initial-value problem’, 
one must specify some physically appropriate initial condition Fa{to,f, k) at starting time to, such that all 
the dynamics prior to this point is effectively embodied in this initial form of Fa- The set of equations 
(I)"® can then be solved in terms of the evolution of the phase-space densities Fa for t > to using Monte 
Carlo methods to simulate the time development of the mixed system of partons, clnsters, and hadrons in 
position and momentnm space |^, P]. 

With the initial state specified as discnssed below, the phase-space distribution of particles at t = to = 0 
can be constructed and then evolved in small time steps At = 0(10“^ fm) forward throughout the par- 
allely evolving stages of parton cascade, parton-cluster formation, and cluster-hadron decays, nntil stable 
final-state hadrons and other particles (photons, leptons, etc.) are left as freely-streaming particles. The 
partons propagate along classical trajectories until they interact, i.e., collide (scattering or fusion process), 
decay (emission process) or recombine to pre-hadronic composite states (cluster formation). Similarly, the 
so-formed pre-hadronic parton-clusters travel along classical paths until they convert into primary hadrons 
(cluster decay), followed by the hadronic decays into stable final state particles. The corresponding proba¬ 
bilities and time scales of interactions are sampled stochastically from the relevant probability distributions 
in the kernels I of eq. (|)-®- 

^ It is worth noting that eq. (^) embodies the momentum space (fe^) evolution of partons through the renormalization of 
the phase-space densities Fp, determined by their change k^dFp(r,k)/dk^ with respect to a variation of the mass (virtuality) 
scale k^ in the usual QCD evolution framework [^, ^ |^. On the other hand, for pre-hadronic clusters and hadrons, renormal¬ 
ization effects are comparatively small, so that their mass fluctuations can be ignored to first approximation, implying 

k^dFcir,k)/dk^ = k^dFh(r,k)/dk^ = 0 . 
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Figure 3: Graphical representation of the equations (2)-(5) for the particle phase-space densities Fp of 
partons, Fc of pre-hadronic clusters, and Fh of hadrons. 
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It is important to realize, that the spatial density and the momentum distribution of the particles are 
intimately connected: The momentum distribution continously changes through the interactions and deter¬ 
mines how the quanta propagate in coordinate space. In turn, the probability for subsequent interactions 
depends on the resulting local particle density. Consequently, the development of the phase-space densities 
is a complex interplay, which - at a given point of time - contains implicitely the complete preceeding history 
of the system. 

It is clear that the description of particle evolution is Lorentz-frame dependent, and a suitable refer¬ 
ence frame (henceforth called global frame) must be chosen (not necessarily the laboratory frame of an 
experiment). When computing Lorentz-invariant quantities, such as cross sections or final-state hadron 
spectra, the particular choice is irrelevant, whereas for non-invariant observables, such as energy distibu- 
tions or space-time-dependent quantities, one must at the end transform from the arbitrarily-chosen frame 
of theoretical description to the actual frame of measurement. 

For calculational convenience, it is most suitable to choose the global center-of-mass (cm) frame of the 
colliding beam particles, with the collision axis in the z-direction. In this global cm-frame, the incoming 
particles A, B (= lepton, hadron, or nucleus) have four-momenta, 

^A,B = {EA,B,0,0,±Pcm) 

Ea,b = ^ (s + Ml ± Ml) , = El^ = (Pl + P|)' (6) 

Pcm = ^ - (Ma + Mb)2] [s - {Ma - Mb)2] , 

where Ma,b are the masses of A and B, and so incoming particles then carry well-defined fractions of Pcm, 
having only a non-vanishing longitudinal momentum along the z-axis. Particularly in the case of nuclei A, 
the daughter nucleons Ni = 1, ... A, have momenta, Pat. = (0,0, ±Pcm/A). 

In the following, the global cm-frame of A and B is assumed to be the reference frame, with the initial 
energy-momentum of the collision system given by (^). Furtermore the terms ‘hadron’, respectively ‘nucleon’, 
are used to distinguish initial states A B with A and/or B being a single hadron, respectively a nucleus 
with A (B) nucleons. 


2.5 Initial state 

If one or both beam particles A, B are leptons, they are considered as point-like objects which carry 
the full beam energy, meaning that any QED or QCD substructure of the leptons, as well as initial-state 
photon radiation by the leptons, is neglected. For lepton-lepton annihilation, it is assumed that the colliding 
leptons produce a time-like 7 or boson of invariant mass = -\-q^ at time t = —Q~^, so that t = 0 
characterizes the point when the 7 (Z^) decays into a quark-antiquark pair. Similarly, for lepton-hadron 
(nucleus) collisions, the lepton is emitting a space-like virtual 7 of invariant mass = —q^ at time 
t = —Q~^, and hence t = 0 is the point when the 7 hits the hadron (nucleus). 

In the general case of collisions involving hadrons and/or nuclei, the incoming particles A and B are 
decomposed into their parton substructure by phenomenological construction of the momentum and spatial 
distributions of their daughter partons on the basis of the known hadron (nucleon) structure functions and 
elastic hadron (nucleon) form-factors. In the cm-frame, where the two incoming particles A, B (= hadron, 
nucleus), are moving close to the speed of light, the parton picture is applicable and the parton substructure 
of the hadrons or nucleons can be resolved with reference to some initial resolution scale Qq. This resolution 
scale generally varies with beam-energy and mass number A, B,m. that it depends on the typical momentum 
and spatial density of partons as well as on their interaction probability. To be specific, it may be identified 
with the statistically estimated expectation value for the interaction scale of all primary parton-parton 
collisions (that are those in which at least one initial state parton is involved) [p(^ ]: 

g§ ^ Ql{.,P,A) = .4“ + kE^) ‘ , (7) 


11 







where {p\)prim is the average relative transverse momentum squared generated by primary parton-parton 
collisions, and Rq = lGeV~^. The pre-factor ^4“ (0<a<4/3a parameter) accounts for the nuclear 
dependence for A,B > 1, and x is the parton’s momentum fraction of the mother hadron or nucleon which 
carries momentum P (= Pcm for single hadron, or = Pcm/A for a nucleon in a nucleus A). The scale Qq 
defines the initial point in momentum space above which the system of partons is treated as an ensemble 
of incoherent quanta. The dynamics prior to this point is contained in the initial parton-structure function 
of the mother hadron (nucleon). Clearly, the convention (|^ yields only an average value dominated by the 
most probable (semi)hard parton collisions with relatively small momentum transfers of a few GeV. However, 
primary parton collisions with a momentum scale 3> Qq, which correspond to relatively rare fluctuations, 
are taken into account individually by the Q^-evolution of the hadron (nucleon) structure functions through 
space-like and time-like radiation processes, discussed later. 

The actual form of the phase-space distribution of partons, eq. (^ is initially to be specified at the 
reference point Qq. It is constructed as the following superposition of the parton distributions in the 
individual hadrons at Qq and at time t = to = 0, the point of collision of A and B: 


Nh 


Fa{r,k)lo=t, = Y.Pa^(k,P;Ql)-R^^{k,r,R) 


i=l 


r°=to 


( 8 ) 


The right hand side is a sum over all 1 < Ni < A^/j hadrons or nucleons contained in the collision system 
A + B. Each term in the sum is a parton phase-space density given by a convolution of an initial momentum 
distribution and a spatial distribution The subscript a = g, qj, qj (j = 1,... ,nf) labels the species 
of partons, and Ni refers to the type of the i-th hadron (for nuclei this is just proton or neutron). The 
four-vectors k = k^^ = {E, k) and r = = {t, r) refer to the partons. The 3-momenta and initial positions 

of the hadrons or nucleons are denoted P, respectively R. All vectors are understood with respect to the 
global cm-frame, at time t = to = t). The partons’ energies E = Ea{k‘^,q^) = \Jk'^ + -|- q^ take into 

account their initial space-like virtnalities < 0 which are distributed around (|g^|) = —Q'q/^, under the 
constraint that for each hadron in the initial state, the total invariant mass of the daughter partons must 
equal the mother hadron mass. This mimics the fact that the initial partons are confined inside their parent 
hadrons or nucleons and cannot be treated as free particles, i.e. they have not enough energy to be on mass 
shell, but are off mass shell by a space-like virtuality q^ (eq. (|l^) below). 

As suggestively illustrated in Fig. 2 and discussed in more detail now, the initial state A + B involving 
hadrons and/or nuclei, appears in the global cm-frame as two approaching ‘tidal waves’ of large-x partons 
(mainly valence quarks), where each ‘tidal wave’ has an extended tail of low—x partons (mostly gluons and 
sea-quarks). 


2.5.1 Initial momentum distribution 

For each hadron or nucleon the number of partons, the distribution of the flavors, their momenta and 
associated initial space-like virtualities, are obtained from the function P^^ in (^). Denoting P = Pcm/Nh^ 
with Nh the total number of hadrons (nucleons) in A -|- H, it is represented in the form 


P^^{k,P-Ql) = g) P^^{x,Ql) p^{x) g{k^) 6{P, - P) d^[Pj_) 


(9) 


with the momentum and energy fractions x = kz/P, x = E/P = x"^ + (A:^ -|- -|- q^)/P"^) and the 

normalizations 


p\ poo 

^ / dx xF^^{x,Ql) = 1 , / d?kxg{ki_) = 1 

„ Jo Jo 


E 


dk^ d?k 


EP^^{k\k,P-Ql) = n^^^iQlP) . 


0 


(27r)32F 

The physics behind the ansatz (I) is the following: 


( 10 ) 

( 11 ) 
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(i) The functions F^^{x, Qq) are the scale-dependent measured hadron (nucleon) structure functions with x 
being the fraction of the parent hadron’s or nucleon’s longitudinal momentum P carried by the parton. 
The transverse momentum distribution g{kx) is specified below. The factor xjx'm. ( p!^) is included to 
form the invariant momentum integral f d^k/[{2'K)'^2E] out of the distribution ||16|| . In (|Tl] ) the 
quantity is the total number of partons in a given nucleon with momentum P = |P| at Qq. 


(ii) The function p^{x) in takes into account nuclear shadowing effects affecting mainly soft (small x) 
partons in a nucleus |^. These shadowing effects are evident in the observations of the European Muon 
Collaboration as a depletion of the nuclear structure functions at small x relative to those of a free 
nucleon. Several mechanisms have been proposed to explain this nuclear shadowing effect on the basis 
of the parton model |^, EM- However, here instead the phenomenological approach of Wang and 
Gyulassy |^] is adopted which is based on the following parametrization for the A dependence of 
the shadowing for both quarks and gluons; 


7AM 


Paix) = 


{x,Q, 


Hp5(x,Q2) 


\l/6 . 

1 + 1.19 ( InAj 


/3a(P±)- , . . Vx 


x^ — 1.5(xo + xl)x‘^ + Sxqxlx 
exp(-xVa:o) ■ 


( 12 ) 


ln(H + 1) 


where F^ represent the structure function of a nucleus A, respectively those of a free 

nucleon, and xq = 0.1, x^ = 0.7. The function f^Air) = 0.1 ^ A^/^ — 1^ | y^l — ^ takes into account 

the impact parameter dependence, with R± labeling the transverse distance of a nucleon from its 
nucleus center and Ra the radius of the nucleus. 


(iii) The distribution g{k±) specifies the primordial transverse momenta of partons according to a normal¬ 
ized Gaussian 


1 




lUI' 


/Cq 


(13) 


independent of the type of parton or hadron (nucleon). It takes into account the uncertainty of 
momentum (Fermi motion) due to the fact that the initial partons are confined within the nucleons. 
This intrinsic k± can be observed in Drell-Yan experiments where it is found that the distribution is 
roughly independent of s and As inferred from these analyses, the width in (|l^) is /cq = 0.42 

GeV, corresponding to a mean {k±) 0.38 GeV. 


The scheme (i)-(iii) to sample the flavor and momentum distribution is carried out independently for 
each nucleon subject to the requirement 



(14) 


where the summation runs over all partons belonging to the same nucleon as determined by (1C) and (pi]), 
and M/j is the mother hadron (nucleon) mass. With the partons’ 4-momenta distributed as outlined above, 
the constraint (M) determines the distribution in the variable = k? — rn^ = — k'^ — < 0, the 

partons’ initial space-like virtualities. The resulting distribution in is a strongly peaked Gaussian with a 
mean value of « Q\l^- 


2.5.2 Initial spatial distribution 

The initial spatial distribution of the partons, R^\ appearing in eq. (j^, depends on the magnitude of their 
momenta, the positions of their parent hadrons (nucleons), as well as on the spatial substructure of the 
latter. It is represented as 


7?: 




(p,r,i7) = b^(R- Rab) \ h^^{r) Hn^R) 


boosted 


^This feature is optional in the simulation, by default it is switched off. 


(15) 


13 




















where the momentum dependence is here purely due to boosting the distributions to the global cm-frame. 

The components of the ansatz have the following meaning: 

(i) The incoming beam particles (a single hadron, or a nucleus) are assumed to have initial cm-positions 
Rab = Bab), corresponding to a chosen impact parameter Bab = {Bx,By) and a min¬ 

imum longitudinal separation AZab- In the case where A and/or B is not a single hadron, but a 
nucleus, the individual nucleons are assigned positions around the centers of the parent nucleus in its 
rest-frame, according to a Fermi distribution for nuclei with mass number A > 12 and a Gaussian 
distribution for nuclei A < 12, 


BnAR) 


^ (l + exp[(i2-c)/a]) ^ {A >12) 
^ ^ exp [-Ryb^] (A < 12) 


(16) 


The parameters are c = tq = 1.14 fm, a = 0.545 fm and b = \/2/3 , where i?™* is the 

mean square radius of the nucleus. 


(ii) Next, the partons are distributed around the centers of their mother hadrons or nucleons, still in the 
rest frame of the A, respectively B, with an exponential distribution 


hy{f) = — — exp [—j/r] 
“ ^ ^ 47r Svr ^ ^ ^ 


(17) 


where v = 0.84 GeV corresponds to the measured elastic formfactor of the mother hadron or nucleon, 
with a mean square radius of = \fV2jv = 0.81 /m. 

(iii) Finally, as indicated by the subscript boosted in eq. (|l^) , the positions of the hadrons or nucleons and 
their associated valence quarks are boosted into the global cm-frame of the colliding beam particles A 
and B. The valence quarks then occupy the Lorentz contracted region (Az)^ ~ 2Ra Mh/P, whereas 
the sea quarks and gluons are smeared out in the longitudinal direction by an amount {Az)g^s ~ 
1/kz < 2 Ra {Ra = r^A^/y behind the valence quarks. This is an important feature of the partons 
when boosting a hadron or nucleus to high rapidities “M, As a consequence, the parton 
positions are correlated in longitudinal direction with their momenta, as required by the uncertainty 
principle. This leads to an enhancement of the densities of gluons and sea quarks with x < 1 /{2RhMh) 
proportional to A^^^, because such partons from different nucleons overlap spatially when the nucleons 
are at the same impact parameter. 


2.6 Parton cascade development 

With the above construction of the initial state, the incoming beam particles A and B are decomposed in 
their associated parton content at the initial resolution scale Qq at time t = to = 0, and the subsequent 
dynamical development of the system for t > to can now be traced according to the kinetic equations (|^- 
(^). In the kinetic approach, the space-time evolution of parton densities may be described in terms of 
parton cascades. Fig. 4 illustrates a typical parton cascade sequence. It is important to realize that, in 
general, there can be many such cascade sequences, being internetted and simultanously evolving (typical 
for hadron-nucleus or nucleus-nucleus collisions). 

Each cascade can be subdivided into elementary 2^2 scatterings, 1 —> 2 branchings (emissions), and 
2^1 fusions (absorptions). In Fig. 4, a primary parton a that originates from one of the incident nuclei, 
collides with another parton b with some momentum transfer characterized by the scale at the collision 
vertex. The parton a has evolved from the initial scale Qq, at which it was resolved in its parent hadron or 
nucleus, up to by successive space-like branchings. From the scattering of a and b the partons c and 
d emerge, both of which can initiate sequences of time-like branchings. These newly produced partons can 
themselves branch, rescatter, or undergo fusion with other partons. 

In practice the dynamical interplay of 2 —> 2, 2 ^ 1 and 1^2 processes can be simulated as follows: 
Since the time evolution of the parton system is described in small discrete time steps At = O{10~^fm), 
the time scale of an interaction process is compared with At to decide whether the interaction occurs within 


14 









>- time 


> 


Figure 4; Schematical illustration of a typical parton cascade development initiated by a collision of two 
partons a and b. The incoming primary parton a evolves through a space-like cascade from the initial 
resolution scale Qq up to the scale at which it collides with parton b. The outgoing partons c and d 
both initiate time-like cascades which are described as a combination of multiple branchings (or emissions), 
fusions (or absorptions), or secondary scatterings (or rescatterings). A particular branch in a cascade tree 
terminates (locally in space-time), when the partons in that branch recombine with neighboring ones to 
pre-hadronic clusters and their subsequent decay into hadrons. 
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this time interval. A virtual parton e* with momentum k'^ ^ 0, that is produced via a + 5 —> e*, has 
a short life-time oc 1 / ^J\kP‘\, if its invariant mass || is large and it is likely to decay within yre* < At 
(the full formulae are derived in [26|). Therefore, the 2^2 process a + h ^ e* c + d preferably occurs. 
On the other hand, if \k‘^\ is small, corresponding to a long Tg*, parton e* is likely not to decay within the 
time step At and the fusion process a + b ^ e* rather happens. In this case, e* propagates as a quasi-stable 
particle until in the following time step it has an increased decay probability ||^. This may then result in 
the 1 —> 2 decay e* c + d. Alternatively, e* might collide with another parton before its decay or may 
propagate freely until the following time step, and so on. In this manner the elementary 2 —> 2, 2 —> 1, and 
1 —> 2 processes are treated on equal footing and double counting is prevented: either a 2 —> 2 , or a 2 —> 1 , 
and possibly a subsequent 1 —> 2 process can occur, but not both. The relative probability is determined by 
the uncertainty principle, i.e. by relating the momentum scale of the process to the time scale as explained. 

For the collisional processes 2 —> 2 and 2 —> I, the statistical occurence is determined by the 2-body 
cross-section in Born apprpoximation, with higher-order inelastic corrections effectively included in the 
parton evolution before and after each collision. This parton-evolution is calculated using the well-known 
based on the ‘modified leading logarithmic approximation’ (MLLA) to the QCD evolution 
Each individual parton-collision factorizes then in an elementary 2-body collision, 
and out-going partons undergo an ordered sequence of elementary branchings a ^ b + c, 


jet calculus [12 


T^|bai 
es |j,| 


of hard processes 
in which the in¬ 
accounting for higher-order radiative corrections to the lowest-order Born approximation. These branchingss 
can be described stochastically as a Markov cascade in position and momentum space. One distinguishes 
initial-state, space-like branchings of the two partons entering a collision vertex, and final-state, time-like 
radiation off the collided partons after the collision. 

The specific feature of the present approach is that, in addition to the definite virtuality and momentum, 
each elementary vertex has a certain space and time position which is obtained by assuming that the partons 
in the cascade propagate on straight-line trajectories in between their interactions. In the MLLA framework, 
the basic properties of both space-like and time-like showers are determined by the DGLAP equations 
[p^], but with essential differences in time ordering, kinematics and the treatment of infrared singularities 
associated with soft gluon emission. 


2.6.1 Elementary parton-parton collisions 

For the elementary parton scatterings a b ^ c + d, and fusion processes a + b ^ c*, two distinct classes 
of processes are considered: 

(i) truly perturbative QCD hard parton collisions with a sufficiently large momentum transfer or invari¬ 

ant mass s; 

(ii) phenomenological treatment of soft parton collisions ^ with low momentum transfer or invariant 

mass s. 


The motivation for this differentiation of parton-parton collisions is to regulate the singular behavior of 
the collision integrals in (^) which results from the divergence of the associated Born-amplitudes squared 
for small momentum transfers (except for s channel processes such as qq annihilation). To render the 
parton-parton cross-sections finite, an invariant hard-soft division scale Pq is introduced such that collisions 
occurring at a momentum scale > Pq are treated as perturbative (semi)hard collisions, whereas for those 
with < Po a soft, non-perturbative interaction is assumed to occur. That is, the total parton-parton 
cross-section for a collision between two partons a and b is represented as 


(^ab[s 



where 


•s = {Pa+Pbf 



fdQ- 

) 

dpi 

f pi - - 

-i{-u) 

( 

= s 

-Pcf , 

u 




\ dQ^ 


'^The inclusion of soft parton collisions is optional in the program, and by default is switched off. 


(18) 


(19) 

( 20 ) 
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That is, the scale of the collision is set by p\ for scatterings a + 6 —> a + 6 , or by s for annihilation 
processes and a + 6 —> c+d and fusion processes a + b ^ c. The sum over c, (d) corresponds to summing over 
all possible reaction channels (processes). The specific value of po is generally dependent on beam energy 
Ecm = {Pa + Pb)^ of colliding hadrons (nuclei) A + B, as well as on the mass numbers A and B. It is 
treated as a parameter that is determined by existing experimental data. A suitable form is 

PO = po{E,^,A,B) = I I^ ^WGeV y („ = 2.0 GeV, 6 = 0.27). (21) 


The complementation of both hard and soft contributions renders the parton cross-section finite and 
well defined for all Q^. It must be emphasized that the effect of soft parton collisions on the global dy¬ 
namics is not essential for collisions involving only leptons and/or hadrons, but plays an important role in 
hadron-nucleus, or nucleus-nucleus collions. The soft parton collisions naturally involve comparably small 
momentum transfer, so that their contribution to transverse energy production is small, but the effect on 
soft particle production is significant. 

In accord with the two-component structure of the cross-section 
guished depending on the momentum transfer; 

(i) Hard collisions above po: 

For the perturbative (semi)hard collisions above po, the momentum scale is determined by the 
corresponding differential cross-sections 


18), these parton collisions are distin- 


dQ^ 


1 

IGvrs^ 


A^ab- 


*cd 


(s,(5^ 


TialiQ'^) 
Q 2 ’ 


( 22 ) 


where |Mp is the process-dependent spin- and color-averaged squared matrix element ^Born approx¬ 
imation, published in the literature (see e.g. Refs. |^| and for massive quarks Refs. 


(ii) Soft collisions below po: 


In the case of a non-perturbative soft collision between two partons it is assumed that a very low 
energy double gluon exchange occurs. This provides a natural continuation to the harder collisions 
above po where the dominant one gluon exchange processes gg —> gg, gq —> gq and qq —> qq have the 
same overall structure [^. A simple, and physically plausible, form for the soft cross-section that 
continues the hard cross-section for below po down to = 0 , may be modelled by introducing a 


screening mass-term in the nominator of (^), 


7 soft 

d<^abAcd 

dQ^ 


oc 


2iTa‘l{p\ 




(23) 


where /x is a phenomenological parameter that governs the overall magnitude of the integrated oc 
ln[(pQ -|- /x^)//x^]. The value of /x is not known precisely, but it can be estimated to be in the range 0.3 
- 1 GeV. 


Notice that both hard and soft scatterings are treated in this approach on completely equal footing. 
That is, with the four momenta of the incoming partons known, the momentum transfer and scattering 
angle are sampled from the respective differential cross-sections da/dt (p^ and (|^), assuming azimuthal 
symmetry of the scattering geometry. 


2.6.2 Space-like parton evolution 

Space-like parton cascades are associated with QGD radiative corrections due to brems-strahlung emitted by 
primary space-like partons, which are contained in the initial state hadrons (nuclei), and which encounter 
their very first collision. Thereby such a primary, virtual parton is struck out of the coherent hadron 
(nucleus) wavefunction to materialize to a real excitation. A typical space-like cascade is illustrated in the 
left part of Fig. 4, where the parton oq, originating from the initial hadron (nucleus) state, undergoes 
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successive space-like branchings oq —> aia'i, ai —> 0202 , fln-i —> OnOn to become the parton a = a-n which 
then actually collides with another parton b. The branching chain proceeds by increasing the virtualities |g?| 
of the partons a* in the cascade, starting from uq with |( 7 q| ~ Q^/4 up to \q‘^\ ~ the space-like virtuality 
of the scattering parton a, where sets the scale of the collision between partons a and b, and hence 

for the evolution from Qq to Q^. The emitted partons on the side branches with momenta ki on the other 
hand have, due to energy-momentum conservation at the branching vertices, time-like virtualities and each 
of them can initiate a time-like cascade, discussed afterwards. 

The proper inclusion of both collinear, hard and coherent, soft branchings is achieved by describing the 
space-like cascade in both space-time and momentum space by using a so-called angular-ordered evolution 
variable (rather than the virtualities |g||) ||32|| 

Qj = Cj+i , 0+1 = - 1 - cos 00 , j+i (0 < j < n) , (24) 

c+o h/j+i 

where qj = {ojj,qj) and O+i = refer to the branching aj aj^iaj_^_^ with momentum 

assignment qj —> qj^ikj^i (see Fig. 4). The space-like cascade is then strictly ordered in the variable 
qj+i > q], which is equivalent to the ordering of emission angles, i+i' < j+ 2 '- 

The space-like cascade terminates with parton an ^ a entering the vertex of collision with parton b, that 
is, in Fig. 4. The history of parton a is however not known until after it has collided with parton 6 , 
because it is this very collision that causes the cascade evolution of parton a. Therefore one must reconstruct 
the cascade backwards in time starting from the time of the collisions at the vertex and trace the history 
of the struck parton a back to the initial state at time to = 0 a-t which it was originally resolved with Qq in its 
hadron (nucleon) mother. The method used here is a space-time generalization of the ‘backward evolution 
scheme’ [^]. To sketch the procedure, consider the space-like branching qn-i —> qnkn which is closest to the 
collision vertex in Fig. 4. The virtualities satisfy |g^| > and qn,qn-i < 0 (space-like) but 

kn > 0 (time-like). The relative probability for a branching to occur between q‘^ and -|- d(p is given by 


dVl%^n{Xn-l,Xn,f-, At) 


dq‘^ dz as ((1 — z)(p) 
(p z 27r 


In-l- 


'{z) 


( F{rn-l-,Xn-l,P) \ 
y F{rn]Xn,P) ) 


T(^)(At) ,(25) 


where xj = [qpz/Pz (j = n,n — 1) are the fractions of longitudinal momentum of the initial mother 
hadron (nucleon), with F{rj‘,Xj,p) = F{rj,qj) the corresponding parton distributions defined by (® , and 
the variables 
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^n-l 


^n—1 

^n—1 


(26) 


specify the fractional energy or longitudinal momentum of parton n and n', respectively, taken away from 
n — 1. The function as/(27r) jIz) is the usual DGLAP branching probability |§, |TJ], with j(z) giving 
the energy distribution in the variable The last factor in (25) determines the time interval 

in the global cm-frame. At = tn — tn-i, that is associated with the branching process a„_i —> 0 ^ 0 )^. It 
accounts for the formation time of from a^-i on the basis of the uncertainty principle: At = Auj/\q‘^\, 
Aiv ~ [xn — Xn-i) Pz- A Very simple form is taken here. 


T(^)(At) = <5 Pz - At) . (27) 

The backwards evolution of the space-like branching qn-i qn + kn is expressed in terms of the probability 
that parton o„_i did not branch between the lower bound ^q, given by the initial resolution scale Qq, and 
(p = p czL In that case, parton n can not originate from this branching, but must have been produced 
otherwise or already been present in the initial parton distributions. This non-branching probability is given 
by the Sudakov form-factor for space-like branchings: 


Sn{xn,pAh At) = exp 



dPi%-iiXn,Z, 



(28) 
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where the sum runs over the possible species a = g,q,q oi parton a„_i. The upper limit of the p^-integration 
is set by associated with the collision vertex with parton b. The limits z± are determined by 

kinematics |Q: Z-{q) = Qo/q and Z-^-{q) = I — Qo/q. 

The knowledge of the space-like formfactor Sn{xn-,q^A q'i is enough to trace the evolution of the 

branching closest to the hard vertex backwards from q'^ at t = tn, the time of collision in the global cm- 
frame, q^-i at tn-i = tn — Xn/\qn\Pz- The next preceding branchings qn -2 —> Qn-ikn-i, etc., are then 
reconstructed in exactly the same manner with the replacements tn tn-i, Xn Xn-i, qn Qn-i^ so 
forth, until the initial point Qq at to = 0 is reached. 


2.6.3 Time-like parton evolution 


Time-like parton cascades are initiated by secondary partons that emerge either from the side-branches of 
a preceding space-like or time-like cascade, or directly from a scattering or fusion process. Consider the 
time-like cascade initiated by the parton c in the right part of Fig. 4, with momentum k = kn- This parton 
has been produced in the collision a -|- 6 ^ c -|- d with a time-like off-shellness k'^ ~ Q^b- 

Again an angular-ordered (rather than virtuality-ordered) evolution in space-time and momentum-space 
of the cascade is employed to incorporate interference effects of soft gluons emitted along the time-like 
cascade tree of Fig. 4, c = Cm —> ci —> coCg. In contrast to (^), the time-like version of the 

angular evolution variable is 




Ci-i = 


kj—i 




i-i 


Ej-iEj-i 


1 — cos ( 




(m > j > 1) 


(29) 


so that the time-like cascade can be described by a fc^-ordered (rather than fe^-ordered) evolution, which 
corresponds to an angular ordering with decreasing emission angles 6jj: > 

FTZI (T) 

Proceeding analogously to the space-like case (c.f. (P5|)), the probability dV!^ [n-i the first branching 
after the 'yq vertex, km —> km-ik'm-i with k‘^_i,k'^_i, is given by the space-time extension [|I|, ^ of the 
usual DGLAP probability distribution i 0 ], 

(iPm^m-liz,k^-, ^t) = ^dz ^ 7m^(m-l),(m-iy(^) , (30) 

where T^E{At) is the probabilty that parton m with virtuality km and corresponding proper lifetime 
Tm oc l/\/^ decays within a time interval At, 

TV)(Ai) = 1 - exp (-j^) . (31) 


The actual lifetime of the decaying parton m in the global cm-frame is then tm{k) = ^/Tm{k), where 
tq{k) « 3E/{2ask‘^) for quarks and tg{k) ^ E/{2ask‘^) for gluons 0]. As before, Ej denotes the local density 
of parton species j = m, m — 1, and Ug/{2 'k)^{z) is the DGLAP branching kernel with energy distribution 
7 ( 2 ). The probability (|30|) is formulated in terms of the energy fractions carried by the daughter partons. 




Em—l 

Em 


1 - Z 


Em-1 


Er, 


(32) 


with the virtuality km of the parton m related to z and ^ through 2E^z{l — z)^, and 

the argument in the running coupling Ug in (|^) is kP' = 2 z^(l — z^E^^ ~ /c^. 

The branching probability (|^ ) determines the distribution of emitted partons in both coordinate and 
momentum space, because the knowledge of four-momentum and lifetime (or At between successive branch¬ 
ings) give the spatial positions of the partons, if they are assumed to propagate on straight paths between 


the vertices. The probability that parton m does not branch between and a minimum value fcg = /Tq is 
given by the exponentiation of (30), yielding the Sudakov form-factor for time-like branchings: 


Tmik'^,^; At) = exp 



At) 


/ 2 . 


(33) 
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which is summed over the species a = g, q, q oi parton m —1. The integration limits fcg and z± are determined 
by the requirement that the branching must terminate when the partons enter the non-perturbative regime 
and begin to hadronize. This condition can be parametrized by the confinement length scale Lc = 0(1 fm) 

with > L~'^ = |Ug, and z+{km) = 1 — Z-(km) = k-o/\J so ^+(^o) — — 1/2 the phase 

space for the branching vanishes. 

The time-like form factor Tmik"^, Kq; At) determines the four-momenta and positions of the partons of a 
particular emission vertex as sketched above for the first branching from k'^ at t = tm, the time of production 
of parton c in the global cm-frame, to at tm-i = tm + Em/\km\- Subsequent branchings are described 

completely analogously by replacing tm —> tm-i, ^ otc.. Hence T{k‘^, k^; At) generates a time-like 

cascade as sequential branchings starting from t = 0 at the hard vertex forward in time, until the partons 
eventually hadronize as discussed below. 


2.7 Cluster formation and hadronization 

In view of lack of knowledge about the details of confinement dynamics and the non-perturbative hadroniza¬ 
tion mechanism, one must rely at present on model-building. In the present approach, the cluster- 
hadronization scheme of Ref. is employed. This phenomenological scheme is inspired by the 

Marchesini-Webber model ||^, however it works in space-time plus color-space. On the other hand it is 
very different from commonly used string-fragmentation models such as the Lund-model [^]. 

In VNI, both the cluster-formation from the collection of quarks and gluons at the end of the perturbative 
phase and the subsequent cluster-decay into final hadrons consist of two components: 

1. The recombination of the secondary time-like partons, their conversion into colorless parton clusters and 

the subsequent decay into secondary hadrons. 

2. The recombination of the primary spaee-like partons that remained spectators throughout the collision 

development into beam clusters and the fragmentation of these clusters. 

The important assumption here is that the process of hadron formation depends only on the local space-, 
time-, and color-structure of the parton system, so that the hadronization mechanism can be modelled as 
the formation of color-singlet clusters of partons as independent entities (pre-hadrons), which subsequently 
decay into hadrons. This concept is reminiscent of the ‘pre-confinement’ property [^] of parton evolution, 
which is the tendency of the produced partons to arrange themselves in color-singlet clusters with limited 
extension in both position and momentum space, so that it is suggestive to suppose that these clusters are 
the basic units out of which hadrons form. 


2.7.1 Cluster formation 


(i) Parton clusters: 

Parton clusters are formed from secondary partons, i.e. those that have been produced by the hard 
interaction and the parton shower development. The coalescence of these secondary partons to color-neutral 
clusters has been discussed in detail in Refs. P, 1R|. Throughout the dynamically-evolving parton cascade 
development, every parton and its nearest spatial neighbour are considered as as potential candidates for 
a 2-parton cluster, which, if color neutral, plays the role of a ‘pre-confined’ excitation in the process of 
hadronization. Within each single time step, the probability for parton-cluster conversion is determined for 
each nearest-neighbor pair by the requirement that the total color charge of the two partons must give a 
composite color-singlet state (if necessary by accompanying gluon emission), and the condition that their 
relative spatial distance L exceeds the critical confinement length scale Lc- The scale L is defined as the 
Lorentz-invariant distance Ljj between parton i and its nearest neighbor j\ 


L — Lij — min(Aji,..., A^j, ■ ■ ■, Am) 


(34) 


where ^{nfil + {nfil + {nfil + {nfil, with 


= r^ — rj, and the probability for the coalescence of 
the two partons i, j to form a cluster is modelled by a distribution of the form 


n 




oc 


1 — exp {—AF Lij) ] ~ I — exp 


Lq — L. 




Ec Eij 


if Ln 


^ Eij ^ Lc ? 


(35) 
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and liij^c = 0 (1) if Lij < Lq {Lij > Lc). Here AH is the local change in the free energy of the system that 
is associated with the conversion of the partons to clusters, and the second expression on the right side is 
a parametrization in terms of Lq = 0.6 fm and = 0.8 fm that define the transition regime. As studied 
in Ref. |R], the aforementioned color constraint, that only colorless 2-parton configurations may produce 
a cluster, can be incorporated by allowing coalescence for any pair of color charges, as determined by the 
space-time separation and the probability (35), however, accompanied by the additional emission of a 
gluon or quark that carries away any unbalanced net color in the case that the two coalescing partons are 
not in a colorless configuration. 


(ii) Beam clusters: 

If one or both of the colliding beam/target particles A and B were a hadron or a nucleus, the one, re¬ 
spectively two beam clusters are formed from the spectator partons that represent the receding beam/target 
remnants of the original particles A and B. More precisely, the remaining fraction of the longitudinal mo¬ 
mentum and energy that has not materialized or been redirected and harnessed during the course of the 
collision, is carried by those primary partons of the initial-state hadrons or nuclei, which remained spectators 
throughout. In the present approach these partons maintain their originally assigned momenta and their 
space-like virtualities. Representing the beam remnants of A and/or H, they may be pictured as the co¬ 
herent relics of the original hadron (or nucleus) wavefunctions. Therefore the primary virtual partons must 
be treated differently than the secondary partons which are real excitations that contribute incoherently 
to the hadron yield. In the global cm-frame, the primary partons are grouped together to form a massive 
beam cluster with its four-momentum given by the sum of the parton momenta and its position given by 
the 3-vector mean of the partons’ positions. 


2.7.2 Hadronization of clusters 
(i) Parton clusters: 

For the decay of each parton cluster into final-state hadrons, the scheme presented in detail in Refs. Q 
is employed: If a cluster is too light to decay into a pair of hadrons, it is taken to represent the lightest 
single meson that corresponds to its partonic constituents. Otherwise, the cluster decays isotropically in its 
rest frame into a pair of hadrons, either mesons or baryons, whose combined quantum numbers correspond 
to its partonic constituents. The corresponding decay probability is chosen to be 


= %{Ec,rnl) N j — exp I-) , (36) 

Jrrih ^ \ ^ 0 / 

where M is a normalization factor, and the integrand is a Hagedorn spectrum p8|| that parametrizes quite 
well the density of accessible hadronic states below rric which are listed in the particle data tables, and 
mo = m^. In analogy to (ID, 7/ is a life-time factor giving the probability that a cluster of mass m^ decays 
within a time interval At in the global cm-frame. 


Ec{Ec,ml) 


1 — exp 


At \ 
tc(Hc,m2)j 


(37) 


with the Lorentz-boosted life time tc = 7 cT c — Ec/rn^. In this scheme, a particular cluster decay mode is 
obtained from ( pBD by summing over all possible decay channels, weighted with the appropriate spin, flavour, 
and phase-space factors, and then choosing the actual decay mode acording to the relative probabilities of 
the channels. 

(ii) Beam clusters: 

The fragmentation of the beam clusters containing the spectator partons mimics in the present model 
what is commonly termed the ‘soft underlying event’, namely, the emergence of those final-state hadrons 
that are associated with the non-perturbative physics which underlies the perturbatively-accessible dynamics 
of the hard interaction with parton shower fragmentation. 

In the spirit of the Marchesini-Webber model a version (suitably modified for the present purposes) 
of the soft hadron production model of the UA5 collaboration |^, which is based on a parametrization of 
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the CERN pp collider data for minimum-bias hadronic collisions. The parameters involved in this model 
are set to give a good agreement with those data. 

Soft hadron production is known to be a universal mechanism ||40|| that is common to all high-energy 
collisions that involve beam hadrons in the initial state, and that depends essentially on the total energy- 
momentum of the fragmenting final-state beam remnant. Accordingly, one may assume that the fragmenta¬ 
tion of the final-state beam cluster depends solely on its invariant mass M, and that it produces a charged- 
particle multiplicity with a binomial distribution p^. 


P{n) 


r(n-|-/c) (n/k)"' 

nW{k) JT+^Jk)^ ’ 


(38) 


where the mean charged multiplicity n = n{M‘^) and the parameter k = k{M‘^) depends on the invariant 
cluster mass | according to the following particle data parametrization 


n{M'^) = 10.68 (M 


2 \ 0.115 


- 9.5 


= 0.029 ln(M2) - 0.064 


(39) 


Adopting the scheme of Marchesini and Webber the fragmentation of a beam cluster of mass M 
proceeds then as follows: First, a particle multiplicity n is chosen from (^), and the actual charged particle 
multiplicity is taken to be n plus the modulus of the beam cluster charge. Next, the beam cluster is split 
into sub-clusters (( 71 ^ 2 ), ( 92 ^ 3 ), • • • {<li = u,d), which are subsequently hadronized in the beam cluster rest 
frame, in the same way as the parton clusters described above. To determine the sub-cluster momenta, the 
following mass distribution is assumed, 


P{M) = c (M - 1) exp [-a(M - 1)] , 


with c a normalization constant and a = 2 GeV resulting in an average value of (M) 
transverse momenta are taken from the distibution 


P{P±) = c p± exp 




(40) 

1.5 GeV. The 

(41) 


with normalization c' and slope parameter b = 3 GeV“^, and the rapidities y are drawn from a simple flat 
distribution P{y) oc const, with an extent of 0.6 units and Gaussian tails with 1 unit standard deviation 
at the ends. Finally, all hadronization products of the sub-clusters are boosted from the rest frame of the 
original beam cluster back into the global cm-frame. 


3 PROGRAM DESCRIPTION 

3.1 The package VNI-3.1 

The program package VNI-3.1 is a completely new written code, using only fragments of its precedessors 
VNI-1.0 (Ref. IQ) and VNI-2.0 (Ref. ||^]). The current program is a hrst release of a long term project 
that aims at unifying the simulation of a wide class of high energy QCD processes, ranging from elementary 
particle physics processes, e.g., e^e“-fragmentation, deep inelastic ep-scattering, via hadronic collisions such 
as pp (pp), to reactions involving (heavy) nuclei, e.g. deep inelastic eA scattering, pA, or AA (AB) collisions. 
The types of collision processes that are available and are discussed in more detail in Section 3.4 below, 
include collisions involving for leptons, p,n,7r, A,Ti... for hadrons, and a, He, O,..., Au, Pb, U for 

nuclei. Some of these collision processes have of course been investigated before; here they are revisited 
within the space-time picture of parton cascades and cluster dynamics in order to provide an alternative 
method for Monte Carlo simulation and also to check the consistency with the wide literature on both 
experimental measurements and theoretical predictions. 

® Notice that in this model M fluctuates statistically, as a result of fluctuations of the initial-state parton configuration in 
the incoming hadrons (or nuclei), as well as due to the fluctuating number of remnant partons during the space-time evolution. 
Hence the distribution ( ^ ) and the mean multiplicity ( pil) vary from event to event. This is in contrast to the original UA5 
model, in which the fixed beam energy ^/s/2 controls the energy dependence of soft hadron production. 
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3.2 Overall structure of the program 


The program package VNI is written entirely in Fortran 77, and should run on any machine with such a 
compiler. The program VNI adopts the common block structure, particle identification codes, etc., from the 
Monte Carlos of the Lund family, e.g., the JETSET/PYTHIA programs [M|. In addition, various program 
elements of the latter are used in the sense of a library, however modified to incorporate the additional 
aspects of the space-time and color label information. VNI also employs certain (significantly modified) 
parts of the HERWIG Monte Carlo [^], used for the final decays of clusters into hadronic resonances and 
stable particles. It is important to stress that the correspondence of VNI with JETSET/PYTHIA and 
HERWIG is rather loose and should not lead to naive identification, in particular since VNI is taylored for 
a larger set of degrees of freedom that is necessary to describe the space-time dynamics and the information 
on color correlations. Nevertheless, by truncating the additional information of the particle record of VNI, 
a direct interface to the particle records of the Lund Monte Carlo JETSET/PYTHIA, and to the program 
HERWIG, is actually straightforward. Instructions are given below in Section 3.10. 

The package VNI-3.1 consists of the following files: 


(i) the documentation vni-3.1.ps (this one) 

(ii) the program VNI-3.1 in 7 parts: vnil.f — vni7.f 

(iii) two include files vnil.inc and vni2.inc 

(iv) two example programs vnixpll.f and vnixpl2.f . 
(vi) a simple, portable histogram package vnibook.f . 


The seven components vnil.f — vni7.f of Fortran source code contain all the subroutines, functions, 
and block data, as categorized above. All default settings and particle data are automatically loaded by 
including the two include files vnil.inc and vni2.inc. The example programs vnixpll.f {a generic one) and 
vnixpl2.f {a more application-oriented one), as well as the histogram package vnibook.f are not integral part 
of the program, i.e. are distributed supplementary as useful guiding tools. 

The program source code is organized in 7 parts, in which subroutines and functions that are related in 
their performance duties, are grouped together (a summary list of all components with a brief description 
of their purpose is given in Appendix E): 

vnil.f: Main steering routines VNIXRIN (initialization), VNIXRUN (event generation), and VNIXFIN 
(finalization), plus other general-duty subroutines. VNIXRIN and VNIXFIN are only to be called 
once, while VNIXRUN is to be called anew for each event, in which the particle system is evolved in 
discrete time steps and 7-dimensional phase-space. 

vni2.f: Initialization routines. They include, e.g. the process dependent event-by-event initialization of the 
chosen type of collision process, and various other slave routines associated with the initial set up of 
an event. 

vniS.f: Package for parton distributions (structure functions). A portable and self-contained collection of 
routines with a number of structure function parametrizations, plus a possible interface to the CERN 
PDFLIB, plus various slave routines. 

vni4.f: Evolution routines. These form the heart of the parton cascade/ cluster-hadronization simulation, 
and include both (process-dependent) particle evolution routines and (process-independent) universal 
routines for space-like and time-like shower evolution, for parton scatterings, parton recombination to 
clusters, and cluster decays to hadrons and other final state particles. 

vni5.f: Diverse utility routines and functions performing the slave-work for the part of perturbative parton 
evolution. 

vniG.f: Diverse utility routines and functions performing the slave-work for the part of cluster fragmentation 
to hadrons. 
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vni7.f: Event study and analysis routines for accumulating, analysing and printing out statistics on observ¬ 
ables. 


3.3 Special features and machine dependence 

Important to remark are the following special features, being not strictly common in standard Fortran 77: 

1. The program uses of a machine-dependent timing routine, called VNICLOC, which measures the elapsed 

CPU time 1^. Two alternatives are provided: a) the TIMEX routine of the CERN library, in which 
case the latter must be linked with the program, or, b) the MCLOCK routine which is specific 
to IBM-AIX compilers. If neither of these options are desired, or does not work in the users local 
environment, the lines involving calls for TIMEX {MCLOCK) calls can easily commented out. [|. 

2. All common block variables and dimensions are defined globally within the include files vnil.inc and 

vni2.inc. This greatly simplifies for instance possible modifications (e.g., enlarging or decreasing) of 
global common block arrays. Consequently, all subroutines contain INCLUDE statements for these 
include files, and therefore when compiling and linking the program, the files vnil.inc and vni2.inc 
must be properly accessible. 

3. Single precision is normally assumed throughout, since most machines are 32 bit processors. For applica¬ 

tions at extremely high energies, single precision for any real variable starts to become a problem. The 
main source of numerical precision-loss arises from multiple Lorentz boosts of particles throughout the 
simulation. Therefore the relevant parts of the code are written in double precision. All double pre¬ 
cision variables have as first character D. As additional precaution, particle energies are recalculated 
from momentum and mass after each time step during the evolution of an event. 

4. Implicit variable-type assignment is assumed, such that variables or functions beginning with letters I - 

M are of integer*^ type, and variables or functions with first letters A - H and O - Z are of type reafrd 
(single precision), except for occasional reahS (double precision) occurrences as stated in item 3. 

5. SAVE statements have been included in accordance with the Fortran standard. Since most ordinary 

machines take SAVE for granted, this part is not particularly well tried out, however. Users on 
machines without automatic SAVE are therefore warned to be on the lookout for any variables which 
may have been missed. 


3.4 The main subroutines 

There is a minimum of three routines that a user must know: 

1. VNIXRIN for the overall initialization of a sample of collision events, 

2. VNIXRUN for the subsequent actual generation of each event, 

3. VNIXFIN for finishing up simulation. 

These three routines, which are briefly descibed below, are to be called in that order by a user’s program as 
is exemplified by the example programs vnixpll.f and vnixpl2.f (see also Section 3.10). 

SUBROUTINE VNIXRIN(NEVT,TFIN,FRAME,BEAM,TARGET,WIN) 

Purpose: to initialize the overall simulation procedure. 


® The usage of the time measurement is not necessary for the performance of the program as-a-whole, although it is a nice 
convenience. 

^ The subroutine VNICLOC is listed as the very first one in vnil.f. 
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NEVT : integer specifying the number of collision events to be simulated for a selected physics process. 

TFIN : final time (in fm) in the global Lorentz frame up to which each collision event is followed. 

FRAME : a character variable used to specify the global Lorentz frame of the experiment. Uppercase 
and lowercase letters may be freely mixed. 

= ’CMS’ : colliding beam experiment in CM frame, with beam momentum in +z direction and target 
momentum in —2: direction. 

= ’FIXT’: fixed target experiment, with beam particle momentum pointing in +z direction. 

= ’USER’ : full freedom to specify frame by giving beam momentum in PSYS(1,1), PSYS(1,2) 
and PSYS(1,3) and target momentum in PSYS(2,1), PSYS(2,2) and PSYS(2,3) in common 
block VNIREC. Particles are assumed on the mass shell, and energies are calculated 
accordingly. 

= ’FOUR’ : as ’USER’, except also energies should be specified, in PSYS(1,4) and 
PSYS(2,4), respectively. The particles need not be on the mass shell; effective 
masses are calculated from energy and momentum. (But note that numerical precision 
may suffer; if you know the masses the option ’FIVE’ below is preferrable.) 

= ’FIVE’ : as ’USER’, except also energies and masses should be specified, i.e the 
full momentum information in PSYS(1,1) - PSYS(1,5) and PSYS(2,1) - PSYS(2,5) should 
be given for beam and target, respectively. Particles need not be on the mass 
shell. Space-like virtualities should be stored as . Four-momentum and mass 

information must match. 

BEAM, TARGET : character variables to specify beam and target particles. 

Uppercase and lowercase letters may be freely mixed. An antiparticle may be denoted 
either by (tilde) or (bar) at the end of the name. It is also possible to 
leave out the charge for neutron and proton. 

= ’e—’: electron e~ . 

= ’e-|-’ : positron e"*". 

= ’mu—’ : muon . 

= ’mu+’ : antimuon . 

= ’gamma’: real photon 7 (not yet implemented) . 

= ’pi+’ : positive pion vr^. 

= ’pi—’ : negative pion 11“ . 

= ’piO’ : neutral pion vr^. 

= ’nO’ : neutron n. 

= ’riO’ : antineutron n. 

= ’p-|-’ : proton p. 

= ’jU’ : antiproton p. 

= ’LambdaO’: A baryon. 

= ’Sigma-|-’ : S"*" baryon. 

= ’Sigma—’ : S“ baryon. 

= ’SigmaO’ : baryon. 

= ’Xi—’ : E~ baryon. 

= ’XiO’ : baryon. 

= ’Omega—’: baryon. 

= ’Alpha’ : fa particle. 

= ’He’ : oHe (Helium) nucleus. 

= ’Ox’ : g O (Oxygen) nucleus. 

= ’Su’ : fgS (Sulfur) nucleus. 

= ’Ag’ : 47® Ag (Silver) nucleus. 

= ’ Au’ : 79'^Au (Gold) nucleus. 

= ’Pb’ : ifspb (Lead) nucleus. 

= ’Ur’: if U (Uranium) nucleus. 

WIN : related to energy of system (in GeV) . Exact meaning depends on FRAME, i.e. for 
FRAME=’CMS’ it is the total 
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energy ^/s of system, and for FRAME=’FIXT' it is the absolute 3-momentum P = of 

beam particle. 


SUBROUTINE VNIXRUN(IEVT,TRUN) 

Purpose: to generate one event of the type specified by the VNIXRIN. This is the main routine, which 
administers the overall run of the event generation and calls a number of other routines for specific tasks. 

lEVT : integer labeling the current event number. 

TRUN : time (in fm) TRUN < TFIN in the global Lorentz frame up to which the space time evolution 
of the current event lEVT is carried out. If TRUN < TFIN, then the particle record contains the 
history up to this time, and repeated calls resume from TRUN of the previous call (c.f. Section 3.10). 


SUBROUTINE VNIXFINO 

Purpose: to finish up overall simulation procedure, to write data analysis results to files, evaluate CPU time, 
close files, etc.. This routine must be called after the last event has been finished (c.f. Section 3.10). 


For a deeper understanding of the physics routines and their connecting structure, I refer to Appendix 
A for a brief description. Generally, for each of the included collision processes, there is an initialization 
routine that sets up the initial state, and an evolution routine that carries out the time evolution of the 
particle distributions in phase-space, starting from the initial state. Within the evolution routines, a number 
of universal slave routines (i.e. independent of the process under consideration) perform the perturbative 
parton evolution in terms of space-like branchings, time-like branchings, and parton collisions, as well as the 
cluster hadronization. 

Finally, the program provides some useful pre-programmed event-analysis routines which collect, analyze 
and print out information on the result of a simulation, and give the user immediate access to calculated 
particle data, spectra, and other observables. A more detailed description of their duties, and how to call 
them can be found in Appendix B. 
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Table 1: Available collision processes. 


IPRO 

type of collision process 

1 

^IZq 2 — jets hadrons 

2 

Ul~ Zq ^ W~^W~ 4 — jets hadrons 

3 

'j + h ^ jets hadrons 

4 

7 -|- A ^ jets hadrons 

5 

1 + h ^ 1 + jets + X ^ hadrons 

6 

1 + A ^ 1 + jets + X ^ hadrons 

7 

h + 1 ^ 1 + jets + X ^ hadrons 

8 

h + h' ^ jets + X ^ hadrons 

9 

h + A ^ jets + X ^ hadrons 

10 

A + 1 ^ 1 + jets + X ^ hadrons 

11 

A4- h ^ jets + X ^ hadrons 

12 

A + A' ^ jets + X ^ hadrons 


3.5 The physics processes 

The program is structured to incorporate different classes of particle collision processes in a modular manner. 
The general classes that are part of the program are summarized in Table 1. All of the generic collision 
processes in Table 1 allow a further specification of beam and/or target particle, .e.g., I = e,^ for leptons, 
h = p,n,7r,... for hadrons, or A = a, He, O,... for nuclei. They are summarized in Table 2. 

As explained before, the actual evolution of any beam/target-particle collision is simulated on microscopic 
level of the system of partons, clusters and hadrons, and their interactions. For the parton interaction 
processes, a selection of elementary hard/soft 2^2 scatterings and 2^1 fusions are included in VNI, 
and future extensions are planned. In addition there are the standard space-like and time-like 1^2 
branching processes, which, when combined with the elementary tree processes, provide the usual parton 
cascade method of including higher order, real and virtual, corrections to the elementary tree processes. The 
parton-cluster formation and hadronization scheme includes a number of 2-parton recombination processes 
and the decay of the formed clusters into final state hadron. 

It is possible to select a combination of partonic subprocesses to simulate. For this purpose, all subpro¬ 
cesses are numbered according to an I SUB code. The list of allowed codes is given below. In the following g 


Table 2: Available beam and target particles. 
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Table 3; Elementary partonic subprocesses. 


Class 

ISUB 

type of subprocess 

a) 2 - 

-^> 2 : 

1 

QiQj - 

qm 



2 

Qidi - 

qkdk 



3 

QiQi 

99 



4 

QiQi 

91 



5 

qm 

77 



6 

did - 

qi9 



7 

qi9 - 

qn 



8 

99 - 

qkqk 



9 

99 - 

99 



10 

soft scattering 

b) 2 - 

1 : 

11 

qm 

^9* 



12 

qtg 

q* 



13 

99 - 

-9* 

c) 1 - 

-> 2 (space-like) : 

21 

9*- 

> qm 



22 

q* - 

^ qi9 



23 

9*- 

99 

d) 1- 

-> 2 (time-like) : 

31 

9 

qm 



32 

qi - 

^ qi9 



33 

9 - 

^ 99 


denotes a gluon, qi represents a quark of flavour i = 1,... ,nf, i.e. for nj = 6 this translates to d, u, s, c, b, t. 
A corresponding antiquark is denoted qi. The notation 7 is for a real photon, i.e. on shell. An asterix * 
denotes an off-shell parton. 


3.6 The particle record 

Each new event generated is in its entity stored in the commonblock VNIREC, which thus forms the event 
record. Here each particle that appears at some stage of the time evolution of the system, will occupy 
one line in the matrices. The different components of this line will tell which particle it is, from where 
it originates, its present status (fragmented/decayed or not), its momentum, energy and mass, and the 
space-time position of its production vertex. The structure of the particle record VNIREC follows closely 
the one of the Lund Monte carlos [^, |^, employing the same overall particle classification of particle 
identification, status code, etc., however with important differences concerning color- and space-time degrees 
of freedom.. 

Note: When in the following reference is made to certain switches and parameters MSTV or PARV, 
these are described in Sec. 3.7 below. 

PARAMETER (NV=100000) 

CDMMDN/VNIREC/N,K(NV,5),L(NV,2),P(NV,5),R(NV,5),V(NV,5) 

Purpose: Contains the event record, i.e. the complete list of all partons and particles in the current event. 


N : number of lines in the K, L, P, R, and V arrays occupied by the current event. N is continuously 
updated as the definition of the original configuration and the treatment of parton cascading, cluster 
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fragmentation and hadron production proceed. In the following, the individual parton/particle num¬ 
ber, running between 1 and N, is called I. The maximum N is limited by the dimension NV of the 
arrays. 

K{I,1) : status code (KS), which gives the current status of the parton/particle stored in the line. The 
ground rule is that codes 1-10 correspond to currently existing partons/particles, while larger codes 
contain partons/particles which no longer exist. 

= 0 ; empty line. 

= 1 : an undecayed particle or an unfragmented parton. 

= 2 : an unfragmented parton, which is followed by more partons in the same color-singlet parton 
subsystem. 

= 3 : an unfragmented parton with special color-flow information stored in K{1, 4) and K(I, 5), such 
that color connected partons need not follow after each other in the event record. 

= 6 : a final state hadron resulting from decay of clusters that formed from materialized (interacted) 
partons. 

= 7 : a final state hadron resulting from soft cluster decay of beam/target remnants formed from 
left-over (non-interacted) initial state partons. 

= 11 : a decayed particle or a fragmented parton, c.f. = 1. 

= 12 : a fragmented jet, which is followed by more partons in the same color-singlet parton subsystem, 
c.f. = 2. 

= 13 : a parton which has been removed when special color-flow has been used to rearrange a parton 
subsystem, cf. = 3. 

= 16 : an unstable and removed hadron resulting from decay of clusters that formed from materialized 
(interacted) partons, c.f. = 6. 

= 17 ; an unstable and removed hadron resulting from soft cluster decay of beam/target remnants 
formed from left-over (non-interacted) initial state partons, c.f. = 7. 

< 0 : Particle entries with negative status codes (e.g. -1, -6, -7) parallel in their meaning those from 
above, but refer to particles which are only virtually present and become active only after a certain 
formation time, upon which their status code is set to the appropriate positive value. 


K{I,2) : Flavor code (KF) for partons, hadrons and electromagnetic particles included in the current 
version of VNI. The particle code is summarized in Tables 3-9. It is based on the flavor and spin 
classification of particles, following the 1988 Particle Data Group numbering conventions. Nuclei are 
attributed a non-standard code that is internal to VNI, given by KF = 1000 • Nmutrons + Kprotons- 
A negative KF code, where existing, always corresponds to the antiparticle of the one listed in the 
Tables 3-9. 

K{I,3) : line number of parent particle or jet, where known, else 0. Note that the assignment of a particle 
to a given jet of a jet system is unphysical, and what is given there is only related to the way the event 
was generated. 

K{I,4) : Special color-flow information (for internal use only) of the form K(I,4) = M.ST\{2)-ICFR + 
ICTO, where ICFR and ICTO give the line numbers of the partons from which the color comes and 
to where it goes, respectively. 

K(/,5) : Special color-flow information (for internal use only) of the form K{I,5) = MSTY(2)-JCFR + 
JCTO, where JCFR and JCTO give the line numbers of the partons from which the anticolor comes 
and to where it goes, respectively. 


: color label {L = 1,... ,NC) of a parton, where NC \s number of colors specified by the parameter 
MSTV(5). Is = 0 for antiquarks and all non-colored particles. 

L{I,2) : anticolor label (L = 1,..., NC) of a parton. Is = 0 for quarks and all non-colored particles. 


: Px, momentum in the x direction in GeV. 
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P(/,2) : Py, momentum in the y direction in GeV. 

P(/,3) : Pz, momentum in the z direction in GeV. 

P(/,4) : E, energy, in GeV. 

P(/,5) : m = p^, mass in GeV. For partons with space-like virtualities, i.e. where = —m? > 0, its value 
is P{I,5) = = -Q. 

: current x position in frame of reference in 1/GeV. 

P(/,2) ; current y position in frame of reference in 1/GeV. 

P(/,3) : current 2 : position in frame of reference in 1/GeV. 

R{I,4) : time (in 1/GeV) of active presence in the system, i.e. tpres = t — tprod- Is set equal to PARW(2) 
at time of production tprod of the particle and then increased by PARW(2) in each timestep until the 
particle fragments or decays. If ( 0, giving the remaining formation time for virtually produced, but 
not yet real particles. 

P(/,5) : number and type of interactions of a particle, and is equal to MSTV(3)^-VCO-I- MSTV {3)- N S B + 
NTB, where - for partons - NCO is the number of 2-body collisions, NSB of space-like branchings, 
and NTB of time-like branchings, undergone up to current time. For clusters, hadrons and unstable 
particles, NSB is zero and NTB counts the number of 2-body decays. 

: X position of production vertex, in 1/GeV. 

V(/,2) : y position of production vertex, in 1/GeV. 

V(/,3) : 2 position of production vertex, in 1/GeV. 

V(/,4) ; time t of production, in 1/GeV. 

F(/,5) ; encodes origin of production as 10-MSTV(4)-/MO -|- 10 • ITM + OR, where IMO is the direct 
mother, ITM the time of production, and OR the generation of particle IP in a cascade tree, i.e. the 
genetical origin with respect to the production vertex of its original ‘Ur’-mother. (This keeps track of 
the production history even if the decayed mother has been removed and allows to follow the genealogy 
of the system.) 

Special cases: 

For the simulation of collisions involving hadrons or nuclei in the initial state, additional initial state 

information on the status and specific variables of ” primary partons” of the incoming hadron (nucleus) are 

stored as follows. Note however, that after the first interaction of a ’’primary parton”, this initial state 

information is lost and the assignments for the components of the arrays K, R, V are the same as above. 

K{I,4) : KF flavour code of nucleus from which a primary parton originates. 

K{I,5) : KF flavour code of hadron (nucleon), which is mother of the primary parton. (For hadrons 
K{I,4) and K{I,5) are equal). 

Ril,5) : Is equal to — |x|, where x is the longitudinal momentum fraction of primary parton that it carries 
off its mother hadron. 

V(I,4) : encodes location IB in the particle record of a primary sea-quark’s brother, i.e. the sea-antiquark 
belonging to the same vacuum polarization loop, as MSTV(4)^ •/R-|-MSTV(4)-ITM + OR, analogous 
to V(I, 5) above. For primary valence quarks and gluons it is 0. 
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Additional remark: The arrays K{I,3) — K{I,5) and the arrays P,R, and V may temporarily take special 
meaning other than the above, for some specific internal use. 


In Section 3.10 below, a typical event listing of the particle record is printed, which serves as an example 
for the organization of the particle record, exhibiting the information for momentum space: It lists the 
particles contained in VNIREC at a certain point during the evolution of an event, where KS and KF 
give status and flavor codes, C and A the color and anticolor index, P the energy-momentum-mass in GeV 
for each particle that appeared at some point in the event history. An analogous listing for the position 
space particle record (not printed here) contains the spatial coordinates and time of active presence for each 
particle, as well as the particle’s rapidity and transverse momentum with respect to the z-axis (jet-axis, or 
beam-axis). Note: both listings can be obtained at any point by calling the routine VNILIST, described in 
Appendix B. 


Table 4; Quark, lepton and gauge boson codes. 


KF 

Name 

Printed 

KF 

Name 

Printed 

KF 

Name 

Printed 

1 

d 

d 

11 

e“ 

e- 

21 

g 

g 

2 

u 

u 

12 

Ve 

nu_e 

22 

7 

gamma 

3 

s 

s 

13 

ti'~ 

mu- 

23 

ZO 

zo 

4 

c 

c 

14 


nu_mu 

24 

W+ 

w+ 

5 

b 

b 

15 

T 

tau- 

25 

HO 

HO 

6 

t 

t 

16 


nu_tau 

26 



7 



17 



27 



8 



18 



28 



9 



19 



29 



10 



20 
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Table 5; Special codes. 


KF 

Printed 

Meaning 

90 

CMsoft 

Center-of-mass of beam/target remnant system before soft fragmentation 

91 

cluster 

Cluster from coalescence of secondary partons 

92 

Beam-REM 

Cluster from remnant primary partons of beam particle 

93 

Targ-REM 

Cluster from remnant primary partons of target particle 

94 

CMshower 

Four-momentum of time-like showering system 

95 

SPHEaxis 

Fvent axis found with VNISPHE 

96 

THRUaxis 

Fvent axis found with VNITHRU 

97 

CLUSj et 

Jet (cluster) found with VNICLUS 

98 

99 

100 

CELLjet 

Jet (cluster) found with VNICELL 
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Table 6: Meson codes, part 1. 


KF 

Name 

Printed 

KF 

Name 

Printed 



pi+ 

213 


rho+ 


K° 


313 

K*° 

K*0 


K+ 


323 

K*+ 

K*+ 


D+ 


413 

D*+ 

D*+ 


D° 


423 

D*o 

D*0 

431 

D+ 

D_s+ 

433 

d:+ 

D*_s+ 

511 

B° 

BO 

513 

B*o 

B*0 

521 

B+ 

B+ 

523 

B*+ 

B*+ 

531 

B° 

B_s0 

533 

B*° 

B*_s0 

541 

Bo+ 

B_c+ 

543 

B,*+ 

B*_c+ 

111 

7r° 

piO 

113 


rhoO 

221 

V 

eta 

223 

UJ 

omega 

331 

v' 

eta’ 

333 

0 

phi 

441 

Vc 

eta_c 

443 

j/v- 

J/psi 

551 

Vh 

eta_b 

553 

T 

Upsilon 

661 

Vt 

eta_t 

663 

e 

Theta 

130 


K_L0 




310 

Kg 

K_S0 
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Table 7: Meson codes, part 2. 


KF 

Name 

Printed 

KF 

Name 

Printed 

10213 

bi 

b_l+ 

10211 

a^ 

a_0+ 

10313 

K? 

K_10 

10311 

KS° 

K*_00 

10323 

Ki- 

K_l+ 

10321 

k;+ 

K*_0+ 

10413 


D_l+ 

10411 

DS+ 

D*_0+ 

10423 

D? 

D_10 

10421 

DS° 

D*_00 

10433 

Di 

D_ls+ 

10431 


D*_0s+ 

10113 

b? 

b_10 

10111 

ao 

a_00 

10223 

h? 

h_10 

10221 

cO 

f_00 

10333 

hf 

h’_10 

10331 

f/0 

o 

o 

1 

10443 

h?c 

h_lcO 

10441 

XOc 

chi_0c0 

20213 

a+ 

a_l+ 

215 


a_2+ 

20313 

Kf 

K*_10 

315 

Kf 

K*_20 

20323 

Ki+ 

K*_l+ 

325 

k ;+ 

K*_2+ 

20413 

D*+ 

D*_l+ 

415 


D*_2+ 

20423 

D.o 

D*_10 

425 

Df 

D*_20 

20433 

DL+ 

D*_ls+ 

435 


D*_2s+ 

20113 

a? 

a_10 

115 

a2 

a_20 

20223 

fO 

f_10 

225 

fO 

f_20 

20333 

If 

f’_10 

335 

r/0 

l2 

f’_20 

20443 

Xlc 

chi_lcO 

445 

X2c 

chi_2c0 

30443 


psi’ 




30553 

T' 

Upsilon’ 





3.7 The general input and control parameters 


The VNIPAR common block contains the status code and parameters that regulate the performance of the 
program. All of them are provided with sensible default values, so that a novice user can neglect them, and 
only gradually explore the full range of possibilities. 

CDMMDN/VNIPAR/MSTV(200),PARV(200),MSTW(200),PARW(200) 

Purpose: to give access to status code and parameters that regulate the performance of the program. If the 
default values, denoted below by (D=...), are not satisfactory, they must in general be changed before the 
VNIXRIN call. Exceptions, i.e. variables that can be changed for each new event, are denoted by (C). 


MSTV(200), PARV(200): control switches and physics parameters 


1. General: 


MSTV(l) : (D=100000) number of lines available in the commonblock VNIREC. Should always be 
changed if the dimensions of the K, P, R, V arrays are changed. 

MSTV(2) : (D=25000) is used to trace internal color flow information in the array K. 
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Table 8; Baryon codes. 


KF 

Name 

Printed 

KF 

Name 

Printed 




1114 

A- 

Delta- 

2112 

n 

nO 

2114 

AO 

DeltaO 

2212 

p 

P+ 

2214 

A+ 

Delta+ 




2224 

A++ 

Delta++ 

3112 

S- 

Sigma- 

3114 

s*- 

Sigma*- 

3122 

AO 

LambdaO 




3212 

SO 

SigmaO 

3214 

S*o 

Sigma*0 

3222 

S+ 

Sigma+ 

3224 


Sigma*+ 

3312 

'w — 

Xi- 

3314 

'w * — 

Xi*- 

3322 

HO 

XiO 

3324 


Xi*0 




3334 

Q- 

Omega- 

4112 


Sigma_cO 

4114 


Sigma*_c0 

4122 

Ac+ 

Lambda_c+ 




4212 


Sigma_c+ 

4214 


Sigma*_c+ 

4222 

S++ 

Sigma_c++ 

4224 


Sigma*_c++ 

4132 

'='0 

‘c 

Xi_cO 




4312 

'c 

X! 

H- 

1 

O 

O 

4314 

‘c 

Xi*_c0 

4232 

‘c 

Xi_c+ 




4322 

'~'c 

Xi'_c+ 

4324 

‘~‘c 

Xi*_c+ 

4332 

OO 

0mega_c0 

4334 

11*0 

Dmega*_c0 

5112 

Sb 

Sigma_b- 

5114 

s- 

Sigma*_b- 

5122 

AS 

Lambda_bO 




5212 

sS 

Sigma_bO 

5214 

v*o 

^b 

Sigma*_b0 

5222 

Sb 

Sigma b+ 

5224 


Sigma*b+ 


Table 9; Nucleus codes. 


KF 

Name 

Printed 

KF 

Name 

Printed 

1001 

2002 

8008 

16016 

"fa 

|He 

lOQ 

32 Q 
16^ 

alpha 

he 

ox 

su 

61047 

118079 

126082 

146092 

47® Ag 

ifAu 

ifPb 

238tt 

92 U 

ag 

au 

pb 

ur 
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MSTV(3) ; (D=10000) is used to store number and type of interactions of a particle in the array R of 
VNIREC. 

MSTV(4) : (D=1000) is used to trace origin of a particle’s production in the array V of VNIREC. 
MSTV(5) : (D=3) number of colors Nc included in the simulation. 

MSTV(6) : (D=6) number of quark flavors Nj included. 

MSTV(7) : (D=2) energy and momentum conservation options. Since both the initial generation of 
initial hadronic or nuclear parton system and the final cluster-hadron decay scheme require occasional 
reshuffling of 4-momentum, an energy and/or momentum imbalance can occur. This can be corrected, 
if desired, by renormalizing particle energies and momenta. 

= 0 ; no explicit conservation of any kind. 

= 1 ; particles share energy imbalance compensation according to their energy (roughly equivalent 
to boosting event to CM frame). 

= 2 : as =1, plus particles share 3-momentum imbalance compensation according to their longitudinal 
mass with respect to the imbalance direction. 

MSTV (8) : (D=0) type of particles to be included in energy and momentum conservation options under 
MSTV(7). 

= 0 ; all living particles are included. 

= 1 ; only partons. 

= 2 : only hadrons. 

MSTV(9) : - not used - 

MSTV(IO) : (D=llllll) initial seed for random number generation. 

MSTV(ll) : (D=250000) maximum number of collision events A + B, after which program is forcibly 
terminated. 

MSTV(12) : (D=2500) maximum number of time steps per collision event, after which a new collision 
event is generated. When increased, the dimension of the arrays in commonblock VNITIM need to 
be altered accordingly. 

MSTV(13) : (D=0) choice of time grid TIME{I) and time increase TINC{I). Generically, t{i) = 
t{i - 1) + [f{i) - f{i - 1)], with f = 0,.. .,MSTV{12), and t(0) = to =PARV(12), tf =PARV(13). 

= 0 : power law increase f{i) = a ■ (i)^, where a =PARV(14) and b =PARV(15). 

= 1 ; logarithmic increase f{i) = a ■ ln(z)^, where a =PARV(14), b =PARV(15). 

MSTV(14) : (D=0) Option to select or switch on/off certain physics processes I SUB in the simulation, 
by specifying the array MSUB (c.f. Sec. 3.8 below), where MSUB{ISUB) = 0 (1) switches a 
subprocess off (on). 

= 0 : the relevant processes are initialized automatically. 

= 1 : user selection of included processes is required. 


PARV(12) : (D=0. fm) choice of initial point of time ti = PARV(12) when evolution of collision event 
begins. 

PARV(13) : (D=0. fm) choice of final point of time tf = PARV(13). If = 0., then it is set automatically 
to T/M£'(MSTW(2)), where MSTW(2) is the selected number of time steps. Otherwise, if > 0., it 
specifies tf such that the minimum mm(t f, tmax) is taken as the point when the evolution of a collision 
event stops, where tmax = TIME{M.STY{12)). 

PARV(14) : (D=0.05) prefactor a of the time-increment function /, see MSTV(13). 
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PARV(15) : (D=l.) exponent b of the time-increment function /, see MSTV(13). For longer collision 
time scales as in hadron-nucleus or nucleus nucleus collisions, it is automatically scaled by PARV(15) 
to 1.5 X PARV(15) to extend the time range. 


2. Initial state of collision system: 


MSTV(20) : (D=0) choice of model for initial parton distributions. 

= 0 : standard structure functions are used. 

= 1 : other parametrization, e.g. a la McLerran et al. (not implemented yet). 

MSTV(21) : (D=l) choice of nucleon structure-function set (c.f. MSTV(22)). 

= 1 : GRV94LO (Lowest order fit). 

= 2 : GRV94HO (Higher order; MSbar fit). 

= 3 : GRV94DI (Higher order; DIS fit). 

= 4 : - not used - 

= 5 : GTEQ2M (best higher order MSbar fit). 

= 6 : GTEQ2MS (singular at small x). 

= 7 : GTEQ2MF (flat at small x). 

= 8 : GTEQ2ML (large A). 

= 9 : GTEQ2L (best lowest order ht). 

= 10 : GTEQ2D (best higher DIS order fit). 

MSTV(22) ; (D=l) choice of nucleon structure-function library. 

= 1 : the internal one with parton distributions chosen according to MSTV(21). 

= 2 : the PDFLIB one with MSTV(21)=1000 • NGROUP + NSET (requires PDFLIB to be linked). 

MSTV(23) ; (D=l) choice of pion structure-function set (c.f. MSTV(24)). 

= 1 : Owens set 1. 

= 2 : Owens set 2. 


MSTV(24) ; (D=l) choice of pion structure-function library. 

= 1 : the internal one with parton distributions chosen according to MSTV(23). 

= 2 : the PDFLIB one with MSTV(23) = WOO-NGROUP + NSET (requires PDFLIB to be linked). 

MSTV(25) ; (D=2) choice of minimum Q-value Qq used in parton distributions of initial hadronic or 
nuclear state. 

= 0 : fixed at Qq =PARV(23). 

= 1 : determined by the average momentum tranfer of primary parton collisions Qo = {Qprim), with 
initial value Qo =PARV(23) before the first event. 

= 2 : as =1, but in addition varying with total energy y/s of the overall collision system, with Qq = 
max[Qo, f • (Eh/GeV)’’], where Qo =PARV(23), a =PARV(24), b =PARV(25), and Eh = V^/hadron 
the average energy per hadron (nucleon). 

MSTV(26) ; (D=l) option to choose x-dependent Qin in parton densities. 

= 0 : switched off, Qo according to selection of MSTV(25) is used. 

= 1 : varying with Bjorken-x of struck initial state partons, as well as with total energy ^/s of collision 

system, with Qo = Q(x) where 1/Q^(x) = I/Qq + Qo/{xEhad)‘^, and Eh = -v/i/hadron the average 
energy per hadron (nucleon). 

MSTV (31) ; (D=0) initial separation along the collision ( 2 :-)axis of beam and target particle in the overall 
CM frame. 

= 0 : fixed at Az =PARV(31). 

> 0 : shifted towards minimal separation, such that beam and target particles almost touch with 
Az =MSTV(31)-PARV(31). 
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MSTV(32) ; (D=2) selection of impact parameter b transverse to the collision (z-) axis of beam and 
target particle. 

= 0 : fixed at 6 = Max[PARV(32),PARV(33)]. 

= 1 ; randomly sampled in between bmin <b < bmax, with bmin =PARV(32) and bmax =PARV(33). 
= 2 ; as 1, but with bmin = PARV(32)-A/cj,irf/7r, and bmax = PARV(33)-a/ow 7^, where and is the 
non-diffractive crossection for given beam/target and beam energy. 

MSTV(33) : (D=2) Choice of nuclear matter distribution for collisions involving one or two nuclei, i.e. 
of nucleons’ positions within nucleus. 

= 0 : uniform distribution with ‘sharp edge’. 

= 1 ; Gaussian distribution for light nuclei A < 12, and ‘smeared edge’ distribution for intermediate 
and heavy nuclei A > 12 (c.f. PARV(27)). 

= 2 : Gaussian distribution for A < 12, and Fermi distribution for A > 12 (c.f. PARV(28)). 

MSTV(34) ; (D=l) spatial distribution of initial partons within their parent hadrons of beam and target 
particles, in the hadron restframe within a sphere of the hadron radius. 

= 0 ; uniform distribution. 

= 1 : exponential distribution with width given by PARV(34). 

= 2 ; Gaussian distribution with width given by PARV(35). 

MSTV(35) ; (D=2) Lorentz-boosted spatial distribution of initial partons. 

= 0 : spatial positions of all partons (valence, sea, glue) inside parent hadron are Lorentz-boosted 
with the parent hadron’s (nucleon’s) y-factor. 

= 1 : only valence quarks are boosted with the parent hadron’s (nucleon’s) y-factor; seaquarks and 
gluons are smeared out symmetrically(!) in hadron (nucleon) direction of motion around the valence 
quark disc. 

= 2 : as = 1, but now seaquarks and gluons are distributed all behind(!) the valence quark disc. 

MSTV(36) : (D=l) choice of primordial fex-distribution of initial partons. 

= 0 : no primordial k±. 

= 1 ; exponential /c_l- distribution with width specified by PARV(37). 

= 2 : Gaussian A:_l- distribution with width specified by PARV(38). 

MSTV(37) : (D=l) choice of masses of initial partons. 

= 0 : on mass shell with quark current masses. 

= 1 : off mass shell with space-like virtuality as determined by the excess of the partons’ momentum 
over their energies. 

MSTV(38) : (D=0) switch for including effects of parton shadowing in nuclear parton structure functions 
by multiplication with effective shadowing factor Ra{x) = fA{x)/fN{x) as a simple parametrization of 
EMG data, where fA and /at are the measured nuclear and nucleon structure functions, respectively. 
= 0 : parton shadowing switched off. 

= 1 ; switched on. 


PARV(21) : (D=l. GeV ^) the minimum resolvable partonic Bjorken-x, denoted Xc, is set by the 
‘confinement radius’ Rc =PARV(21), which defines the maximum possible longitudinal spread of a 
confined parton within Az = 1/xP < l/{xcP) = Ro where P is the longitudinal momentum of the 
mother hadron. 

PARV(22) : (D=0.95) the minimum summed x-value when ‘filling’ a hadron with partons of 

energy (momentum) fractions Xi of the hadron. The value depends, slightly on the value of Qo (c.f. 
PARV(23)-PARV(25)), or if cuts on x are made. R should be adjusted such that when adding one 
more parton, the average (X^i Xi) should be = 1. If the sum is smaller than PARV(22), another parton 
is added to the hadron, i.e. PARV(22) = 1 — {J2iXi)- 

PARV(23) : (D=l. GeV) minimum value Qo that is used in initial parton distributions (structure 
functions or alternative distributions). I.e., Q < Qo, then the value PARV(23) is used. 
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PARV(24) , PARV(25) : (D=2.50,0.25) possibility of using effective, energy-dependent Qin value in 
structure functions if MSTV(25)=2, according to the function Qq = max[Qo, f • {Eh/GeV)^], where 
Qo =PARV(23), a =PARV(24), b = PARV(25), and = -y/i/hadron the average energy per hadron 
(nucleon). 

PARV(26) : (D=l .18 fm) nuclear radius parameter vq, with R{A) = tqA^/3. 

PARV(27) : (D=3.0 fm) parameter Rd in ‘smeared-edge’ distribution of nucleons’ positions in nucleus 
with A > 12 (for MSTV(33) = 1). 

PARV(28) : (D=0.545 fm) parameter Ra in Fermi distribution of nucleons’ positions in nucleus with 
A > 12 (for MSTV(33) = 2). 

PARV(31) : (D=0.25 fm) initial separation along the collision ( 2 ;-)axis of beam and target particle in the 
overall CM frame for MSTV(31)> 1. 

PARV(32) : (D=0. fm) determines value of minimum impact parameter bmini transverse to the colli¬ 
sion ( 2 ;-)axis, such that bmin =PARV(32) for MSTV(32) = 0 or 1, or bmin =PARV(32)-A/ovid7^ for 
MSTV(32)=2. 

PARV(33) : (D=l. fm) determines value of maximum impact parameter bmax, transverse to the colli¬ 
sion (z-)axis, such that bmax =PARV(33) for MSTV(32) = 0 or 1, or bmax =PARV(33)-A/(Tnd/vr for 
MSTV(32)=2. 

PARV(34) : (D=1.19 GeV“^) width w in exponential distribution exp[—r/tc] of partons’ positions around 
its parent-hadron center-of-mass (c.f. MSTV(33)). (Note; w = l/i^, where n = 0.84 GeV is the 
measured nucleon form factor.) 

PARV(35) : (D=2.38 GeV“^) sigma s in Gaussian distribution exp[—r^/(2 s^)] of partons’ positions 
around its parent-hadron center-of-mass (c.f. MSTV(33)). (Note: s = 2/z^, where n = 0.84 GeV is 
measured nucleon form factor.) 

PARV(36) : (D=2.) cut-off for sampling initial partons’ spatial distributions from exponential or Gaus¬ 
sian above, such that r < aRh, where a =PARV(36) and Rh is the radius of the mother hadron. 

PARV(37) : (D=0.44 GeV) width w in exponential distribution exp[—fex/rc] of primordial k± distribution 
of initial partons (c.f. MSTV(36)). 

PARV(38) : (D=0.44 GeV) sigma s in Gaussian distribution exp(—/c^/(2 s^)] of primordial k± distribu¬ 
tion of initial partons (c.f. MSTV(36)). 

PARV(39) : (D=2. GeV) upper cut-off for sampling primordial k± from above exponential or Gaussian 
distribution. 


3. Parton scatterings: 


MSTV(40) : (D=3) master switch for parton collisions. 

= 0 : generally switched off. 

= 1 : only primary-primary collisions, i.e. those involving 2 partons that both have not interacted 
before. 

= 2 : as = 1 plus primary-secondary collisions, i.e. those involving 2 partons, of which one has had 
already at least one previous interaction. 

= 3 : as = 2 plus secondary-secondary interactions, i.e. all multiple collisions included. 
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MSTV(41) : (D=l) choice of (semi)hard QCD collision cross-sections. 

= 1 : standard 0{a‘^) perturbative cross-sections for 2 ^ 2 and 2 —> 1 processes. 

= 2 : alternative of 2 —> n cross-sections (not implemented yet!). 

MSTV(42) : (D=l) choice of po to separate hard and soft parton collisions, and to act as a regulator for 
the parton cross-sections da/dQ'^, that are divergent for vanishing momentum transfer (or invariant 
mass) Q. 

= 0 : fixed at po = max[PARV{A2),max{VKIN{3),VKIN{5)]. 

= 1 : energy and mass dependent, with pQ = pq(s,A), according to the parametrization given by 
PARV(42), PARV(43). 

= 2 : initially, before the first event, set to po{s,A), c.f. =1, but in subsequent events set equal 
to the average momentum transfer of primary parton-parton collisions po = {Qprim) accumulated by 
statistics. 

MSTV(43) : (D=2) Q^-definition in 2 —> 2 parton collision processes; for s-channel or fusion processes, 
is always chosen to be m? = s. Note that the parton-parton collisions invariants are denoted as 
s = (Pl +P2)^ i = {pi - p;)2, u = {p^- p'^f. 

= 1 •. = 2siu/{s^ + P + M^). 

= 2 : = 0.5 (m5_ i 

= 3 : = min(—t, —u). 

= 4 : Q'^ = s. 

MSTV(44) : (D=2) choice of scheme for parton collision-time estimate. 

= 0 : zero collision-time, i.e. instantaneous process. 

= 1 : finite collision-time sampled from an exponential distribution exp(—x), where x = t/r, t the 
time in the lab, and r the mean life-time of the parton in the lab, see PARV(44). 

= 2 : finite collision-time sampled from Gyulassy-Wang distribution (2/r) x/(l-|-x^)^, where x = t/r, 
t the time in the lab, and r the mean life-time of the parton in the lab, see PARV(44). 

MSTV(45) : (D=l) selection of maximum allowed distance between two partons in their cm-frame in 
order to be able to collide (used for ruling out widely separated pairs to speed up simulation). 

= 0 : fixed at a value given by ri 2 max = PARV(45). 

= 1 : variable at a value r 12 max = PARV(45)/(\/I/GeP). 

MSTV(46) : (D=0) choice of probability distribution W(x) from which collision probability of 2 partons 
is sampled, where x = jr^, with r[^ 2 ^^ the transverse 2-parton separation at closest approach 

and the transverse radius of interaction as given by the parton-parton cross-section a. 

= 0 : flat distribution, W{x) = 9{a — x), with a = PARV(46). 

= 1 : exponential distribution, W{x) = exp(—x/a), a = PARV(46). 

MSTV(47) : (D=2) special treatment of ‘flavor excitation’ by parton collisions involving one or two 
primary hadronic or nuclear (anti)quarks by requiring a minimum momentum transfer of the resolving 
parton with the heavy quark. 

= 0 : no special treatment of heavy quarks over light quarks. 

= 1 : requirement of minimum momentum transfer for liberation of heavy quark out of initial parton 
distribution by scattering, where the minimum required momentum trasnfer is given by PARV(47) 
times the mass of the heavy quark. 

= 2 : as = 1, but now both the struck heavy quark as well as its initial antiquark sibling are liberated 
(i.e. get on mass shell). 

MSTV(48) : (D=0) switch for including soft parton collisions according to a phenomenological cross- 
section dajdQ’^ oc as(Po)/(Q^ + parton collisions (c.f. eq. (HD), where po = PARV(42), 

p = PARV(48), and = p\ for t, it channel and = s for s channel. 

= 0 : switched off. 

= 1 : soft collisions are treated complementary to hard collisions, i.e. occur only if Q < po- 
= 2 ; soft collisions are treated supplementary to hard collisions, i.e. compete with the latter, 
whichever cross-section dominates. 
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MSTV(49) ; (D=l) switch to allow for 2^1 parton fusion processes a + h ^ c* in competition with 
2^2 scattering processes a+ h ^ c* ^d + e, such that for a given pair a,b either of those processes 
is selected according to the relative probability given by the respective cross-sections and depending 
on the ‘life-time’ of the particle c*. 

= 0 : fusion processes switched off. 

= 1 : switched on, with relative weight determined by PARV(49). 


PARV(41) : (D=0.25 GeV) nominal A-value used in running ag for parton-parton collisions. 

PARV(42) : (D=2.0 GeV) cut-off for perturbative QGD cross-sections that are divergent when the mo¬ 
mentum transfer Q ^ 0. Here the 2-body collsion scale Q is taken equal to p±, the transverse 
momentum exchanged in a parton scattering, or equal to s for annihilation or fusion processes. For 
MSTV(42)=0, the value po = a with a = PARV(42) defines hard {Q > po) and soft {Q < po) collisions. 
For MSTV(42)=1, see PARV(43). 

PARV(43) : (D=0.27) parameter for effective po value, if MSTV(42)=1, instead of po = const for 
MSTV(42)=0. The effective po assumes a dependence on the total collision energy ^/s and the mass 
of the (nuclear) collision system. R is parametrized as po(\/s, A, H) = | • [EhjGeV'f ^ where a = 
PARV(42), b = PARV(43), and Eh = 2^/s/{A + B) with A {B) the mass number of beam (target) 
particle, if it is a hadron or nucleus. 

PARV(44) : (D=0.1) if MSTV(44)=1 or MSTV(44)=2, the proportionality between mean collision time 
T of a 2-parton collision process, and the associated transverse momentum generated (for s-channel 
processes, the invariant mass of the pair). R is parametrized as r = c{\/Q), with c = PARV(45) and 
Q = Pa_ or s. 

PARV(45) : (D=5. GeV“^) maximum separation of two partons to be able to interact through a 2-body 
collision. For MSTV(45)=0, it is fixed at ri 2 = PARV(45), whereas for MSTV(45)=1, it is taken 
s-dependent as ri 2 = PARV(45)/(\/I/GeV). 

PARV(46) : (D=1.0) parameter in the collision-probability distribution W{x) = 6{a — x) (for 

MSTV(46)=0) and W{x) = exp(—x/a) for MSTV(46)=1, with x = relative 

transverse parton separation of closest approach and the transverse radius of interaction as given 
by the parton-parton cross-section d. 

PARV(47) : (D=1.0) for ‘flavor excitation’ processes where a primary hadronic or nuclear (anti)quark is 
struck out by a collision, an effective minimum resolution scale is required, depending on the flavor, 
such that the momentum transfer > a ■ nif, with a = PARV(47) and nif the flavor-dependent 
quark mass. This parameter acts in addition to the constraint set by PARV(42) above. 

PARV(48) : (D=0.5 GeV) parameter in the phenomenological cross-section dd/dQ"^ oc a^(pg)/((5^ -|- p?) 
for soft parton collisions, where p = PARV(48), po is given by PARV(42), and = p^ for t, u channel 
and = s for s channel. 

PARV(49) : (D=10.) parametric strength of 2 ^ 1 parton fusion processes a -|- 6 ^ c competing with 
2^2 scattering processes a + b ^ c + d. R is given by the ratio of the two processes which is 
proportional to l/(i?^s), where 1/i?^ =PARV(49). 


4. Space-like parton branchings: 


MSTV(50) : (D=l) master switch for space-like branchings. 
= 0 : switched off 

= 1 : on. 
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MSTV(51) ; (D=3) level of coherence imposed on the space-like parton shower evolution. 

= 1 : none, i.e. neither values nor angles need be ordered. 

= 2 : values at branches are strictly ordered, increasing towards the hard interaction. 

= 3 : values and opening angles of emitted (on-mass-shell or time-like) partons are both strictly 

ordered, increasing towards the hard interaction. 

MSTV(52) : (D=2) structure of associated time-like parton evolution, i.e. showers initiated by emission 
off the incoming space-like partons. 

= 0 ; no associated showers are allowed, i.e. emitted partons are put on mass-shell. 

= 1 : a shower may evolve, with maximum allowed time-like virtuality set by the phase space only. 
= 2 : a shower may evolve, with maximum allowed time-like virtuality set by phase space or by the 
scale Q^, the virtuality of the space-like parton created at the same vertex, whichever is the stronger 
constraint. 


PARV(51) : (D=0.25 GeV) A-value used in ag and in structure functions for space-like parton evolution. 

PARV(52) : (D=1.0 GeV) invariant mass cutoff qq of space-like parton showers, below which par- 
ton branchings re not evolved. For consistency it is taken as qo = max((5o,PARV(52)) where 
Qo =PARV(23) is the initial resolution scale in the hadron structure functions. 

PARV(53) : (D=4.) Q^-scale of the hard scattering is multiplied by PARV(53) to define the maximum 
parton virtuality allowed in space-like branchings associated with the hard interaction. This does not 
apply to s-channel resonances, where the maximum virtuality is set by m? = s. 


5. Time-like parton branchings: 


MSTV(60) : (D=2) master switch for time-like branchings. 

= 0 : no branchings at all, i.e. switched off. 

= 1 : QCD type branchings of quarks and gluons. 

= 2 : also emission of photons off quarks; the photons are assumed on mass-shell. 

MSTV(61) : (D=2) branching mode for time-like parton evolution. 

= 1 : conventional branching, i.e. without angular ordering. 

= 2 : coherent branching, i.e. with angular ordering. 

MSTV(62) ; not used. 

MSTV(63) : (D=l) choice of formation-time scheme for parton emission. 

= 0 : no formation time, i.e. instantanous emission. 

= 1 : finite formation time sampled from an exponential distribution exp(—x), where x = t/r, t the 
time in the lab, and r the mean life-time of the parton in the lab, see PARV(63). 

= 2 : finite formation time sampled from Gyulassy-Wang distribution (2/r) x/{l -|-x^)^, where 
X = t/r, t the time in the lab, and r the mean life-time of the parton in the lab, see PARV(63). 

MSTV(64) : not used. 

MSTV(65) ; (D=0) selection of kinematics reconstruction for branchings initiated by a single off-shell 
parton. 

= 0 : conservation of energy and jet direction, but longitudinal momentum is not conserved. 

= 1 : with full energy-momentum conservation at each branching, but the jet direction is not conserved 
due to recoil. 
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PARV(61) : (D=0.29 GeV) A-value used in running ag for time-like parton showers. 

PARV(62) : (D=1.0 GeV) invariant mass cutoff hq of timeloke parton evolution, below which partons are 
not assumed to radiate. Exception: for the case MSTV(84)=3, for reasons of consistency, the value is 
automatically scaled to 1.5 x PARV(62). 

PARV(63) : (D=0.01) if MSTV(63)=1 or MSTV(63)=2, the proportionality between mean-life time r of 
an off-shell parton and its virtuality /c^, is parametrized as r = c(l/\/^), with c = PARV(63). (Is 
inversely related to PARV(65)). 

PARV(64) : (D=4.) Q^-scale of the hard scattering is multiplied by PARV(64) to define the maximum 
parton virtuality allowed in time-like branchings associated with the hard interaction. This does not 
apply to s-channel resonances, where the maximum virtuality is set by = s. 

PARV(65) : (D=4.) invariant virtuality k‘^ of off-shell partons is multiplied PARV(65) to define the 
maximum parton virtuality allowed in time-like branchings of single parton decays. (Is inversely 
related to PARV(63)). 


6. Parton-cluster formation and cluster-hadron decay: 

MSTV(80) : (D=l) master switch for hadronization of parton clusters. 

= 0 : off. 

= 1 : on. 

MSTV(81) : (D=l) choice of definition of 2-parton spatial separation ri 2 = ri — r 2 between two partons 
1,2, for measuring probability of pre-hadronic cluster formation of the two partons. 

= 0 : distance ri2 measured in the global frame of the collision event. 

= 1 : distance ri2 measured in the center-of-mass frame of 1 and 2. 

MSTV(82) : (D=4) choice of definition of 2-parton momentum separation pi 2 of two partons 1, 2, for 
pre-hadronic cluster search. 

= 1 : three-momentum, i.e, pr, = n + pj u + pf kp 
= 2 : transverse momentum, i.e. pi 2 = 12 + Py 12 - 

= 3 : ‘pseudo-mass’, i.e. pi 2 = \/2 E 1 E 2 (1 — P 1 P 2 )/{'^ 11712 )■ 

= 4 : invariant mass', i.e, pm = - piEp^^ + pI u- 

MSTV(83) : (D=l) choice of probability distribution from which cluster formation of 2 partons is sam¬ 
pled. Simulates the conversion process as a ‘tunneling’ of the partons through a potential barrier set 
by PARV(81) and PARV(82), separating perturbative and non-perturbative regimes. 

= 0 : flat distribution, i.e. 0-function. 

= 1 : exponential distribution. 

MSTV (84) : (D=3) level of including color correlations among pairs of partons in the process of cluster 
formation. 

= 0 : none, except that clustering of two partons originating from the same mother is vetoed. 

= 1 : probability of two cluster candidates being in a color singlet state is sampled uniformly with 
factors Cgg = 1/9, Cqq = 1/3, Cgq = 1/3. 

= 2 : exact matching of colors and anticolors of two partons is required to give color singlet(s). 

= 3 : general case, where arbitrary color configuration of two partons is allowed, and color singlet 
formation is balanced by additional gluon emission(s) to conserve color locally at each vertex. 
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MSTV(85) : (D=0) option allowing for diquark clustering, i.e. sequential clustering of two (anti)quarks 
plus a third (anti)quark (not yet working properly!). 

= 0 : off. 

= 1 : on. 

MSTV(87) : (D=0) choice of effective invariant masses of partons to ensure that kinematical thresholds 
for cluster-hadron decay can be satisfied. 

= 0 : ‘curent’ masses are used, i.e. rrif = (0.05, 0.01, 0.01, 0 . 2 , 1.5, 5.0, 150) GeV for 
f = g,d, u, s, c, b, t, respectively. 

= 1 : ‘constituent’ masses are used, i.e. m/ = (0.5, 0.32, 0.32, 0.5, 1.8, 5.2, 150) GeV for 
f = g,d, u, s, c, b, t, respectively. 

MSTV(88) : (D=0) choice of measure for rough distinction between ‘slow’ and ‘fast’ clusters in statistics 
accumulation of exogamous cluster production with numerical value of boundary given by PARV( 88 ). 
= 0 : energy Ed- 
= 1 : 3-momentum Pd 
= 2 : cluster velocity (3 = Pd/Ed- 
= 3 : Lorentz factor 7 = Ed/Md- 

MSTV(89) : (D=l) smearing of cluster rapidities within an interval |y| < yc, where yc is determined by 
PARV(89). 

MSTV(90) : (D=l) master switch for hadronization of beam/target remnants. 

= 0 : off. 

= 1 : on. 

MSTV(91) : (D=l) switch for the decay of unstable hadrons, emerging either directly from cluster- 
decays, or from previous unstable hadron decays, 

= 0 ; off. 

> 1 : the decay chain is iterated MSTV(91) times, but at most 5 times. 

MSTV(92) : (D=l) switch for simulating baryon stopping effect in nucleus-nucleus collisions of protons 
and neutrons (and associated A’s) by redistributing them in transverse momentum p± according to 
a 1/(1 -|- bp^)^ distribution (approximately an exponential for p± < 2 GeV), and in longitudinal mo¬ 
mentum pz (or rapidity) according to a Gaussian distribution exp[—P^^g^^/(2c^)]. The parameters b 
and c are given by PARV(97) and PARV(98), respectively. 

= 0 : off. 

= 1 : on. 


PARV(81) : (D=3.6 GeV ^) minimum separation of 2 partons, below which a cluster cannot be formed 
(denoted in Ref. |^). 

PARV(82) : (D=4. GeV“^) critical separation of 2 partons that sets the space-time scale for cluster 
formation (denoted Rcrt in Ref. |^). Is related to the average cluster size and typically about 10 % 
above R^. 

PARV(83) : (D=1.0) ‘slope’ of the exponential probability distribution W{R) = 6 {R — 

R^) ex.p{—FARM{83)R/Rcrt) for MSTV(83)> 0, describing the probability of 2 partons to form a 
pre-hadronic cluster by tunneling through the ‘potential barrier’ marked by Ry. and Rcrt, with average 
cluster size Rd = 2Rcrt/PARV(83). 

PARV(84) : (D=0.3 GeV) minimum invariant mass of 2 partons required to form a cluster being neces¬ 
sarily above threshold for hadron decay. 

PARV(85) : (D=1000. GeV) maximum allowed invariant mass of 2 partons to form a cluster with mass 
Md , in order to provide option to suppress production of heavy clusters - in addition to the requirement 
of PARV(81),PARV(82). Note: default corresponds to no suppression. 
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PARV(86) : (D=5. GeV“^) maximum allowed separation of 2 partons that are unambigous cluster 
candidates. Cluster formation is enforced above this value, in order to avoid unphysical separation. 
In this case the requirement M^i < PARV(85) is overridden. 

PARV(87) : (D=l.) ‘slope’of exponential Hagedorn-type distribution for sampling the masses of clusters 
in their 2-body decays, given by exp(—a • M^i/Mq), where a = PARV(87), M^i is the cluster mass and 
Mq = 0.2 GeV is the ‘Hagedorn temperature’ (see Ref. [^). 

PARV(88) : (D=l. GeV) boundary for distinguishing ‘slow’ and ‘fast’ clusters in statistics accumulation 
of exogamous cluster production (c.f. MSTV(88)). 

PARV(89) : (D=0.2) determines value of tJc = PARV(89) Vheam for smearing of cluster rapidities within 
an interval \y\ < i/c, if MSTV(89) > 0. 

PARV(91) : (D=0.6) Parametric factor for charged multiplicity N^h from soft fragmentation of 

beam/target remnants with respect to pp collisions at CM energy ^/s, such that Nchi\^) = 

• ^/s), where a = PARV(91) (if MSTV(90) = 1). 

PARV(92) : (D=1.0) Nuclear dependence in eA or yA collisions of charged multiplicity Nch, resulting from 
soft fragmentation of nuclear remnant of the original nucleus with mass number A. It is parametrized 
as: N^\y/s,A) = A" N^\^/s), where a = PARV(92). 

PARV(93) : (D=1.0) Nuclear dependence in hA or AB collisions o the charged multiplicity N^h, re¬ 
sulting from soft fragmentation of nuclear beam or/and target remnants of the original nuclei with 
mass number A [B) (case of hadron h corresponds to A = 1 or R = 1). It is parametrized as: 

A B) = [{A + B)/2f where {3 = PARV(93). 

PARV(94) : (D=25.GeV“^) Time in the global collision frame when soft fragmentation starts to become 
gradually active (default corresponds to 5 fm.) 

PARV(95) : (D=250.GeV~^) The positions of particles produced in fragmentation of both parton clusters 
and beam/target clusters, are assigned such they become active only after an effective formation time 
At = PARV(95) X min(R/M^,PARV(95)) (default value corresponds At < 50 fm). 

PARV(96) : (D=10.) proportionality factor in the uncertainty relation At = PARV{96) x R/M^, which 
determines the mean formation time of particles produced in soft beam/target fragmentation (c.f. 
PARV(95)). 

PARV(97) : (D=0.35 GeV) parameter b in distribution for simulating baryon stopping effect 

in nucleus-nucleus collisions of protons and neutrons (and associated A’s) by redistributing them in 
transverse momentum p_i_ (c.f. MSTV(92)). 

PARV(98) : (D=0.14 GeV) parameter c in Gaussian distribution exp[—(2c^)] for simulating baryon 
stopping effect in nucleus-nucleus collisions of protons and neutrons (and associated A’s) by redistribut¬ 
ing them in longitudinal momentum p^ (c.f. MSTV(92)). 


7. Other settings: 


• Handling of errors and warnings: 

MSTV(IOO) : (D=2) check on possible errors during program execution. 

= 0 : errors do not cause any immediate action. 

= 1 : possible errors are checked and in case of problem, it is excited from the subprogram, but the 
simulation is continued from there on. For the first MSTV(lOl) errors a message is printed; after that 
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no messages appear. 

= 2 : possible errors are checked and in case of problem, the the simulation is forcibly terminated. 
For the first MSTV(lOl) errors a message is printed; after that no messages appear. 

MSTV(lOl) : (D=10) max number of errors that are printed. 

MSTV(102) : (D=l) printing of warning messages. 

= 0 : no warnings are written. 

= 1 : first MSTV(103) warnings are printed, thereafter no warnings appear. 

MSTV(103) : (D=10) max number of warnings that are printed. 

• Handling of input/output: 


MSTV(llO) : (D=l) direction of program output. 

= 0 ; all output is directed to standard output unit 6. 

= 1 ; output is directed to units MSTV(111)-MSTV(113). 


MSTV(lll) 


(D=10) unit number to which general output is directed. 


MSTV(112) 


(D=6) unit number for writing ‘on-line’ information. 


MSTV(113) 


(D=20) unit number for listing warnings and error messages. 


MSTV(114) : (D=2) switch for general information on unit MSTV(lll), which is directed automatically 
to file VNIRUN.DAT. Allows to select thr amount of output provided on the global performance of a 
simulation consisting of of a sample of collision events; 

= 0 ; only minimal output concerning selected process and main parameters is written out. 

= 1 ; as =0, plus a listing the 1st event in momentum and position space. 

= 2 : as =1, plus listing of number and properties of elementary subprocesses that occurred during 
the simulation. 


MSTV(115) : (D=2) switch for ‘on-line’ information on unit MSTV(112), i.e. writing of initialazation 
and termination info, and initial and final entries for each event as the simulation goes on. 

= 0 ; no ‘online output’. 

= 1 : only initialization and finalization info is printed. 

= 2 : continous listing of minimum event information to keep control of the program performance. 


• Miscellanous control switches for subroutines: 

MSTV(120) : (D=0) flag to indicate initialization of certain HERWIG common blocks and default values 
that are necessary for using parts of the HERWIG program for the hadronization of clusters. Is set 
equal to 1 after first intitialization call. 

MSTV(121) : (D=0) is set to 1 before a VNIROBO call, the V vectors are reset (in the particle range to 
be rotated/boosted) are set to 0 before the rotation/boost. MSTV(121) is inactive during a VNIROBO 
call and is set back to 0 upon return. 

MSTV(122) : (D=0) specifies in a VNIROBO call the type of rotation of the P, R, and V vectors. 
The rotation is clockwise (active rotation) for MSTV(122)=0 and anticlockwise (passive rotation) 
otherwise. The value of MSTV(122) is set back to 0 upon return. 

MSTV(123) : (D=0) specifies in a VNIROBO call for performing a boost, whether the vectors R are 
boosted (MSTV(123)=0) or not (MSTV(123)=1). The vectors P are always boosted. The value of 
MSTV(123) is set back to 0 upon return. 

MSTV(124) : (D=l) pointer to lowest entry in particle record VNIREC to be included in data analysis 
routines VNIANA1-VNIANA5. 
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MSTV(125) : (D=100000) pointer to highest entry in particle record VNIREC to be included in data 
analysis routines VNIANA1-VNIANA5. 

MSTV(126) : (D=20) number of lines in the beginning of the particle record that are reserved for internal 
event history information. The value should not be reduced, but can be increased if necessary. 

MSTV(127) : (D=l) in listing of the current state of the particle record by calling the subroutine 
VNILIST, it can be chosen between listing either color/anticolor C A and the space time origin of a 
particle I according to the information contained V{I,5) explained above, ( = 0 ), or, alternatively 
the mother KMO and color flow information contained in R(I, 3) — R(/, 5), ( = 1 ). All other listed 
quantities are the same for either option. 

MSTV(128) : (D=0) specifies in calls to VNIEDIT and VNIROBO the classes of particles to be included. 
The default = 0 eacludes all inactive or decayed particles with R(1,1) < 0 or K{1, 1) > 10, whereas 
the setting = 1 includes any entry listed in the particle record. MSTV(128) is set back to 0 upon 
return. 

• Version, date of last ehange: 

MSTV(140) : (D=l) Print-out of VNI logo on first occasion; MSTV(140) is reset to 0 afterwards. 

MSTV(141) : VNI version number. 

MSTV(142) : VNI subversion number. 

MSTV(143) : year of last change of VNI version 

MSTV(144) : month of last change of VNI version. 

MSTV(145) : day of last change of VNI version. 


46 



8. Statistics, event study and data analysis: 

• Subroutine VNILIST: Full listing of first event. 

MSTV(150) : (D=0) Gives an example listing of the evolution history of the first sucessful event in 
momentum and position space. Note; if the general output switch MSTV(114) is set > 1, then 
MSTV(150) is automatically set to 1. 

• Subroutine VNlANAl: General event statisties. 

MSTV(151) : (D=0) Statistics on initial parton state. 

MSTV(152) : (D=0) Number and momenta of produced particles. 

MSTV(153) : (D=0) Factorial moments. 

MSTV(154) : (D=0) Energy-energy correlations. 

MSTV(155) : (D=0) Decay channels. 

• Subroutine VN1ANA2: Statistics on hadronic observables. 

MSTV(161) : (D=0) Distributions of particles in y = ln(l/x). 

MSTV(162) : (D=0) Bose-Einstein correlation analysis. 

• Subroutine VN1ANA3: Statistics on pre-hadronic clusters. 

MSTV(171) : (D=0) Rapidity spectra of clusters diV/dy. 

MSTV(172) : (D=0) Longitudinal space-time spectra dA^/dz. 

MSTV(173) : (D=0) Transverse momentum spectra (1/p_l) dN/dp± 

MSTV(174) : (D=0) Transverse space-time spectra (l/r_L) dN/dr±. 

MSTV(175) : (D=0) Distributions of cluster sizes and masses. 

MSTV(176) : (D=0) Polarization profile of cluster density. 

• Subroutine VN1ANA4: Statistics on partons and produced hadrons. 

MSTV(181) : (D=0) Rapidity spectra of partons and hadrons dN/dy. 

MSTV(182) : (D=0) Longitudinal space-time spectra dA'/dz. 

MSTV(183) : (D=0) Transverse momentum spectra (1/p_l) dN/dp±. 

MSTV(184) ; (D=0) Transverse space-time spectra (l/rj,) dN/dr±. 
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• Subroutine VNIANA5: Global thermodynamic multipartiele properties. 

MSTV(191) : (D=0) Flavor, energy and momentum composition. 

MSTV(192) : (D=0) Flow velocity profiles. 

MSTV(193) : (D=0) Particle densities, energy densities, pressures. Note: in order to use this option, 
the flow-velocity profiles must! be obtained beforehand in a seperate run, because they are assumed to 
be read in as input (see MSTV(192)). 


MSTW(200), PARW(200): generated quantities and statistics 


1. General: 

MSTW(l) : total number of collision events to be generated. 

MSTW(2) : total number of time steps per collision event. 

MSTW(3) : physics process A B that is simulated. The current available beam (A) and target {B) 
particles and the collision processes are the ones listed before in Table 1, with MSTW(3)=IPRO. 

MSTW(4) : Overall Lorentz frame of reference for event siumulation. 

= 1 : center-of-momentum frame (CMS) of beam and target particle. 

= 2 : fixed-target frame (FIXT) with target particle at rest. 

= 3 : user-defined frame (USER) with given 3-momentum of beam and target. Particles are assumed 
on the mass shell. 

= 4 : user-defined frame (FOUR) with given f-momentum ( i.e., 3-momentum and energy) of beam 
and target particle. The particles need not to be on the mass shell. 

= 5 : user-defined frame (FIVE) with given 5-momentum ( i.e., 3-momentum, energy and mass) of 
beam and target particle. The particles need not to be on the mass shell, but 4-momentum and mass 
information must(!) match. 

MSTW(5) : KE flavour code for beam particle A. 

MSTW(6) : KE flavour code for target particle B. 

MSTW(7) : type of incoming beam particle A: 1 for lepton, 2 for hadron, and 3 for nucleus. 

MSTW(8) : type of incoming target particle B: 1 for lepton, 2 for hadron, and 3 for nucleus. 

MSTW(9) : combination of incoming beam and target particles. 

= 1 : lepton on lepton 
= 2 : lepton on hadron 
= 3 : lepton on nucleus 
= 4 : hadron on lepton 
= 5 : hadron on hadron 
= 6 : hadron on nucleus 
= 7 : nucleus on lepton 
= 8 : nucleus on hadron 
= 9 : nucleus on nucleus 
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MSTW(IO) : performance flag for current event. Is =0 at beginning, and set =1 if during the evolution 
the event turns out to be rejectable, either due to unphysical kinematics, particle combinations, etc., 
or due to numerical errors. 

MSTW(ll) : current collision event. 

MSTW(12) : current time step in this collision event. 

MSTW(13) : current number of successful collision events, i.e. those that completed gracefully with 
MSTW(10)=0. 

MSTW(14) : counter for ”non-diffractive” collisions events, i.e. those that involved at least one parton 
collision (in hadron or nucleus collisions). 

MSTW(15) : Status of Lorentz transformation between different global Lorentz frames in which the 
collision event is simulated (which may be different from the initially specified frame, c.f. MSTW(4)). 
The value of MSTW(15) saves the last performed transformation; 

= 1 ; from global cm-frame to fixed-target or user-specified frame; 

= 2 : from global cm-frame to hadronic center-of-mass in DIS (photon-hadron cm-frame); 

= -1 : from fixed-target or user-specified frame to global cm-frame; 

= -2 : from hadronic center-of-mass in DIS (photon-hadron cm-frame) to global cm-frame. 

PARW(l) : time increment TINC{I) in current time step I. 

PARW(2) : time TIME{I) in current time step I. 

PARW(3) : local machine time (in CPU) passed when simulation started. 

PARW(4) : local machine time (in CPU) when simulation ended. 

PARW(5) : conversion factor CPU ^ seconds (machine-dependent: for most processors, it is equal to 
0.01, i.e. 1 sec = 100 CPU. 

PARW(ll) : Total invariant -y/i = Eqm: i-e. the total CM energy of the collision system. 

PARW(12) : Invariant s = E^j^ mass-square of complete system. 

PARW(13) : CM energy per nucleon Ecm/{A + B) (where A {B) is the mass number of beam (target), 
and A + B is the total number of nucleons in the system) in collisions involving nuclei. Is equal to 
Ecm for elementary particle- and hadronic collisions. 

PARW(14) : Invariant s^ff = E‘^j^/{A + B) mass-square per nucleon in collisions involving nuclei. Is 
equal to s for elementary particle- and hadronic collisions. 

PARW(15) : mass Ma of beam particle. 

PARW(16) : mass Mb of target particle. 

PARW(17) : longitudinal momentum P^a of beam particle in specified global frame. 

PARW(18) : longitudinal momentum P^b of target particle in specified global frame. 

PARW(19) : angle 9 of rotation from CM frame to user-defined frame (also, in e^e via W^W , the 
rotation of the lU-pair along the z-axis). 

PARW(20) : azimuthal angle (j) of rotation from CM frame to user-defined frame, corresponding to 
PARW(18). 
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2. Initial state of collision system: 


MSTW(21) : number of neutrons of beam particle A. 

MSTW(22) : number of neutrons of target particle B. 

MSTW(23) : number of protons of beam particle A. 

MSTW(24) : number of protons of target particle B. 

MSTW(25) : number of cZ-valence quarks of beam particle A. 

MSTW(26) : number of d-valence quarks of target particle B. 

MSTW(27) : number of rt-valence quarks of beam particle A. 

MSTW(28) : number of rt-valence quarks of target particle B. 

MSTW(29) : number of s-valence quarks of beam particle A. 

MSTW(30) : number of s-valence quarks of target particle B. 

MSTW(31) : total number of initial state partons in collisions involving one or more hadron or nucleus. 
MSTW(32) : number of initial gluons. 

MSTW(33) : 10^ X number of d quarks + number of d antiquarks. 

MSTW(34) : 10^ x number of u quarks + number of u antiquarks. 

MSTW(35) : 10^ x number of s quarks + number of s antiquarks. 

MSTW(36) : 10^ X number of c quarks + number of c antiquarks. 

MSTW(37) : 10^ X number of b quarks + number of b antiquarks. 

MSTW(38) : 10^ X number of t quarks + number of t antiquarks. 


PARW(21) : absolute value of velocity (Sa of beam particle A. 

PARW(22) : absolute value of velocity Pb of target particle^. 

PARW(23) ; Lorentz factor 7^4 of beam particle A. 

PARW(24) : Lorentz factor 75 of target particle B. 

PARW(25) ; rapidity Ya of beam particle A. 

PARW(26) ; rapidity Yb of target particle B. 

PARW(27) : radius Ra (in 1/GeV) of beam particle in its restframe. 

PARW(28) : radius Rb (in 1/GeV) of target particle in its restframe. 

PARW(29) ; radius Rn (in 1/GeV) of neutron within in nucleus. 

PARW(30) : radius Rp (in 1/GeV) of proton within in nucleus. 

PARW(31) ; accumulated average Qq of initial hadron (nucleus) parton distributions (equals the average 
momentum scale of primary parton collisions {Qprim))- 
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PARW(32) : accumulated average momentum fraction xa_ of beam-side initial hadron (nucleus) parton 
distributions. 

PARW(33) : accumulated average momentum fraction xb oi target-side initial hadron (nucleus) parton 
distributions. 

PARW(34) ; 3 times the total charge of beam-side particle. 

PARW(35) : 3 times the total charge of target-side particle. 

PARW(36) : ratio b/bmax of actual impact parameter bmin < & < bmax to maximum allowed impact 
parameter of beam and target particles, if variable impact parameter is chosen (c.f. MSTV(32) and 
PARV(32), PARV(33)). 

PARW(37) ; azimuthal angle of colliding beam and target particles around the beam axis, if variable 
impact parameter is chosen. 

PARW(38) : sign of r^j-coordinate of beam-particle A in global frame for variable impact parameter 
selection (= — sign of target-particle B). 

PARW(39) : sign of r^^-coordinate of beam-particle A in global frame for variable impact parameter 
selection (= — sign of target-particle B). 


3. Parton scatterings: 


MSTW(40) : Total number of 2-parton collisions a + b ^ N {N = 1,2). 
MSTW(41) : Number of hard 2^2 collisions oi q + q,q + q,q + q. 

MSTW(42) : Number of hard 2^2 collisions oi q + g,q + g. 

MSTW(43) : Number of hard 2^2 collisions oi g + g. 

MSTW(44) : Number of soft 2^2 collisions oi q + q,q + q,q + q. 

MSTW(45) : Number of soft 2^2 collisions oi q + g,q + g. 

MSTW(46) : Number of soft 2^2 collisions of -|- ^f. 

MSTW(47) : Number of 2 ^ 1 fusions oi q + q. 

MSTW(48) : Number of 2 ^ 1 fusions oi q + g,q + g. 


MSTW(49) : Number of 2 ^ 1 fusions of -|- ^f. 


PARW(40) : total number of ‘primary’ (i.e., first) collisions among primary partons. 
PARW(41) : accumulated average of the ‘primary’ collisions. 

PARW(42) : accumulated average \/l of hard 2^2 collisions. 

PARW(43) : accumulated average rapidity y* = yl — 2/2 hard 2^2 collisions. 
PARW(44) : accumulated average p± of hard 2^2 collisions. 

PARW(45) : accumulated average \/I of soft 2^2 collisions. 
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PARW(46) 

PARW(47) 

PARW(48) 

PARW(49) 

PARW(50) 


accumulated average rapidity y* = yX — ^2 of soft 2^2 collisions, 
accumulated average p_\_ of soft 2^2 collisions, 
accumulated average \/l of hard 2^1 collisions, 
accumulated average rapidity y* = y\ — y^, of hard 2^1 collisions, 
integrated 2 -parton cross-section (t(s) = / dp\{da{s,p\)/dp\). 


4. Space-like parton branchings: 


MSTW(50) 

MSTW(51) 

MSTW(52) 

MSTW(53) 

MSTW(54) 

MSTW(58) 


Total number of space-like branchings a ^ h + c. 
Number of space-like processes CMshower ^ b + c. 
Number of space-like processes q{q) q{q) + g. 
Number of space-like processes g ^ g + g. 

Number of space-like processes g ^ q + q- 
Number of space-like processes q/q q{q) + 7 . 


PARW(50) 

PARW(51) 

PARW(52) 

PARW(53) 

PARW(54) 

PARW(55) 

PARW(56) 

PARW(57) 


average a/ Q'^ of starting scale for space-like branching evolution. 

average \/—q^ of virtuality of branching particle. 

average value of x in x ^ x'x” branching. 

average fraction z = x' jx of space-like branchings. 

average relative transverse momentum q±^ of space-like branchings. 

average longitudinal momentum q^ of branching particle. 

average energy of q^ branching particle. 

average invariant mass A/|g^| of branching particle. 
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5. Time-like parton branchings: 


MSTW(60) 

: Total number of time-like branchings a —> 6 -|- c. 

MSTW(61) 

: Number of time-like processes CMshower h + c. 

MSTW(62) 

: Number of time-like processes q{q) —> q{q) + g. 

MSTW(63) 

: Number of time-like processes g ^ g + g- 

MSTW(64) 

: Number of time-like processes g ^ q + q- 

MSTW(68) 

: Number of time-like processes q{q) —> q{q) + 7 . 

PARW(60) : 

average of starting scale for time-like branching evolution. 

PARW(61) : 

average of virtuality of branching particle. 

PARW(62) : 

average value of x in x ^ x'x” branching. 

PARW(63) : 

average energy fraction z = x'/x of time-like branchings. 

PARW(64) : 

average relative transverse momentum k± of time-like branchings. 

PARW(65) : 

average longitudinal momentum kz of branching particle. 

PARW(66) : 

average energy of branching particle. 

PARW(67) : 

average invariant mass Vm? of branching particle. 


6. Parton-cluster formation and clnster-hadron decay: 


MSTW(80) 

: Total number of 2-parton cluster formations. 

MSTW(81) 

: Total number of gg clusterings, gg —> CC + A, where X is either ‘nothing’, or g, or gg. 

MSTW(82) 

: Number of processes, gg —> CC + g. 

MSTW(83) 

: Number of processes, gg —> CC + gg. 

MSTW(84) 

: Total number of qq clusterings, qq —> CC + X, where X is either ‘nothing’, or g, or gg. 

MSTW(85) 

: Number of processes, qq —> CC + g. 

MSTW(86) 

: Number of processes, qq —> CC + gg. 

MSTW(87) 

: Total number of gq {gq) clusterings, gq ^ C + X, where X is either q or gq. 

MSTW(88) 

: Number of processes, gq —> CC + gq. 

MSTW(89) 

: Total number of {qq)q —> CC, {qq)q CC clusterings. 


MSTW(90) - MSTW(99) : Same as MSTW(80) - MSTW(89), but now for the corresponding numbers 
of ‘exogamously’ produced clusters, with ‘exogamy index’ Sd = (cj + ej)/2 unequal to 0 or 1 (see Ref. 
11 ). 
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MSTW(IOO) : Total number of ‘fast’ clusters, classified according to the choices of MSTV(88)] and 
PARV(88). 

MSTW(lOl) : Total number of ‘slow’ clusters, the fraction complimentary to the ‘fast’ one. 

MSTW(llO) : Total number of primary ’’neutral” particles per event, produced directly by cluster-hadron 
decays. 


MSTW(lll) 

MSTW(112) 

MSTW(113) 

MSTW(114) 

MSTW(115) 

MSTW(116) 

MSTW(117) 

MSTW(118) 

MSTW(119) 


Number of primary leptons and gauge bosons. 
Number of primary light mesons. 

Number of primary strange mesons. 

Number of primary charm and bottom mesons. 
Number of primary tensor mesons. 

Number of primary light baryons. 

Number of primary strange baryons. 

Number of primary charm and bottom baryons. 
Number of other particles. 


MSTW(120) - MSTW(129) : Same as MSTW(llO) - MSTW(119), but now for charged particles only. 

MSTW(130) - MSTW(139) : Same as MSTW(llO) - MSTW(119), but now for the numbers of sec¬ 
ondary neutral particles per event (i.e. those produced by decays of unstable particles). 

MSTW(140) - MSTW(149) : Same as MSTW(130) - MSTW(139), but now for charged particles only. 


PARW(81) : accumulated average space-time separation of clustered partons. 

PARW(82) : accumulated average energy-momentum separation of clustered partons. 

PARW(83) : minimum encountered separation of any two clustered partons with respect to chosen space- 
time measure (c.f. MSTV(81)). 

PARW(84) : maximum encountered separation of any two clustered partons with respect to chosen 
space-time measure (c.f. MSTV(81)). 

PARW(85) : minimum encountered separation of any two clustered partons with respect to chosen energy- 
momentum measure (c.f. MSTV(82)). 

PARW(86) : maximum encountered separation of any two clustered partons with respect to chosen 
energy-momentum measure (c.f. MSTV(82)). 

PARW(91) - PARW(99) : Gives the ‘exo(endo)-gamy distribution’ (see Ref. |M) of clusters formed 
from partons with ‘exogamy’ index ei,ej, such that Cd = {ei + ej)/2. Inital beam/target particles have 
e = 0(1), so that ed lies in the interval [0,1], which is binned as 0,0.05,0.2,0.3,0.45,0.55,0.7,0.8,0.95,1, 
with PARW(90)-PARW(99) giving the numbers of clusters in the bins. Note: only those clusters are 
counted that are classified as ‘fast’ (c.f. MSTW(88)). 

PARW(lOl) - PARW(109) : Same as PARW(91)-PARW(99), but for ‘slow’ clusters only. The sums 
PARW(91)-|-PARW(101), etc., hence give the total numbers. 

PARW(llO) - PARW(149) : not used. 
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8. Statistics, event study and data analysis: 


MSTW(151) - MSTW(160) : (D=l,21,41,61,91,121,151,181,221,301) Sequential numbers of time steps, 
at which a snapshot analysis during the space-time evolution of the collision system may be performed 
in routines VNIANA4 - VNIANA5. 

MSTW(161) : (D=40) Number of bins for energy distribution in ln(l/a;) (wher x is the energy fraction) 
in routine VNIANA2. 

MSTW(162) : not used. 

MSTW(163) : (D=40) Number of bins for pair mass distribution of Bose-Einstein correlations in routine 
VN1ANA2. 

MSTW(164)-MSTW(170) : not used. 

MSTW(171) : (D=20) Number of bins for cluster y-distribution in routine VNIANA3. 

MSTW(172) : not used. 

MSTW(173) : (D=20) Number of bins for cluster r^-distribution in routine VNIANA3. 

MSTW(174) : not used. 

MSTW(175) : (D=20) Number of bins for cluster p_L-distribution in routine VNIANA3. 

MSTW(176) : not used. 

MSTW(177) : (D=20) Number of bins for cluster r_i_-distribution in routine VNIANA3. 

MSTW(178) -MSTW(180) : not used. 

MSTW(181) : (D=20) Number of bins for parton/hadron y-distribution in routine VNIANA4. 
MSTW(182) : not used. 

MSTW(183) : (D=20) Number of bins for parton/hadron r^-distribution in routine VNIANA4. 
MSTW(184) : not used. 

MSTW(185) : (D=20) Number of bins for parton/hadron p_i_-distribution in routine VNIANA4. 
MSTW(186) : not used. 

MSTW(187) : (D=20) Number of bins for parton/hadron r_L-distribution in routine VNIANA4. 

MSTW(188) - MSTW(200) : not used. 


PARW(150) - PARW(160) : not used. 

PARW(161) ,PARW(162) : (D=0.,6.) Lower and upper bound for energy distribution in ln(l/x) (where 
X is the energy fraction of particles), in routine VNIANA2. 

PARW(163) ,PARW(164) ; (D=0.,3. GeV) Lower and upper bound for pair mass distribution of 
Bose-Einstein correlations in routine VNIANA2. 

PARW(165) - PARW(170) : not used. 
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PARW(171) ,PARW(172) : (D=-8.,8.) Lower and upper bound for cluster ^-distribution in routine 
VNIANA3. 

PARW(173) ,PARW(174) : (D=-40.,40. GeV“^) Lower and upper bound for cluster rz-distribution in 
routine VNIANA3. 

PARW(175) ,PARW(176) : (D=0.,10. GeV) Lower and upper bound for cluster p^-distribution in 
routine VNIANA3. 

PARW(177) ,PARW(178) : (D=0.,40. GeV“^) Lower and upper bound for cluster r^-distribution in 
routine VNIANA3. 

PARW(179) - PARW(180) : not used. 

PARW(181) ,PARW(182) : (D=-8.,8.) Lower and upper bound for parton/hadron y-distribution in 
routine VNIANA4. 

PARW(183) ,PARW(184) : (D=-100.,100. GeV“^) Lower and upper bound for parton/hadron re¬ 
distribution in routine VNIANA4. 

PARW(185) ,PARW(186) : (D=0.,10. GeV) Lower and upper bound for parton/hadron p^-distribution 
in routine VNIANA4. 

PARW(187) ,PARW(188) : (D=0.,100. GeV“^) Lower and upper bound for parton/hadron re¬ 
distribution in routine VNIANA4. 

PARW(189) - PARW(190) : not used. 

PARW(191) ,PARW(192) : (D=-20.,20. GeV“^) Lower and upper bound for r^ in (r^ x re) grid for 
density profiles in routine VNIANA5. ^-distribution in routine VNIANA4. 

PARW(193) ,PARW(194) : (D=0.,10. GeV“^) Lower and upper bound for re in (re x re) grid for 
density profiles in routine VNIANA5. 

PARW(195) - PARW(200) : not used. 


3.8 Kinematics cuts and selection of subprocesses 

The commonblock VNIKIN contains the arrays VKIN and MSUB, which may be used to impose specific 
kinematics cuts and to switch on/off certain subprocesses, respectively. The notation used in the following 
is that ’hat’ denotes internal variables in a partonic interaction, while ’*’ is for variables in the global CM 
frame of the system as-a-whole. All dimensionful quantities are in units of GeV or GeV“^. 

Kinematics cuts can be set by the user before the VNIXRUN call. Most of the cuts are directly related 
to the kinematical variables used in the evolution of the system in each event, and in the event data analysis. 
Except for the entries VKIN(7) - VKIN(12), the assignments follow the convention of JETSET/PYTHIA, 
in order to allow an easy interfacing. 

Note, that in the current version of VNI, only the values VKIN(7)-VKIN(12) are actively used. All 
other quantities are related to parton collision processes involving one or more partons, which are only of 
relevance in reactions with hadrons (nuclei) in the initial state. As stated before, the latter are presently 
worked out, but not yet included in the package. 
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C0MM0N/VNIKIN/VKIN(200),MSUB(200),ISET(200) 


Purpose: To impose kinematics cuts on the particles’ interactions and their space-time evolution. Also, to 
allow the user to run the program with a desired subset of processes. 

VKIN(l), VKIN(2) : (D=2.,-l.) range of allowed m = y/§ values. If VKIN(2) < 0, the upper limit is 
inactive. 

VKIN(3), VKIN(4) ; (D=2.,-l.) range of allowed (D=0.,-1.) range of allowed p± values for parton 
collisions, with transverse momentum p± is defined in the rest-frame of the 2-body interaction. If 
VKIN(4) < 0, the upper limit is inactive. For processes which are singular in the limit p_l ^ 0 (see 
VKIN(6)), VKIN(5) provides an additional constraint. The VKIN(3) and VKIN(4) limits can also be 
used in 2 ^ 1 ^ 2 processes. Here, however, the product masses are not known and hence assumed 
vanishing in the event selection. The actual p± range for massive products is thus shifted downwards 
with respect to the nominal one. 

VKIN(5) ; (D=l.) lower cutoff on p± values, in addition to the VKIN(3) cut above, for processes which 
are singular in the limit p± ^ 0 (see VKIN(6)). 

VKIN(6) : (D=l.) 2-body parton collision processes, which do not proceed only via an intermediate 
resonance (i.e. are 2 ^ 1 ^ 2 processes), are classified as singular in the limit p± 0, if either or 
both of the two final state products has a mass m < VKIN(6). 

VKIN(7), VKIN(8) : (D=-100.,100.) range of allowed particle rapidities y* in the CM frame of the 
overall system. 

VKIN(9), VKIN(IO) : (D=0.,1000.) range (in GeV) of allowed particle transverse momenta per¬ 
pendicular to the z-axis in the global CM frame of the event. 

VKIN(ll), VKIN(12) ; (D=-1000.,1000) range (in GeV“^) of allowed longitudinal positions z* of 
particles from the CM of the global collision system. 

VKIN(13), VKIN(14) : (D=0.,1000.) range (in GeV“^) of allowed transverse distances perpendic¬ 
ular to the z-axis in the CM frame of the global collision system. 


VKIN(21), VKIN(22) : (D=0.,1.) range of allowed xa (Bjorken-x) values for the parton of beam-side 
A that enters a parton collision. 

VKIN(23), VKIN(24) : (D=0.,1.) range of allowed xb (Bjorken-x) values for the parton of target-side 
B that enters a parton collision. 

VKIN(25), VKIN(26) : (D=-l.,l.) range of allowed Feynman xp values, where xp = xa — xb- 

VKIN(27), VKIN(28) : (D=-l.,l.) range of allowed cos(0) values in a 2 ^ 2 parton collision, where 6 
is the scattering angle in the rest-frame of the 2-body collision. 


VKIN(31), VKIN(32) : (D=2.,-l.) range of allowed m' values, where m' is the mass of the complete 
three- or four-body final state in 2 ^ 3 or 2 ^ 4 processes (while m (without a prime), constrained 
in VKIN(l) and VKIN(2), here corresponds to the one- or two-body central system). If VKIN(32) < 
0 , the upper limit is inactive. 

VKIN(35), VKIN(36) : (D=0.,-1.) range of allowed abs(t) = —i values in 2 ^ 2 processes. Note that 
for deep inelastic scattering this is nothing but the scale of the photon vertex, in the limit that 
initial and final state radiation is neglected. If VKIN(36) < 0, the upper limit is inactive. 
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VKIN(37), VKIN(38) : (D=0.,-1.) range of allowed abs(ii) = —u values in 2 ^ 2 processes. If 
VKIN(38) < 0, the upper limit is inactive. 


VKIN(39) - VKIN(IOO) : currently not used. 


VKIN(lOl) - VKIN(200) : reserved for internal use of storing kinematical variables of an event. 


MSUB(ISUB) : array to be filled when MSTV(14) = 1 (for MSTV(14) = 0, the array MSUB is initialized 
in VNIXRIN automatically) to choose which subset of subprocesses to include in the generation. The 
ordering follows the ISUB code given before in Table 3. 

= 0 : the subprocess ISUB is excluded. 

= 1 : the subprocess ISUB is included. 


ISET(ISUB) : gives the type of kinematical variable selection scheme used for subprocess ISUB. 
= 1:2^1 processes (parton fusion processes). 

= 2:2^2 processes (hard and soft parton scattering processes). 

= 3 : 1^2 processes (parton branching processes). 

= -1 : process is not yet implemented, but space is reserved. 

= -2 : process is not defined. 


3.9 Instructions on how to use the program 

Two example programs are included in the package VNI-3.1: vnixplel.f and vnixple2.f. The first one, 
vnixplel.f is disguised as a generic skeleton program that calls the steering routines VNIXRIN, VNIXRUN 
and VNIXFIN (c.f. Sec. 3.4) on a black-box-level. The second example vnixpl2. f illustrates more specifi¬ 
cally how to extract data from a simulation by using histograms. Below the simple example vnixplel .f is 
printed. This example generates 1000 pp events at ^/s = 546 GeV. Each event is evolved for a time range 
of 0 - 30 fm in the pp center-of-mass frame (the global CM frame). 

PROGRAM VNIXPLl 

C...Purpose: example for p+p~ collisions at CERN collider with 546 GeV. 

INCLUDE ’vnil.inc’ 

CHARACTER*(8) FRAME,BEAM,TARGET 


C...Required user input: 

C.number of events and final time (in fm). 

NEVT = 1000 
TFIN = 30. 

C.colliding particles, global Lorentz frame, and total energy. 

BEAM = ’p+ 

TARGET= ’p- 
FRAME = ’cms ’ 

WIN = 546. 

C.Specify your change of default parameters here: 
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C... Initialize simulation procedure. 

CALL VNIXRIN(NEVT,TFIN,FRAME,BEAM,TARGET,WIN) 


C...Begin loop over of events. 

DO 500 IEVT=1,NEVT 

C.Generate collision event. 

CALL VNIXRUN(IEVT,TFIN) 

C.Analysis of event as-a whole. 

C . 

C...End loop over events. 

500 CONTINUE 


C...Finish up simulation procedure. 
CALL VNIXFINO 


END 

This example is rather self-explanatory, and will become further clear with the following subsections. 
The ”...” indicate where one can interface on a black-box-level with the program as-a-whole without needing 
to dig deeper into its subroutine structure. In particular, one can insert here histograming (either using the 
included portable package ’’vnibook.f’, or any other package such as the HBOOK or PAW of the CERN 
library). Note also the pre-programmed histogram options discussed in Sec. 3.7 (switches MSTV(150)- 
MSTV(200)) and Appendix B (subroutines VNIANA1-VNIANA5). 

In the above example only the final state at the end of each event is collected. If the user is interested 
in the actual space-time development of the collision events, he/she can access this by inserting in the loop 
over events the following extension: 


C...Generate collision event; analyze evolution every 5 fm. 
TSTEP=5. 

NTIM=20 


C.Begin loop over time intervals. 

DO 400 ITIM=1,NTIM 
TRUN=ITIM*TSTEP 

IF(TRUN.GT.TFIN.OR.ITIM.EQ.NTIM) TRUN=TFIN 
CALL VNIXRUN(IEVT,TRUN) 


C.Analysis of time development. 

C . 

C.End loop over time intervals. 

400 IF(TRUN.EQ.TFIN) GOTO 450 
450 CONTINUE 


Remarks: 

(i) In case one desires to interface with the JETSET/PYTHIA programs, one should write a little program 
with common block 
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C0MM0N/LUJETS/N,K(4000,5),P(4000,5),V(4000,5) 


and which copies the variables N, K, P, and V of the particle record VNIREC directly. Be aware 
however, that VNIREC is dimensioned for 100000 entries, whereas LUJETS has space for only 4000 
entries. 

(ii) In case one desires to interface with the HERWIG program, one should write a little program with 
common block 

PARAMETER (NMXHEP=2000) 

COMMON/HEPEVT/NEVHEP,MHEP,ISTHEP(NMXHEP),IDHEP(MMXHEP), 

&JM0HEP(2,NMXHEP),JDAHEP(2,NMXHEP),PHEP(5,NMXHEP),VHEP(4,NMXHEP) 

and which calls the routine HWHEPC(MCOA^V, IP, NP) with MCONV = 1, IP = l,and NP = N. 
Again, watch the dimensions, 100000 versus 2000. 


3.10 Example for a typical collision event 

As a result of running the above example program vnixplel.f, a number of files are created in the directory, 
all of which are called VNI???.DAT, where ”???” are either 3 letters, or 3 digits. The latter are data files 
that contain event statistics and other information on the performance of a run. (see Sec. 3.7, switches 
MSTV(150)-MSTV(200) and Appendix B for further options). By default only a few files are created; 

VNIRUN.DAT containing the main summary of a simulation with print-out of results, 

VNICOL.DAT containing a detailed summary of all parton collisions that occurred, 

VNITIM.DAT containing the generated time-step grid and its real-time translation, 

VNIRLU.DAT containing the status of the random number generator, 

VNIERR.DAT containing error messages in case problems occur. 

Most important is the file VNIRUN.DAT, which gives a summary of what occurred during a simulation 
on the parton cascade level, the cluster formation level and the hadron decay level. One can obtain a listing 
of an event by simply calling VNILIST(I) at the end or during an event, which gives a first impression (for 
a complete listing see Appendix D): 

(i) The event record begins with the beam/target particles, here p and p, which are then transformed to 

the center-of-mass frame (here it coincides with the global frame of reference). 

(ii) Then follow the three valence quarks of p, and after that its intrinsic gluons and seaquark-pairs. After 

that the initial-sate valence quarks and gluons and seaquarks of p are listed. 

(iii) Then comes the history of the cascade evolution, including hard scatterings with shower 
(bremsstrahlung), subsequent cluster-formation of the materialized partons and hadron decay of these 
parton clusters. 

(iv) Finally, all non-materialized (non-interacted) partons are recombined to a beam- and a target cluster 
that undergo a soft fragmentation, again via cluster-hadron decay. 

(v) At the very end the sum of charge, three-momentum, energy, and the total invariant mass of the collision 

system is printed. 
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Example of the particle record for a pp collision event at ^/s = 546 GeV. 


Particle listing (summary: momentum space P in GeV) 
event no. 1 time step 200 at time 34.75 fm 


I 

particle 

KS 

KF 

C 

A 

P. 

_x 

P-y 

P. 

_z 

E 

M 

1 

(p+) 

15 

2212 

0 

0 


000 

.000 

272 

998 

273 

000 

.938 

2 

(p~-) 

15 

-2212 

0 

0 


000 

.000 

-272 

998 

273 

000 

.938 

3 

(p+) 

15 

2212 

0 

0 


000 

.000 

272 

998 

273 

000 

.938 

4 

(p~-) 

15 

-2212 

0 

0 


000 

.000 

-272 

998 

273 

000 

.938 

5 

(u) 

13 

2 

1 

0 


650 

.048 

1 

786 

2 

319 

.006 

6 

(d) 

13 

1 

1 

0 

- 

593 

-.165 

2 

092 

2 

587 

.010 

7 

(d~) 

13 

-1 

0 

3 

- 

135 

-.452 

6 

724 

6 

886 

-.490 

8 

(d) 

13 

1 

3 

0 


413 

-.211 


408 


935 

-.484 

9 

(g) 

14 

21 

1 

1 

1 

751 

.556 

2 

232 

3 

398 

-.682 

48 

(d~) 

13 

-1 

0 

2 


526 

-.150 

-16 

659 

15 

157 

.010 

49 

(u~) 

13 

-2 

0 

3 

- 

121 

-.114 

-4 

941 

4 

106 

-.174 

50 

(u) 

13 

2 

3 

0 

- 

252 

-.044 

-1 

576 


935 

-.249 

51 

(g) 

13 

21 

3 

2 


110 

.229 

-12 

733 

11 

449 

-.234 


83 

(g) 

14 

21 

2 

1 

2 

.052 

2 

.858 

1 

.999 

1 

.045 

-3 

,909 

84 

(d) 

14 

1 

1 

0 


.000 


.000 

-8 

.008 

8 

.008 


.000 

85 

(d) 

14 

1 

2 

0 


.000 


.000 

8 

.008 

8 

.008 


.000 

86 

(CMshower) 

11 

94 

0 

0 

2 

.661 

3 

.728 

-3 

.710 

8 

. 110 

5 

,571 

87 

(d) 

14 

1 

2 

0 

1 

.128 


.073 

-1 

.364 

3 

.504 

3 

,023 

88 

(cluster) 

11 

91 

0 

0 

1 

.116 

1 

.615 

-2 

.687 

3 

.412 


.750 

89 

(cluster) 

11 

91 

0 

0 

- 

.450 

- 

.605 

- 

.228 

1 

.213 


.922 

90 

(cluster) 

11 

91 

0 

0 

-1 

.528 

- 

.662 

1 

.381 

2 

.181 


.280 

91 

(cluster) 

11 

91 

0 

0 

1 

.252 


.689 

-1 

.272 

3 

. 136 

2 

,484 

92 

(cluster) 

11 

91 

0 

0 


.110 

- 

.650 

5 

.219 

5 

.796 

2 

,433 

93 

piO 

6 

111 

0 

0 


.660 


.523 

- 

.880 

1 

.171 


.135 

94 

piO 

6 

111 

0 

0 


.461 

1 

.090 

-2 

.094 

2 

.377 


.135 

95 

rho- 

6 

-213 

0 

0 

- 

.334 

- 

.562 

- 

.332 

1 

.068 


.769 

96 

piO 

6 

111 

0 

0 

- 

.110 

- 

.046 

- 

. 182 


. 186 


.135 

97 

pi- 

6 

-211 

0 

0 

- 

.678 

- 

.312 


.548 

1 

.066 


.140 

98 

pi+ 

6 

211 

0 

0 


.845 


.352 


.547 

1 

.203 


.140 

109 

(Beam-REM) 

11 

92 

0 

0 


.315 


.737 

274. 

.428 

281, 

.871 

70. 

, 164 

110 

(Targ-REM) 

11 

93 

0 

0 


.907 

-1. 

.587 

-276. 

.508 

247, 

.044 

-116, 

,181 

111 

(CMF) 

15 

100 

0 

0 

-1, 

.222 


.850 

-2. 

.080 

528, 

.915 

528, 

,909 

112 

(cluster) 

11 

91 

0 

0 

- ^ 

.315 


.737 

274. 

.428 

281, 

.871 

70, 

, 164 
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113 

(cluster) 

11 

91 

0 

0 

- 

.907 

-1 

.587 

-276, 

.508 

247, 

.044 

-116 

. 181 

114 

(Beam-REM) 

11 

92 

0 

0 


.411 


.648 

127, 

.485 

127, 

.543 

3 

.762 

115 

(a_l-) 

17 - 

-20213 

0 

0 


.493 


.228 

19, 

.200 

19. 

,250 

1 

.275 

116 

(Delta+) 

17 

2214 

0 

0 

- 

.081 


.420 

108, 

.286 

108, 

.293 

1 

.232 

117 

piO 

7 

111 

0 

0 


.490 


.124 

4. 

,473 

4. 

,831 


.135 

118 

(rho-) 

17 

-213 

0 

0 


.005 


.102 

14, 

.584 

14. 

,605 


.769 

119 

piO 

7 

111 

0 

0 

- 

.160 


.015 

35, 

.475 

37. 

,042 


.135 

120 

P+ 

7 

2212 

0 

0 


.085 


.402 

72, 

.525 

75. 

,579 


.938 

121 

piO 

7 

111 

0 

0 

- 

.122 

- 

.268 

8. 

.519 

9. 

.014 


.135 

122 

pi- 

7 

-211 

0 

0 


.133 


.367 

5. 

,779 

6. 

.174 


.140 

123 

(Targ-REM) 

11 

93 

0 

0 

-1 

.252 

- 

.510 

-63, 

.766 

63. 

,825 

2 

.384 

124 

(Delta'-) 

17 

-2214 

0 

0 

- 

.779 

- 

.303 

-21, 

.201 

21. 

.253 

1 

.232 

125 

(eta) 

17 

221 

0 

0 

- 

.473 

- 

.208 

-42, 

.565 

42. 

,572 


.549 

126 

piO 

7 

111 

0 

0 

- 

.072 

- 

.254 

-3, 

.205 

3. 

.199 


.135 

127 

P~- 

7 

-2212 

0 

0 

- 

.701 

- 

.052 

-18, 

.281 

18. 

,903 


.938 

128 

piO 

7 

111 

0 

0 

- 

.179 


.069 

-12, 

.277 

12. 

,622 


.135 

129 

piO 

7 

111 

0 

0 

- 

.060 

- 

.060 

-12, 

.901 

13. 

,269 


.135 

130 

piO 

7 

111 

0 

0 

- 

.225 

- 

.220 

-17, 

.816 

18. 

,383 


.135 


336 

(rho-^) 

17 

213 

0 

0 

-.263 

.404 

-5, 

.474 

5 

.549 

.769 

337 

(omega) 

17 

223 

0 

0 

.067 

-.237 

-31, 

.485 

31 

.496 

.783 

338 

piO 

7 

111 

0 

0 

-.342 

.522 

-5, 

.018 

5 

.113 

.135 

339 

Pi+ 

7 

211 

0 

0 

.085 

-.120 


,742 


.657 

.140 

340 

piO 

7 

111 

0 

0 

-.309 

.029 

-10, 

.513 

10 

.790 

.135 


sum: .00 .000 .000 .000 546.000 546.000 


More information on pre-programmed print-out of results concerning for particle distributions, time 
evolution, observables, etc., can be switched on with switches MSTW(151) - MSTW(200), as explained in 
the previous Sec. 3.7. 
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A Further physics routines 

Aside from the three main routines VNIXRIN, VNIXRUN, and VNIXFIN, the physics routines described 
below form the heart of the program. In general, for each of the included collision processes of Table 1, there 
is an initialization routine VNI??IN that sets up the initial state, and an evolution routine VNIEV?? that 
carries out the time evolution of the particle distributions in phase-space, starting from the initial state. The 
routines VNICOLL, VNISBRA, VNITBRA, etc., therein perform the perturbative parton evolution in terms 
of parton collisions, space-like branchings, and time-like branchings, respectively. The routines VNICHAD, 
VNIBHAD, etc., deal with the cluster hadronization. These routines are universal and independent of the 
process under consideration. 

SUBROUTINE VNIC0L2(NCOL) 

Purpose: to test for all possible pairings of ‘alive’ partons, whether a collision is allowed, and, if it is allowed, 
to evaluate the quantities characterizing the 2-body collision at the parton level according to the cross-section 
and to choose one of the possible subprocesses (channels). 

NCOL : total number of found collisions. 


SUBROUTINE VNICOLL(ISUB,Q2,IPl,IP2,IPS,IP4) 

Purpose: to perform a 2-body parton collision, find outgoing flavors and colors and calculate kinematics 
and parton-parton collision. 

ISUB : Subprocess of 2-body collision (c.f. Table 3). 

Q2 : Momentum transfer scale of subprocess (c.f. MSTV(43), Sec. 3.7). 

IPl, IP2 : Line numbers of paerons entering the interaction vertex. 

IPS, IP4 : Line numbers of particles emerging from the interaction vertex. For 2^1 processes I PA is 
equal to zero. 


SUBROUTINE VNISBRA(IP1,IP2,QMAX) 

Purpose: to generate space-like parton evolution, by backward tracing of branching processes down to some 
initial resolution scale associated with a non-perturbative cut-off. The performance is regulated by the 
switches MSTV(50)-MSTV(52) and parameters PARV(51)-PARV(53) (c.f. Sec. 3.7). 

IPl, IP2 : partons in lines IPl and IP2 in the commonblock VNIREC, are evolved by space-like parton 
branchings from initial mass scale QMAX. If one of the particles is a non-radiating one, then only one 
branch is evolved, and the spectator particle role is then only to take up recoil due to emitted particles 
by the radiating one. 

QMAX : initial mass scale, from which the two particles are evolved from Q = QMAX backwards to 
lower Q < QMAX until the cut-off Q = qo =PARV(52) is reached. 


SUBROUTINE VNITBRA(IP1,IP2,QMAX) 

Purpose: to generate time-like parton evolution, by coherent branching processes within a finite time interval, 
Q ^ Q9, Q ^ 7, 9 ^ 99 1 9 QQ- The performance is regulated by the switches MSTV(61)-MSTV(65) and 
parameters PARV(61)-PARV(65) (c.f. Sec. 3.7). 
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IPl > 0 IP2 = 0 : generate a time-like parton branching for the parton in line IPl in commonblock 
VNIREC, with maximum allowed mass QMAX. 

IPl > 0 IP2 > 0 : generate time-like parton branchings for the two partons in lines IPl and IP2 in the 
commonblock VNIREC, with maximum allowed mass for each parton QMAX. 

QMAX : the maximum allowed mass of a radiating parton, i.e. the starting value for the subsequent 
evolution. (In addition, the mass of a single parton may not exceed its energy, the mass of a parton 
in a system may not exceed the invariant mass of the system.) 


SUBROUTINE VNICHAD(I1,12) 

Purpose: to search at a given point of time the particle record for possible clustering of materialized partons, 
which are those which have been struck out from the initial state or have been produced as secondaries by 
parton collisions or radiative processes. It consists of checking for each possible pair of two partons, the 
invariant space-time distance between them, their flavor and origin, the matching of color and anticolor, 
their common invariant mass, etc.. If a pair satisfies all necessary conditions for the formation of one or two 
color singlet clusters, then they are combined and converted into a composite cluster with internal flavor 
content and invariant mass determined by the partons, which then decays either into a pair of color singlet 
clusters, or into a singlet cluster plus an additional gluon or quark. 

71, 72 : first and last entry of particle record VNIREC in between the cluster formation is carried out by 
looping over all particles 7 with 71 < 7 < 72. 


SUBROUTINE VNICFRA(I1,12) 

Purpose: to fragment color singlet clusters with given qq substructure into final state hadrons. Eor each 
cluster its invariant mass and flavor content are checked, and it is determined whether it decays isotropically 
into a pair of hadrons, or converts directly into a single hadron of corresponding flavor. The hadron decay 
is described by a phase-space model, in which the probability for forming a certain hadron state is given by 
the density of hadronic states with mass below the cluster mass. Once a hadron state is chosen, it either 
propagates on undisturbed if it is a stable particle, or it decays further according to the particle data tables, 
if it is a resonance. 

71 > 0 : entry of particle record VNIREC from which analysis starts. 

72 > 72 : entry of particle record up to which analysis extends. 


SUBROUTINE VNIBHAD(I1,12) 

Purpose: to perform beam cluster formation and hadronization. The remnant partons of beam/target 
particles, those that originate from the partonic coherent initial state in collisions involving hadrons or 
nuclei and have not interacted with decohereing effect. These primary partons are reassembled to clusters 
along the beam axis, which are then fragmented based on the minimum bias model of the UA5 collaboration, 
such that color, flavor and energy-momentum of the collision system as-a-whole are conserved. 

71, 72 : first and last entry of particle record VNIREC in between the soft beam/target fragmentation is 
carried out by looping over all particles 7 with 71 < 7 < 72. 
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SUBROUTINE VNIBFRA(I1,12) 


Purpose: to hadronize beam/target clusters consisting of remnant beam partons, using parts of the HERWIG 
5.7 program. Beam remnants are those that emerge from the coherent partonic initial state without having 
undergone any decohering interaction, and carry the remaining energy, flavor and color of the system. The 
beam remnants are lumped to beam clusters which are fragmented similarly as the clusters arising from the 
parton cascade evolution. 

71 > 0 : entry of particle record VNIREC from which analysis starts. 

72 > 71 : entry of particle record up to which analysis extends. 


SUBROUTINE VNIPADC(I1,12,NUST) 

Purpose: to decay unstable particles resulting from cluster decays to hadrons. Out of these primary hadrons, 
those that are unstable resonances are decayed according to the Particle Data, which gives lower mass 
hadrons plus possibly gauge bosons and leptons. These secondary decay products may decay further. 

71 > 0 : entry of particle record VNIREC from which analysis starts. 

72 > 71 : entry of particle record up to which analysis extends. 

NUST : returns number of decay products that are still unstable. 
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B Event analysis routines 


The following routines collect, analyze and print out information on the result of a simulation, and give the 
user immediate access to calculated particle data, spectra, and other observables. The following routines 
are adopted (with some extensions and modifications) from the JETSET/PYTHIA program. 


• The routine VNILIST provides a useful tool to obtain a quick overall view of the state of the particle 
record at any time during the simulation of one or many events, and (in parts) the space-time history 
of the system. 

• The routine VNIEDIT allows to edit the particle record, e.g. to compress, or to exclude certain particle 
species from the record. It is used frequently during the simulation procedure to clean up used but 
further on unecessary particle information. 

• Via the functions KVNI and PVNI, values of some frequently appearing variables may be obtained 
more easily. 

• The routines VNIANAL and VNIANAl - VNIANA5 carry out various data accumulation duties, 
and print out tables at the end of the simulation, which are automatically directed to files named 
VNI???.DAT where ??? is a 3-digit number, specified below. 


Einally, the six routines VNISPHE, VNITHRU, VNICLUS, VNICELL, VNIJMAS and VNIEOWO 
allow the possibility to find some global event shape properties as discussed in detail in the JET- 
SET/PYTHIA documentation [41|. 


SUBROUTINE VNILIST(MSEL) 

Purpose: to list the current state of the particle record during the evolution of an event, or certain specified 
excerpts of it. Also, to list parton or particle data, or current parameter values. 

MSEL : determines what is to be listed. 

= 0 : writes a logo with program version number and last date of change; is mostly for internal use. 
= 1 : gives a simple list of current event record, in an 80 column format suitable for viewing directly 
on the computer terminal. Eor each entry, the following information is given; the entry number I, 
the parton/particle name, the status code KS (K(I,1)), the flavour code KF (K(I,2)), the color and 
anticolor, C (L(I,1)), and A (L(I,2)), and the three-momentum, energy and mass (P(I,1)-P(I,5)). A 
final line contains information on total charge, momentum, energy and invariant mass. 

= 2 : gives a more extensive list of the current event record, in a 132 column format, suitable for 
printers or workstations. Eor each entry, the following information is given: the entry number I, the 
parton/particle name (with padding as described for MSEL = 1), the status code KS (K(I,1)), the 
flavour code KF (K(I,2)), the color and anticolor, C (L(I,1)) and A (L(I,2)), the origin of the particle 
in a 3 digit N1N2N2>, where N1 is the event and N2 the time step in this event, when the particle 
was first produced, and N2> is the number of the entry in the particle record that the particle occupied 
when it first appeared. This labeling allows a unique specification of its ‘historical’ origin. Einally, it 
follows the three-momentum, energy and mass (P(I,1) - P(I,5)). A final line contains information on 
total charge, momentum, energy and invariant mass. 

= 3 : gives the same basic listing as = 2, but with an additional line for each entry containing infor¬ 
mation on production vertex position and time (V(I,1) - V(I,4)). 

= -1 : gives the format of listing as = 1 , but instead of momentum, energy and mass, it lists now the 
particle’s position in the overall reference frame, (R(I,1) - R(I,3)), its momentum-rapidity (Y(P)) and 
its transverse momentum (PT) with respect to the z-axis, which is usually the jet-axis, or beam axis. 
= -2 : gives the format of listing as = 2 , but as for = - 1 , instead of momentum, energy and mass, 
now the particle’s positions, rapidities and transverse momenta in the overall Lorentz frame are listed. 
As for = 2 , additional information on the ‘historical’ origin of the particle is listed. 

= -3 : gives the format of listing as = 3, again, instead of momentum, energy and mass, now the 
particle’s positions, rapidities and transverse momenta are listed. 
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= 11 : provides a simple list of all parton/particle codes defined in the program, with KF code 
and corresponding particle name. The list is grouped by particle kind, and only within each group in 
ascending order. 

= 12 : provides a list of all parton/particle and decay data used in the program. Each parton/particle 
code is represented by one line containing KF flavour code, KC compressed code, particle name, an¬ 
tiparticle name (where appropriate), electrical and color charge, mass, resonance width and maximum 
broadening, average invariant lifetime and whether the particle is considered stable or not. 

= 13 : gives a list of current parameter values for MSTV, PARV, and for MSTW and PARW. This 
is useful to keep check of which default values were changed in a given run. 

= 14 : gives a list of time grid and size of time slizes used in the current run. 


SUBROUTINE VNIEDIT(MEDIT) 

Purpose: to exclude unstable or undetectable partons/particles from the event record. One may also use 
VNIEDIT to store spare copies of events (specifically initial parton configuration) that can be recalled to 
allow e.g. different fragmentation schemes to be run through with one and the same parton configuration. 
Einally, an event which has been analysed with VNISPHE, VNITHRU or VNICLUS (see below) may be 
rotated to align the event axis with the z-direction. 

MEDIT : tells which action is to be taken. 

= 0 ; empty {KS = 0) and documentation {KS > 20) lines are removed. The partons/particles 
remaining are compressed in the beginning of the VNIREC commonblock and the N value is updated 
accordingly. The event history is lost, so that information stored in K(I,3), K(I,4) and K(I,5) is no 
longer relevant. 

= 1 : as = 0, but in addition all partons/particles that have fragmented/decayed {KS > 10) are 
removed. 

= 2 : as = 1, but also all neutrinos and unknown particles (i.e. compressed code KC = 0) are 
removed. 

= 3 : as = 2, but also all uncharged, color neutral particles are removed, leaving only charged, stable 
particles (and unfragmented partons, if fragmentation has not been performed). 

= 5 ; as = 0, but also all partons which have branched or been rearranged in a parton shower and 
all particles which have decayed are removed, leaving only the fragmenting parton configuration and 
the final state particles. 

= 11 ; remove lines with K(I,1) < 0. Update event history information (in K(I,3) - K(I,5)) to refer 
to remaining entries. 

= 12 : remove lines with K(I,1) = 0. Update event history information (in K(I,3) - K(I,5)) to refer 
to remaining entries. 

= 13 ; remove lines with K(I,1) = 11, 12 or 15, except for any line with K(I,2) = 94. Update event 
history information (in K(I,3) - K(I,5)) to refer to remaining entries. In particular, try to trace origin 
of daughters, for which the mother is decayed, back to entries not deleted. 

= 14 : remove lines with K(I,1) = 13 or 14, and also any line with K(I,2) = 94. Update event history 
information (in K(I,3) - K(I,5)) to refer to remaining entries. In particular, try to trace origin of 
rearranged partons back through the parton shower history to the shower initiator. 

= 15 : remove lines with K(I,1) > 20. Update event history information (in K(I,3) - K(I,5)) to refer 
to remaining entries. 

= 21 : all partons/particles in current event record are stored (as a spare copy) in bottom of com¬ 
monblock VNIREC. 

= 22 ; partons/particles stored in bottom of event record with = 21 are placed in beginning of record 
again, overwriting previous information there (so that e.g. a different fragmentation scheme can be 
used on the same partons). Since the copy at bottom is unaffected, repeated calls with = 22 can be 
made. 

= 23 : primary partons/particles in the beginning of event record are marked as not fragmented or 
decayed, and number of entries N is updated accordingly. Is simple substitute for = 21 plus = 22 
when no fragmentation/decay products precede any of the original partons/particles. 
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= 31 : rotate largest axis, determined by VNISPHE, VNITHRU or VNICLUS, to sit along the z- 
direction, and the second largest axis into the xz plane. For VNICLUS it can be further specified to 
+z axis and xz plane with x > 0, respectively. Requires that one of these routines has been called 
before. 

= 32 : mainly intended for VNISPHE and VNITHRU, this gives a further alignment of the event, in 
addition to the one implied by = 31. The ’’slim” jet, defined as the side (z > 0 or z < 0) with the 
smallest summed p±_ over square root of number of particles, is rotated into the +z hemisphere. In 
the opposite hemisphere (now z < 0 ), the side of x > 0 and x < 0 which has the largest summed pz 
absolute is rotated into the z < 0, x > 0 quadrant. Requires that VNISPHE or VNITHRU has been 
called before. 


FUNCTION KVNI(I,J) 

Purpose: to provide various integer-valued event data. Note that many of the options available (in particular 
/ > 0, J > 14) which refer to event history will not work after a VNIEDIT call. 

I = 0,J = : properties referring to the complete event. 

= 1 : N, total number of lines in event record. 

= 2 ; total number of partons/particles remaining after fragmentation and decay. 

= 6 : three times the total charge of remaining (stable) partons and particles. 


7 > 0, J = : properties referring to the entry in line no. I of the event record. 

= 1 - 5 : K(I,1) - K(I,5), i.e. parton/particle status KS, flavour code KF and origin/decay 
product/color-flow information. 

= 6 : three times parton/particle charge. 

= 7:1 for a remaining entry, 0 for a decayed/fragmented/ documentation entry. 

= 8 : KF code (K(I,2)) for a remaining entry, 0 for a decayed/fragmented/documentation entry. 

= 9 ; KF code (K(I,2)) for a parton (i.e. not color-neutral entry), 0 for a particle. 

= 10 : KF code (K(I,2)) for a particle (i.e. color-neutral entry), 0 for a parton. 

= 11 : compressed flavour code KC. 

= 12 ; color information code, i.e. 0 for color neutral, 1 for color triplet, -1 for antitriplet and 2 for 
octet. 

= 13 : flavour of’’heaviest” quark or antiquark (i.e. with largest code) in hadron or diquark (including 
sign for antiquark), 0 else. 

= 14 : generation number. Beam particles or virtual exchange particles are generation 0, original 
partons/particles generation 1 and then 1 is added for each step in the fragmentation/decay chain. 

= 15 ; line number of ancestor, i.e. predecessor in first generation (generation 0 entries are disre¬ 
garded) . 

= 16 : rank of a hadron in a jet system it belongs to. Rank denotes the ordering in flavour space, with 
hadrons containing the original flavour of the jet having rank 1 , increasing by 1 for each step away 
in flavour ordering. All decay products inherit the rank of their parent. Whereas the meaning of a 
first-rank hadron in a quark jet is always well-defined, the definition of higher ranks is only meaningful 
for independently fragmenting quark jets. In other cases, rank refers to the ordering in the actual 
simulation, which may be of little interest. 

= 17 : generation number after a collapse of a jet system into one particle, with 0 for an entry not 
coming from a collapse, and -1 for entry with unknown history. A particle formed in a collapse is 
generation 1 , and then one is added in each decay step. 

= 18 : number of decay/fragmentation products (only defined in a collective sense for fragmentation). 
= 19 : origin of color for showering parton, 0 else. 

= 20 : origin of anticolor for showering parton, 0 else. 

= 21 : position of color daughter for showering parton, 0 else. 

= 22 : position of anticolor daughter for showering parton, 0 else. 

= 23 : number of 2-body collisions of parton/particle. 

= 24 : number of space-like decays of parton/particle. 
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= 25 : number of time-like decays of parton/particle. 

= 27 : entry of direct mother of parton/particle. 

= 28 : timestep of production of parton/particle. 

= 29 : index of generation in the cascade tree of parton/particle. 


FUNCTION PVNI(I,J) 

Purpose: to provide various real-valued event data. Note that some of the options available (/ > 0, J = 
20 — 25) are primarily intended for studies of systems in their respective CM frame. 

/ = 0,J = : properties referring to the complete event. 

= 1 - 4 : sum of Px,Py,Pz and E, respectively, for the stable remaining entries. 

= 5 : invariant mass of the stable remaining entries. 

= 6 ; sum of electric charge of the stable remaining entries. 

/ > 0, J = : properties referring to the entry in line no. I of the event record. 

= 1 - 5 : P(I,1) - P(I,5), i.e. noimally Px,Py,Pz, E and m for parton/particle. 

= 6 : electric charge q. 

= 7 ; momentum squared = Px + Py + pi- 
= 8 : momentum p = 

= 9 ; transverse momentum squared p\ = p‘^+ Py. 

= 10 : transverse momentum p_i_. 

= 11 : transverse mass squared = m? +p^+Py- 
= 12 : transverse mass m±. 

= 13 - 14 ; polar angle 9 in radians (between 0 and tt) or degrees, respectively. 

= 15 - 16 : azimuthal angle (j) in radians (between — tt and tt) or degrees, respectively. 

= 17 : true rapidity y = ^ln((£' +pz)l{E — pz))- 

= 18 : rapidity i/jr obtained by assuming that the particle is a pion when calculating the energy E, 
to be used in the formula above, from the (assumed known) momentum p. 

= 19 : pseudorapidity r? = ^ln{{p + pz)/{p — Pz))- 

= 20 : momentum fraction Xp = 2pjW, where W is the total energy of initial parton/particle 
configuration. 

= 21 : Xp = 2pz/W (Feynman-re if system is studied in CM frame). 

= 22 : rE_L = 2p±/W. 

= 23 : XE = 2E/W. 

= 24: •. z+ = {E + Pz)/W. 

= 25 : Z- = {E — Pz)/W. 


SUBROUTINE VNIANAL(MSEL) 

Purpose: to administrate in each event, which of the available analysis duties (described below) are carried 
out, and when during the time evolution of the particle system data samples are taken and accumula- ted 
to perform a statistical analysis at the end. It calls the routines VNIANAl - VNIANA5 to do this job, for 
the options that are switched on by default or by the user. The switches and their effect are described in 
Sec. 3.7. After the last event of a run, the results are automatically printed out in tabular form to files 
called VNI???.DAT, with numbers ???, corresponding to the selected analysis options. 

MS EL : gives the type of action to be taken. 

= 0 ; all arrays for data analysis in VNIANAl - VNIANA5 are reset. 

= 1 : data samples are accumulated, stored and analyzed in VNIANAl - VNIANA5. 

= 10 : final results are written to output files called VNI.777.DAT, where ‘777’ is a 3-digit number 
specihed in the descriptions below. 
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SUBROUTINE VNIANAl(MSEL) 


Purpose: to provide a number of global analyses referring to time- independent quantities of events as- 
a-whole, as e.g. total charged multiplicities, etc.. Is adopted from JETSET/PYTHIA and only slightly 
altered. Accumulated statistics is written out at the end of a run. When errors are quoted, these refer 
to the uncertainty in the average value for the event sample as a whole, rather than to the spread of the 
individual events, i.e. errors decrease like one over the square root of the number of events analyzed. 

MSEL : determines which action is to be taken. Generally, a last digit equal to 0 indicates that the 
statistics counters for this option is to be reset. Last digit 1 leads to an analysis of current event with 
respect to the desired properties. The statistics accumulated is output in tabular form with last digit 
2 , and written to disc. 

= 10 : statistics on parton multiplicity is reset. 

= 11 : the parton content of the current event is analyzed, classified according to the flavour content 
of the hard interaction and the total number of partons. 

= 12 : gives a table on parton multiplicity distribution, written to file VNI110.DAT. 

= 20 : statistics on particle content is reset. 

= 21 : the particle/parton content of the current event is analyzed, also for particles which have 
subsequently decayed and partons which have fragmented. Particles are subdivided into primary and 
secondary ones, the main principle being that primary particles are those produced directly, while 
secondary come from decay of other particles. 

= 22 : gives a table of particle content in events, written to file VNI120.DAT. 

= 30 : statistics on factorial moments is reset. 

= 31 : analyzes the factorial moments of the multiplicity distribution in different bins of rapidity and 
azimuth. 

= 32 : prints a table of the first four factorial moments for various bins of pseudorapidity and azimuth, 
written to file VNI130.DAT. The moments are properly normalized so that they would be unity (up 
to statistical fluctuations) for uniform and uncorrelated particle production according to Poissonian 
statistics, but increasing for decreasing bin size in case of ’’intermittent” behaviour. 

= 40 : statistics on energy-energy correlation is reset. 

= 41 : the energy-energy correlation (EEC) of the current event is analyzed. Events are assumed 
given in their CM frame. The weight assigned to a pair i and j is 2EiEj/W'^, where W is the sum of 
energies of all analyzed particles in the event. Energies are determined from the momenta of particles 
and masses. Statistics is accumulated for the relative angle 6ij, ranging between 0 and 180 degrees, 
subdivided into 50 bins. 

= 42 : prints a table of the energy-energy correlation (EEC) and its asymmetry (EECA), with er¬ 
rors, written to file VNI140.DAT. The definition of errors is not unique. Each event is viewed as one 
observation, i.e. an EEC and EECA distribution is obtained by summing over all particle pairs of an 
event, and then the average and spread of this event-distribution is calculated in the standard fashion, 
with the quoted error containing a factor of 1 /Vno. events. 

= 50 : statistics on complete final states is reset. 

= 51 : analyzes the particle content of the final state of the current event record. During the course 
of the run, statistics is thus accumulated on how often different final states appear. Only final states 
with up to 8 particles are analyzed, and there is only reserved space for up to 200 different final states. 
Most high energy events have multiplicities far above 8 , so the main use for this tool is to study the 
effective branching ratios obtained with a given decay model for e.g. charm or bottom hadrons. 

= 52 : gives a list of the (at most 200) channels with up to 8 particles in the final state, with their 
relative branching ratio, written to file VNI150.DAT. The ordering is according to multiplicity, and 
within each multiplicity according to an ascending order of KF codes. The KF codes of the particles 
belonging to a given channel are given in descending order. 
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SUBROUTINE VNIANA2(MSEL) 

Purpose: to analyze final state hadron distributions and calculate global distributions for the total (charged) 
ensemble of hadronic final states. 

MSEL : determines which action is to be taken, with the convention noted before. 

= 10 : statistics on y = ln{l/x) distribution of final hadrons is reset, where x = ‘lEj^fs is the 
energy fraction of a particle. 

= 11 : hadron content at the end of each event is binned in y = ln(l/a;) individually for all, and for 
charged, hadrons. Covered range in y is given by PARW(161) and PARW(162). 

= 12 : gives a table of the obtained y = ln{l/x) spectra, written to file VNI210.DAT. 

= 20 : statistics on Bose-Einstein correlations of final state pions and kaons is reset. 

= 21 : Bose-Einstein effects are simulated according to a simple parametrization, with an algorithm 
adopted from JETSET/PYTHIA. The distribution of pair masses is stored for both cases, without and 
with Bose-Einstein correlations taken into account. Covered range in pair mass is given by PARW(163) 
and PARW(164). 

= 22 ; gives a table of the obtained pair mass spectrum of ‘like-sign’ pions and kaons, without and 
with Bose-Einstein effects, written to file VNI220.DAT. 


SUBROUTINE VNIANA3(MSEL) 

Purpose: to analyze the dynamics of pre-hadronic clusters formed from partons during the parton-hadron 
conversion stage of the system. The cluster size- and mass-distributions, and momentum and spatial spectra 
are accumulated at various time steps during the time evolution of the system, so that the time change of 
the spectra gives a history of the cluster formation processes, locally in phase-space. 

MSEL : determines which action is to be taken, with the convention noted before. The statistics accu¬ 
mulated is output in tabular form with bins, ideally covering the range of momenta and coordinates 
of relevant particles. Depending on the kinematics, however, the limits of this range may require an 
adjustment, and/or a finer resolution of bins may be necessary. 

= 10 : statistics on cluster sizes and their invariant mass spectrum is reset. 

= 11 ; the cluster content of the system is analyzed each single time step during an an event, and 
cluster sizes and masses are stored in bins. Covered range of analysis is 0 < i? < 1.5 fm, and 0 < M < 7 
GeV. 

= 12 : gives a table on the obtained cluster size- and mass-spectra, written to file VNI310.DAT. 

= 20 : statistics on rapidity (y) spectra of clusters is reset. 

= 21 : the cluster content at certain pre-selected time steps during an an event is analyzed. Covered 
range in y is given by PARW(171) and PARW(172). 

= 22 : gives a table on the obtained y-spectra, written to file VNI320.DAT. 

= 30 : statistics on longitudinal (r^) spatial distribitions of clusters is reset. 

= 31 : the cluster positions at certain pre-selected time steps during an event are analyzed, with 
respect to the overall CM frame and to the z-axis, usually the jet-axis, or beam-axis. Covered range 
in rz is given by PARW(173) and PARW(174). 

= 32 : gives a table on the obtained r^-spectra, written to file VNI330.DAT. 

= 40 : statistics on transverse momentum {p±) distribitions of clusters is reset. 

= 41 : the cluster momenta p± at certain pre-selected time steps during an an event are analyzed, 
transverse with respect to the z-axis. Covered range in p± is given by PARW(175) and PARW(176). 
= 42 : gives a table on the obtained y_L-spectra, written to file VNI340.DAT. 
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= 50 : statistics on transverse (rj,) spatial distribitions of clusters is reset. 

= 51 : the cluster content at certain pre-selected time steps during an an event is analyzed, with 
respect to the distance transverse to the z-axis. Covered range in r± is given by PARW(177) and 
PARW(178). 

= 52 : gives a table on the obtained r_i_-spectra, written to file VNI350.DAT. 


SUBROUTINE VNIANA4(MSEL) 

Purpose: to analyze momentum and spatial distributions of partons and hadrons at various time steps during 
the time evolution of the system, so that the time change of the spectra gives a history of the dynamics 
from the initial state all the way to the final state of hadrons. 

MS EL : determines which action is to be taken, with the convention noted before. The statistics accu¬ 
mulated is output in tabular form with bins, ideally covering the range of momenta and coordinates 
of relevant particles. Depending on the kinematics, however, the limits of this range may require an 
extension. 

= 10 : statistics on rapidity (y) spectra of partons and hadrons is reset. 

= 11 ; the particle content at certain pre-selected time steps during an an event is analyzed, for gluons, 
quarks, mesons and baryons seperately. Covered range in y is given by PARW(181) and PARW(182). 
= 12 : gives a table on the obtained y-spectra, written to file VNI410.DAT. 

= 20 : statistics on longitudinal (r^) spatial distribitions of partons and hadrons is reset. 

= 21 ; the particle positions at certain pre-selected time steps during an event are analyzed, with 
respect to the overall CM frame and to the z-axis, usually the jet-axis, or beam-axis. Covered range 
in rz is given by PARW(183) and PARW(184). 

= 22 : gives a table on the obtained r^-spectra, written to file VNI420.DAT. 

= 30 : statistics on transverse momentum {p±) distribitions of partons and hadrons is reset. 

= 31 : the particle momenta p± at certain pre-selected time steps during an an event are analyzed, 
transverse with respect to the z-axis. Covered range in p± is given by PARW(185) and PARW(186). 
= 32 : gives a table on the obtained pT-spectra, written to file VNI430.DAT. 

= 40 : statistics on transverse (r_i_) spatial distribitions of partons and hadrons is reset. 

= 41 : the particle content at certain pre-selected time steps during an an event is analyzed, with 
respect to the distance transverse to the z-axis. Covered range in r± is given by PARW(187) and 
PARW(188). 

= 42 : gives a table on the obtained r_L-spectra, written to file VNI440.DAT. 


SUBROUTINE VNIANA5(MSEL) 

Purpose: to analyze the time evolution of global quantities and ‘thermodynamic’ properties of particles at 
various time steps during the time evolution of the system. 

MSEL : determines which action is to be taken, with the convention noted before. 

= 10 : statistics on time evolution of particle numbers, energy partitions, longitudinal momentum 
fractions and average transverse momentum is reset. 

= 11 : the particle content is analyzed in each single time step during an event, for the different 
particle species individually. 

= 12 : gives a table on the obtained time evolution profiles, written to file VNI510.DAT. 

= 20 : statistics on particle flow-velocity profiles is reset. 


72 



= 21 ; at certain pre-assigned time steps, the phase-space densities of the different particle species are 
analyzed, from this the local flow-velocities are calculated in cellular space, and stored in 2 -dimensional 
bins in z and r±. 

= 22 : gives a table on the obtained time evolution profiles, written to file VNI520.DAT. 

= 30 : statistics on particle densities, energy densities, and pressures of the various particle species 
is reset. 

= 31 : at certain pre-assigned time steps, the phase-space densities the different particle species are 
analyzed, convoluted with the previously obtained flow-velocity profile, and Lorentz invariant parti¬ 
cle densities, energy densities, diagonal and off-diagonal components of the pressures are stored in 
2-dimensional bins in z and r_\_. IMPORTANT; In order to use this option, the flow velocity profiles 
must! be obtained beforehand in a seperate run, with the same parameters and conditions, because 
they are needed to calculate the densities, and are read in as input data. 

= 32 : gives a table on the obtained density evolution profiles, written to file VNI530.DAT. 


SUBROUTINE VNISPHE(SPH,APL) 

Purpose: to diagonalize the momentum tensor, i.e. find the eigenvalues Ai > A 2 > A 3 , with sum unity, and 
the corresponding eigenvectors. Momentum power dependence is given by PARU(41); default corresponds 
to sphericity, PARU(41)=1. gives measures linear in momenta. Which particles (or partons) are used in the 
analysis is determined by the MSTU(41) value. 

SPH : 3 (A 2 -|-A 3 )/ 2 , i.e. sphericity (for PARU(41)=2.). Is returned =-1, if analysis not performed because 
event contained less than two particles (or two exactly back-to-back particles, in which case the two 
transverse directions would be undefined). 

APL : 3 / 2 A 3 , i.e. aplanarity (for PARU(41)=2.). Is returned = -1, as SPH = —1 

Remark: the lines N + I through A^ -|- 3 {N — 2 through N for MSTU(43) = 2) in VNIREC will, after a 
call, contain the following information: 

K(A^ -|- i, 1) = 31; 

K(A^-Pi,2) = 95; 

K(A^ + 3) : i, the axis number, i = 1, 2,3; 

K{N + i,4), K(iV + i,5) =0; 

P(A^ + A 1) - P(A' -P i, 3) : the i-th eigenvector, x, y and z components; 

P(A^ + ^4:) : Aj, the i-th eigenvalue; 

P(iV + i,5) = 0; 

V(A^ + i,l) - V(Ai + i,5) = 0. 

Also, the number of particles used in the analysis is given in MSTU(62). 


SUBROUTINE VNITHRU(THR,OBL) 

Purpose: to find the thrust, major and minor axes and corresponding projected momentum quantities, in 
particular thrust and oblateness. The performance of the program is affected by MSTU(44), MSTU(45), 
PARU(42) and PARU(48). In particular, PARU(42) gives the momentum dependence, with the default value 
1. corresponding to linear dependence. Which particles (or partons) are used in the analysis is determined 
by the MSTU(41) value. 

THR : thrust (for PARU(42)=1.). Is returned = -1, if analysis cannot be not performed, because event 
contained less than two particles. Is returned = -2, if remaining space in VNIREC (partly used as 
working area) is not large enough to allow analysis. 

OBL : oblateness (for PARU(42)=1.). Is returned = -1, respectively = -2, as for THR. 
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Remark: the lines + 1 through + 3 {N — 2 through N for MSTU(43) = 2) in VNIREC will, after a 
call, contain the following information: 

K(iV + i, 1) = 31; 

K(A^ + i,2) = 96; 

K(A^ + 3) : i, the axis number, i = 1, 2,3; 

K(iV + i,4), K(Af + l,5) =0; 

P(A^ + A 1) - P(A^ + b 3) : the thrust, major and minor axis, respectively, for i = 1, 2 and 3; 

P(A^ + i,4:) : corresponding thrust, major and minor value; 

P(Ar + i,5) =0; 

V(A^ + i,l) - V(Ar + z,5) = 0. 

Also, the number of particles used in the analysis is given in MSTU(62). 


SUBROUTINE VMICLUS(MJET) 

Purpose: to reconstruct an arbitrary number of jets using a cluster analysis method based on particle 
momenta. Two different distance measures are available, the traditional one roughly corresponding to 
relative transverse momentum, and one based on the JADE method, which roughly corresponds to invariant 
mass (in both cases with some important modifications). The choice is controlled by MSTU(46). The 
distance scale djoin, above which two clusters may not be joined, is normally given by PARU(44). In general, 
djoin may be varied to describe different ”jet resolution powers”. With the mass distance, PARU(44) can 
be used to set the absolute maximum cluster mass, or PARU(45) to set the scaled one, i.e. in y = m? jW"^, 
where IT^ is the total invariant mass-squared of the particles being considered. It is possible to continue the 
cluster search from the configuration already found, with a new higher djoin scale, by selecting MSTU(48) 
properly. In MSTU(47) one can also require a minimum number of jets to be reconstructed; combined with 
an artificially large djoin this can be used to reconstruct a predetermined number of jets. Which particles 
(or partons) are used in the analysis is determined by the MSTU(41) value, whereas assumptions about 
particle masses are given by MSTU(42). The parameters PARU(43) and PARU(48) regulate more technical 
details (for events at high energies and large multiplicities, however, the choice of a larger PARU(43) may 
be necessary to obtain reasonable reconstruction times). 

N JET : the number of clusters reconstructed. Is returned = -1, if analysis cannot performed, because event 
contained less than MSTU(47) (normally 1) particles, or analysis failed to reconstruct the requested 
number of jets. Is returned = -2, if remaining space in VNIREC (partly used as working area) is not 
large enough to allow analysis. 

Remark: if the analysis does not fail, further information is found in MSTU(61) - MSTU(63) and PARU(61) 
- PARU(63). In particular, PARU(61) contains the invariant mass for the system analyzed, i.e. the 
number used in determining the denominator oiy = re? . PARU(62) gives the generalized thrust, 
i.e. the sum of (absolute values of) cluster momenta divided by the sum of particle momenta (roughly 
the same as multicity). PARU(63) gives the minimum distance d (in p_i_ or m) between two clusters 
in the final cluster configuration, 0 in case of only one cluster. Further, the lines N + 1 through 
N + N JET (N — N JET + 1 through N for MSTU(43) = 2) in VNIREC will, after a call, contain the 
following information; 

K(V -T f, 1) = 31; 

K(iV + f,2) = 97; 

K(N + f, 3) : i, the jet number, with the jets arranged in falling order of absolute momentum; 

K{N + f, 4) : the number of particles assigned to jet i; 

K(V + f, 5) = 0; 

P(V + i, 1) - P(V + f, 5) : momentum, energy and invariant mass of jet i; 

V(V + f,l) - V(iV + z,5) = 0. 

Also, for a particle which was used in the analysis, K(/,4) = i, where I is the particle number and i 
the number of the jet it has ben assigned to. Undecayed particles not used then have K(1,4) = 0. An 
exception is made for lines with K(/, I) = 3 (which anyhow are not normally interesting for cluster 
search), where the color flow information stored in K(/, 4) is left intact. 
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SUBROUTINE VNICELL(NJET) 


Purpose: to provide a simpler cluster routine more in line with what is currently used in the study of high-p_i_ 
collider events. A detector is assumed to stretch in pseudorapidity between -PARU(51) and +PARU(51) 
and be segmented in MSTU(51) equally large rj (pseudorapidity) bins and MSTU(52) (p (azimuthal) bins. 
Transverse energy E± for undecayed entries are summed up in each bin. For MSTU(53) nonzero, the energy 
is smeared by calorimetric resolution effects, cell by cell. This is done according to a Gaussian distribution; 
if MSTU(53) = 1 the standard deviation for the E± is PARU(55) -v/ E±, if MSTU(53) = 2 the standard 
deviation for the E is PARU(55)\/^, E_i and E expressed in GeV. The Gaussian is cut off at 0 and at 
a factor PARU(56) times the correct E± or E. All bins with E± >PARU(52) are taken to be possible 
initiators of jets, and are tried in falling E± sequence to check whether the total E± summed over cells no 
more distant than PARU(54) in a/(A r/)^ + (A(/))^ exceeds PARU(53). If so, these cells dehne one jet, and 
are removed from further consideration. Gontrary to VNIGLUS, not all particles need be assigned to jets. 
Which particles (or partons) are used in the analysis is determined by the MSTU(41) value. 

NJET : the number of jets reconstructed (may be 0). Is returned = -2, if remaining space in VNIREG 
(partly used as working area) is not large enough to allow analysis. 

Remark: the lines A" + 1 through N + NJET {N — NJET + 1 through N for MSTU(43) = 2) in VNIREG 
will, after a call, contain the following information: 

K(iV + i, 1) = 31; 

K(A + i,2) = 98; 

K(A + i, 3) : i, the jet number, with the jets arranged in falling order in E_\_; 

K(A + i, 4) ; the number of particles assigned to jet i; 

K(A + i, 5) = 0; 

V(A + i,l) - V(A + f,5) = 0; 

Further, for MSTU(54) = 1: 

P(A + i, 1), P(A + f, 2) = position in rj and p of the center of the jet initiator cell, i.e. geometrical 
center of jet; 

P(A + f, 3), P(N+i,4) = position in rj and p of the -weighted center of the jet, i.e. the center of 
gravity of the jet; 

P(A + f, 5) = X) of the jet; 

while for MSTU(54) = 2; 

P(A + f, 1) - P(A + i, 5) ; the jet momentum vector, constructed from the summed E_i_ and the r/ and 
4> of the E_i_-weighted center of the jet as {px-,Py,Pz-,E,m) = E_i_(cos(/>, sin i;/), sinhr/, cosh 77 , 0 ); 

and for MSTU(54) = 3; 

P(A + i, 1) - P(A + i, 5) : the jet momentum vector, constructed by adding vectorially the momentum 
of each cell assigned to the jet, assuming that all the E±_ was deposited at the center of the cell, and 
with the jet mass in P(A + i, 5) calculated from the summed E and p as = E'^ — p'^ — Py — p^. 

Also, the number of particles used in the analysis is given in MSTU(62), and the number of cells 
hit in MSTU(63). 


SUBROUTINE VNIJMAS(PMH,PML) 

Purpose: to reconstruct high and low jet mass of an event, by dividing the collision event into two hemi¬ 
spheres, transversely to the sphericity axis. Then one particle at a time is reassigned to the other hemisphere 
if that reduces the sum of the two jet masses- squared rnjj + m\. The procedure is stopped when no further 
significant change (see PARU(48)) is obtained. Often, the original assignment is retained as it is. Which 
particles (or partons) are used in the analysis is determined by the MSTU(41) value, whereas assumptions 
about particle masses are given by MSTU(42). 
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PMH : heavy jet mass (in GeV). Is returned = -2m, if remaining space in VNIREC (partly used as 
working area) is not large enough to allow analysis. 

PML : light jet mass (in GeV). Is returned = -2, if remaining space in VNIREC (partly used as working 
area) is not large enough to allow analysis. 

Remark: light jet mass (in GeV). After a successful call, MSTU(62) contains the number of particles used 
in the analysis, and PARU(61) the invariant mass of the system analyzed. The latter number is helpful 
in constructing scaled jet masses. 


SUBROUTINE VMIFOWO(HIO,H20,H30,H40) 

Purpose: to do an event analysis in terms of the Eox-Wolfram moments The moments Hi are normalized to 
the lowest one, Hq. Which particles (or partons) are used in the analysis is determined by the MSTU(41) 
value. 

HIO : Hi/Hq. Is = 0 if momentum is balanced. 

H20 : H 2 /H 0 . 

H30 : H 3 /H 0 . 

H40 : 

Remark: the number of particles used in the analysis is given in MSTU(62). 
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C Further commonblocks 


The subroutines and commonblocks that one normally comes in direct contact with have already been de¬ 
scribed. For the sake of completeness, two further commonblocks are here described, which might be of 
interest to a user that desires a more detailed insight. The common block VNIINT serves as a tempo¬ 
rary storage of quantities during the simulation of a specific parton collision plus its associated branching 
evolution. It is very similar to a corresponding array in the JETSET/PYTHIA program. 

COMMON/VNIINT/MINT(400),VINT(400) 


Purpose: to provide an intermediate storage for integer and real valued variables as work area, used in the 
program to carry out the evolution of a single partonic subsystem at a time, consisting of a parton collision 
and its associated space-like (initial state) and time-like (final state) branchings. These variables must not 
be changed by the user. 


MINT(l) : specifies the general type of partonic subprocess that has occured, according to the ISUB code 
given in Table 3. 

MINT(2) : whenever MINT(l) (together with MINT(15) and MINT(16)) are not enough to specify the 
type of a 2 —> 2 parton collision uniquely, MINT(2) provides an ordering of the different possibilities, 
related to different color-flow topologies. The different possibilities are discussed in Refs. |5^. With i 
= MINT(15), j = MINT(16) and k = MINT(2), they are: 


ISUB = 1, i = j, qiQi qiqf, 
k = 1 : color configuration A. 
k = 2 : color configuration B. 

ISUB = 1, z = j, qiqj^ q^qj; 
k = 1 : only possibility. 

ISUB = 2, qtqt qiqp 
k = 1 : only possibility. 

ISUB = 3, qiqi gg] 

k = 1 : color configuration A. 
k = 2 : color configuration B. 

ISUB = 6, qig qig; 

k = 1 : color configuration A. 
k = 2 : color configuration B. 

ISUB = 8 , gg ^ qiqr, 

k = 1 : color configuration A. 
k = 2 : color configuration B. 

ISUB = 9 , gg ^ gg- 

k = 1 : color configuration A. 
k = 2 : colour configuration B. 
k = 3 : color configuration C. 


MINT(3) : number of partons produced in 
elements used, where currently only n = 


the hard interactions, i.e. the number n of the 2 
1 or n = 2 processes are included. 


n matrix 


MINT(4) : number of documentation lines in common block VNIREC that are reserved and used as work 
space for simulating a parton collision. 


MINT(5) 

MINT(6) 

MINT(7) 


not used. 

current frame of event as initialized, cf. MSTW(4). 

line number for documentation of outgoing particle 1 from a parton collision 2 —> 2 or 2 —> 1. 
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MINT(8) ; line number for documentation of outgoing particle 2 from a parton collision 2^2. Is 0 for 
2^1. 

MINT(9) : not used. 

MINT(IO) : not used. 

MINT(ll) : KF flavour code of initial beam (side 1) particle (c.f. MSTW(5)). 

MINT(12) : KF flavour code of initial target (side 1) particle (c.f. MSTW(6)). 

MINT(13) : KF flavour code of side 1 parton that initiates a space-like cascade (only for primary partons, 
is 0 else). 

MINT(14) : KF flavour code of side 2 parton that initiates a space-like cascade (only for primary partons, 
is 0 else). 

MINT(15) : KF flavour code for side 1 incoming parton to collision subprocess. 

MINT(16) : KF flavour code for side 2 incoming parton to collision subprocess. 

MINT(17) - MINT(20) : not used. 

MINT(21) : KF flavour code for parton 3 outgoing from a parton collision. 

MINT(22) : KF flavour code for parton 4 outgoing from a parton collision (is 0 for 2^1 processes). 

MINT(23) : KF flavour code for parton 5 outgoing from a parton collision - not used currently (relevant 
only for 2 ^ 3, 2 ^ 4 processes). 

MINT(24) : KF flavour code for parton 6 outgoing from a parton collision - not used currently (relevant 
only for 2^4 processes). 

MINT(25) - MINT(30) : not used. 

MINT(31) - MINT(40) : not used. 

MINT(41) : type of initial beam particle; 1 for lepton and 2 for hadron and 3 for nucleus (c.f. MSTW(7)). 

MINT(42) : type of initial target particle; 1 for lepton and 2 for hadron and 3 for nucleus (c.f. MSTW(8)). 

MINT(43) : combination of incoming beam and target particles (c.f. MSTW(9)). 

MINT(44) : as MINT(43), but a photon counts as a lepton. 

MINT(45), MINT(46) : structure of incoming beam and target particles. 

= 1 : no internal structure, i.e. lepton carries full beam energy. 

= 2 : defined with structure functions, i.e. hadrons and nuclei. 

MINT(47) : combination of incoming beam and target particle structure function types. 

= 1 ; no structure function either for beam or target. 

= 2 : structure functions for target but not for beam. 

= 3 ; structure functions for beam but not for target. 

= 4 ; structure functions both for beam and target. 

MINT(48) : total number of subprocesses switched on. 
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MINT(49), MINT(50) : not used. 

MINT(51) : internal flag that the simulation of a 2 -parton collision system, including space-like initial 
state and time-like hnal state cascades, has failed cuts. 

= 0 ; no problem. 

= 1 ; failed and to be skipped. 

MINT(52) - MINT(60) : not used. 

MINT(61) - MINT(70) : not used. 

MINT(81), MINT(82) : not used. 

MINT(83) : number of lines in the event record already filled by previously considered pileup events. 

MINT(84) : MINT(83) -|- MSTP(126), i.e. number of lines already filled by previously considered events 
plus number of lines to be kept free for event documentation. 

VINT(l) ; Ecu, CM energy. 

VINT(2) ; s = Eq]^, mass-square of complete system. 

VINT(3) ; mass of beam particle. 

VINT(4) ; mass of target particle. 

VINT(5) ; momentum of beam (and target) particle in CM frame. 

VINT(6) ; angle 0 for rotation from CM frame to fixed target frame. 

VINT(7) ; angle ^ for rotation from CM frame to fixed target frame. 

VINT(8) ; [3x veiocity for boost from CM frame to fixed target frame. 

VINT(9) ; [3y velocity for boost from CM frame to fixed target frame. 

VINT(IO) : j3z velocity for boost from CM frame to fixed target frame. 

VINT(ll) : (3x velocity for boost from overall CM frame to ‘microscopic’ cm frame of a parton-parton 
collision subsystem. 

VINT(12) : Py velocity for boost from overall CM frame to parton-parton cm frame. 

VINT(13) : Pz velocity for boost from overall CM frame to parton-parton cm frame. 

VINT(14) : angle 6 for rotation from overall CM frame to ‘microscopic’ cm frame of a parton-parton 

collision subsystem. 

VINT(15) : angle p for rotation from overall CM frame to parton-parton cm frame. 

VINT(16) : Lorentz boost factor 7 to parton-parton collision cm frame. 

VINT(17) : Rapidity ycm of 2-parton cm sytem with respect to global CM frame. 

VINT(18) : Transverse momentum p_\_ cm of 2-parton cm sytem in global CM frame. 

VINT(19) : Longitudinal separation Af^ of 2 partons in their c.m. frame of collision, at the point of 
‘closest approach’. 
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VINT(20) : Transverse separation Ar± of 2 partons in their c.m. frame of collision, at the point of ‘closest 
approach’. 

VINT(21) : r = X 1 X 2 S with xi ^2 the longitudinal momentum fractions (Bjorken-x) of 2 colliding partons. 
VINT(22) : y* = 1/2 ln(xi/x 2 ), the ‘relative’ rapidity of 2 colliding partons. 

VINT(23) : cos( 0 )) (scattering angle) in 2 ^ 2 parton collision. 

VINT(24) : (j) (azimuthal angle) in in 2 ^ 2 parton collision. 

VINT(25) - VINT(30) : not used. 

VINT(31) : ti, lower value of kinematic limit on t in parton collision. 

VINT(32) : upper value of kinematic limit on t in parton collision. 

VINT(33) : |t|i, upper absolute value of limit of integration over dt of 2 -parton cross-section. 
VINT(34) : |t| 2 , lower absolute value of limit of integration over dt of 2 -parton cross-section. 

VINT(35) : cos{e)max for 0 < 0. 

VINT(36) : cos{9)max for 0 > 0. 

VINT(37) : integrated 2 -parton cross-section d'(s) = / dp‘]_{dd{s,p‘j_)/dp‘j_. 

VINT(38) : differential 2 -parton cross-section dd{s,p‘]_)/dp‘]_. 

VINT(39) : characteristic proper collision time Tcol{p\) oc l/p± in the 2 -body c.m.-frame. 

VINT(40), not used. 

VINT(41) : the momentum fraction xi taken by parton 1 at the 2 -parton vertex. 

VINT(42) : the momentum fraction X 2 taken by parton 2 at the 2 -parton vertex. 

VINT(43) : m = ^/I, mass of 2 -parton collision subsystem. 

VINT(44) : s of the collision subprocess (2 ^ 2 , 2 ^ 1 ). 

VINT(45) : i of the collision subprocess (2 ^ 2 , 2 ^ 1 ). 

VINT(46) : u of the collision subprocess (2 ^ 2 , 2 ^ 1). 

VINT(47) : p± of the collision subprocess, i.e., transverse momentum evaluated in the rest frame of the 
collision. 

VINT(48) : p\ of the collision subprocess; see VINT(47). 

VINT(49) : p±min of process, i.e. PARV(42), or max(VKIN(3),VKIN(5)), depending on which is larger, 
and whether the process is singular in p± 0 , or not. 

VINT(50) : not used. 

VINT(51) : Q of the hard subprocess. The precise definition is process-dependent, see MSTV(43). 
VINT(52) : Q'^ of the hard subprocess; see VINT(51). 

VINT(61) - VINT(65) : 9, ip and (3 for rotation and boost from c.m. frame to hadronic c.m. frame of 
a lepton-hadron or lepton-nucleus event. 
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CDMMDN/VNICDL/MCDL(1000,0:10),VCOL(1000,0:10) 


Purpose: to provide an intermediate storage for a minimal set of integer and real valued variables for all 
detected parton collisions that occur during a single timestep during an event. The first index I in the 
arrays MCOL and VCOL labels the sequence number of a collision and the second index J is a pointer to 
interface with the above MINT and VINT arrays wherein the calculated quantities are written to. 


MCOL(I,l) : 

line number in 

MCOL(I,2) : 

line number in 

MCOL(I,3) : 

= MINT(l) 

MCOL(I,4) : 

= MINT(2) 

MCOL(I,5) : 

= MINT(3) 

MCOL(I,6) : 

= MINT(15) 

MCOL(I,7) : 

= MINT(16) 

MCOL(I,8) : 

= MINT(21) 

MCOL(I,9) : 

= MINT(22) 

MCOL(I,10) 

not used. 

VCOL(I,l) : 

= VINT(23) 

VCOL(I,2) : 

= VINT(24) 

VCOL(I,3) : 

= VINT(26) 

VCOL(I,4) : 

= VINT(41) 

VCOL(I,5) : 

= VINT(42) 

VCOL(I,6) : 

= VINT (44) 

VCOL(I,7) : 

= VINT(45) 

VCOL(I,8) : 

= VINT(46) 

VCOL(I,9) : 

= VINT(48) 

VCOL(I,10): 

= VINT(52) 


VNIREC of parton 1 entering the 
VNIREC of parton 2 entering the 


collision 

collision 


subprocess. 

subprocess. 
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D Sample output 

1. General information on the simnlation as-a-whole 

The following output stems from the example program vnixpll.f, simulating 1000 pp collisions at CERN 
collider energy ^/s = 546 GeV. Unless switched of by setting MSTV(llA) = 0, this type of output informa¬ 
tion is provided automatically and written to the file VNIRUN.DAT at the end of a simulation. The output 
includes a summary of the initialization overall collision system, the gluon and quark content of the initial 
state, the type and number of occured interaction processes during the parton-cascade development, the 
number and type of produced pre-hadronic clusters from partons, and the numbers of resulting final-state 
hadrons. 
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Main references: 


VNIXRIN: initialization Of VNI rOUtiUeS ****************** 


1 1 

I VNI will be initialized for a I 

I I 

I p+ on p- collider I 

I I 

I at 546.00 GeV center-of-mass energy I 

I I 

I corresponding to 546.00 GeV cm-energy per hadron I 

I I 
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I 

h + h -> jets 

+ X -> hadrons 


I 

I 





I 

I 

Number of collision events: 

10000 

I 

I 

Number of time steps per event 

: 193 

I 

I 





I 

I 

Initial time: .00 1/GeV 

= .00 

fm = 

.OOE+00 sec 

I 

I 

Final time: 166.75 1/GeV 

= 32.90 

fm = 

.22E-22 sec 

I 

I 





I 

I 





I 

I 

Initial parton distributions from VNI 


I 

I 

nucleon structure functions GRV 

set 

1 

I 

I 

pion structure functions DO 

set 

1 

I 

I 





I 

I 

Minimum initial resolution scale Q_0: 


3.02 GeV 

I 

I 

Minimum sqrt(s) in parton collisions: 


2.00 GeV 

I 

I 

Minimum pT_0 in parton collisions: 


2.07 GeV 

I 

I 

Minimum q_0 in space-like branchings: 


3.02 GeV 

I 

I 

Minimum mu_0 in time-like branchings: 


1.48 GeV 

I 

I 





I 

I 





I 

I 

Initial Collision System: 



I 

I 





I 

I 


Beam 


Target 

I 

I 

Particle: 

P+ 


P- 

I 

I 





I 

I 

Number of neutrons: 

0 


0 

I 

I 

Number of protons: 

1 


-1 

I 

I 

Number of d_v quarks: 

1 


-1 

I 

I 

Number of u_v quarks: 

2 


-2 

I 

I 

Number of s_v quarks: 

0 


0 

I 

I 

Radius R_ms (fm): 

.810 


.810 

I 

I 





I 

I 

Velocity: 

.999994 


.999994 

I 

I 

Gamma factor: 

290.505 


290.505 

I 

I 

Rapidity: 

6.365 


-6.365 

I 

I 





I 

I 

p_x (GeV): 

.000 


.000 

I 

I 

p_y (GeV): 

.000 


.000 

I 

I 

p_z (GeV): 

272.998 


-272.998 

I 

I 

E (GeV): 

273.000 


273.000 

I 

I 

M (GeV): 

.938 


.938 

I 

I 





I 

I 

charge px_tot py_tot 

pz_tot 

E_tot 

M_tot 

I 

I 

.0 0. 0. 

0. 

546. 

546. 

I 

I 





I 

I 

Initial separation parallel 

to beam: Delta Z 

= .25 fm 

I 

I 

Min/max impact transverse to 

beam: .00 < B 

< 2.29 fm 

I 

I 





I 

I 





I 


:*!:*!******************** VNIXRIN: initialization completed ********************* 
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:t!*********************:i= VNIXFIN: run finished normally *********************** 


I I 
I Total CPU time needed: 7 hours 8 minutes 55 seconds I 
I I 


INITIAL STATE 


Internal parton structure functions selected: 

set 1(1) for nucleons (pions) at initial QO = 2.870 GeV 

Minimum resolvable x-value resulting from spread of confined 
partons with 1/xP < R_conf = .20 fm: 


Beam p+ : x_min = .366E-02 

Target p~- : x_min = .366E-02 


Average sum of momentum fractions per hadron (nucleon): 

Beam p+ : <sum x> = .9940 

Target p~- : <sum x> = .9958 


Average number of initial state partons: 


gluons 

46.7329 

u+u~+d+d~ 

17.7952 

u-u~+d-d~ 

16.0016 

s+s~ 

3.4325 

c+c~ 

.6430 

b+b~ 

.0000 

t+t~ 

.0000 

total 

87.6017 


PARTON CASCADE EVOLUTION 


Total number of events: 10000 

Total number of time steps: 193 

Initial/final time: .0 fm / 32.9 fm 

Analysis based on 10000 event samples 
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All numbers per non-diffractive event: fraction N_nd/N_tot = .3870 

Actual minimum momentum transfer for hard collisions: 2.0699 GeV 


hard 2 -> 2 collisions: 
soft 2 -> 2 collisions: 
fusions 2 -> 1: 
space-like branchings 1 
time-like branchings 1 


switched on 
switched off 
switched on 

> 2: switched on 

> 2: switched on 


Average number of parton collisions a+b->c (+d): 


hard 2 -> 2 


q+q, q+q~ 

. q~+q~ 

.0399 

q+g. 

q~+g 

.3146 


g +g 

.7230 

soft 2 -> 2 

q+q, q+q~ 

. q~+q~ 

.0000 

q+g. 

q~+g 

.0000 


g +g 

.0000 

fusion 2 -> 1 


q+q~ 

.0004 

q+g. 

q~+g 

.0000 


g +g 

.0343 

total a+b -> 

c(+d) 

1.122 


Average number of space-like branchings a -> b + c: 


( CMshower -> b+c : 1.6852 ) 

q/q~ -> q/q~+g : .2068 

g -> g+g : -5309 

g • .0226 

q/q~ -> q/q~+gamma: .0000 

total a -> b+c : .7603 


Average number of time-like branchings a -> b + c: 


CMshower -> b+c : 

4.3426 

q/q~ -> q/q~+g : 

.3128 

g -> g+g : 

2.0833 

g -> q+q~ : 

.1996 

q/q~ -> q/q~+gamma: 

.0007 

total a -> b+c : 

2.5964 
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PARTON CLUSTER FORMATION 


NOTE: numbers refer to clusters from interacted partons only; 

beam/target remnant clusters (if present) are not included. 

Total number of events: 10000 

Total number of time steps: 193 

Initial/final time: .0 fm / 32.9 fm 

Analysis based on 10000 event samples 

All numbers per non-diftractive event: fraction N_nd/N_tot = .3870 


Total numbers of produced clusters: 


all: 
fast: 
slow: 

exogamous: 


3.1235 

2.3974 ( > PARV(88)= 

.7261 ( < PARV(88)= 

.6583 


1.000 MSTV(88)= 0 ) 
1.000 MSTV(88)= 0 ) 


Number 

of elementary 

cluster 

formation 

processes 




type 






total 

exogamous 

g 

+ 

g 

-> 

c 

+ 

C 



.2244 

.0287 

g 

+ 

g 

-> 

c 

+ 

g 



.5929 

.0757 

g 

+ 

g 

-> 

c 

+ 

g 

+ 

g 

.9401 

.1436 

q 

+ 

q~ 

-> 

c 

+ 

C 



.0946 

.0357 

q 

+ 

q~ 

-> 

c 

+ 

g 



.1752 

.0646 

q 

+ 

q~ 

-> 

c 

+ 

g 

+ 

g 

.0000 

.0000 

g 

+ 

q 

-> 

c 

+ 

q 



.3629 

.0964 

g 

+ 

q 

-> 

c 

+ 

q 

+ 

g 

.7334 

.2136 

qq 4 

■ q 

-> 

c 

+ 

c 



.0000 

.0000 

sum: 








3.1235 

.6583 


Distribution of exogamous cluster processes: 


origin 

all 

fast 

slow 

0.00 - 

0.05 

.3509 

.2902 

.0607 

0.05 - 

0.20 

.0427 

.0262 

.0165 

0.20 - 

0.30 

.0587 

.0430 

.0157 

0.30 - 

0.45 

.0401 

.0247 

.0154 

0.45 - 

0.55 

.2006 

.1666 

.0340 

0.55 - 

0.70 

.0521 

.0304 

.0217 

0.70 - 

0.80 

.1275 

.0906 

.0369 

0.80 - 

0.95 

.1357 

.0731 

.0626 

0.95 - 

1.00 

2.1152 

1.6526 

.4626 


CLUSTER HADRONIZATION AND UNSTABLE PARTICLE DECAYS 


NOTE: numbers refer to both, clusters from interacted partons 
as well as beam/target remnant clusters (if present). 

Total number of events: 10000 
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Total number of time steps: 193 

Initial/final time: .0 fm / 32.9 fm 

Analysis based on 10000 event samples 

parton-cluster fragmentation: switched on 

beam/target-cluster fragmentation: switched on 


Average number of primary particles from cluster decays: 

type charged + neutral charged only 


leptons, gauge bosons: 

.0000 

.0000 

light mesons: 

15.2396 

7.3508 

strange mesons: 

4.5977 

1.4296 

charm, bottom mesons: 

.0442 

.0216 

tensor mesons: 

1.4084 

.7034 

light baryons: 

2.2769 

1.5013 

strange baryons: 

1.7426 

.3953 

charm, bottom baryons: 

.0004 

.0004 

other: 

.0012 

.0004 

total: 

40.1700 

11.4028 


Average number of secondary particles from unstable particle decays: 


type 

charged + neutral 

charged only 

leptons, gauge bosons: 

1.1555 

.0018 

light mesons: 

32.2698 

20.2790 

strange mesons: 

3.1257 

1.5487 

charm, bottom mesons: 

.0000 

.0000 

tensor mesons: 

.0000 

.0000 

light baryons: 

4.6264 

2.3977 

strange baryons: 

.0000 

.0000 

charm, bottom baryons: 

.0000 

.0000 

other: 

.0000 

.0000 

total: 

41.0902 

24.1830 
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2. A typical event listing 


The event listing below is part of the standard output obtained with the example program vnixpll.f 
for 1000 pp collisions at CERN collider energy ^/s = 546 GeV. The listing, which is is written to the file 
VNIRUN.DAT, embodies the complete history of the first event from the initial beam particles p and p, 
their original gluon and quark content, the newly produced partons by hard scattering and radiation during 
the parton-cascade development, the formed pre-hadronic clusters and resulting final state particles. The 
latter emerge from both the parton clusters that result from parton-cascade products, and the beam/target 
clusters that consist of the remnant initial partons which have not participated in the parton cascades. 
As explained in Sec. 3.6, each entry contains name, status, flavor, and origin of a particle, as well as it’s 
three-momentum, energy and mass. Particles in brackets are inactive, decayed or fragmented particles. The 
final line of the listing shows the total charge, three-momentum, anargy and invariant mass of the particle 
system, given by the sums over all active, stable particles. 


Particle listing (summary: momentum space P in GeV) 
event no. 1 time step 200 at time 34.75 fm 


1 

particle 

KS 

KF 

C 

A 

P_x 

'<! 

P_z 

E 

M 

1 

(p+) 

15 

2212 

0 

0 

.000 

.000 

272.998 

273.000 

.938 

2 

(p~-) 

15 

-2212 

0 

0 

.000 

.000 

-272.998 

273.000 

.938 

3 

(p+) 

15 

2212 

0 

0 

.000 

.000 

272.998 

273.000 

.938 

4 

(p--) 

15 

-2212 

0 

0 

.000 

.000 

-272.998 

273.000 

.938 

5 

(u) 

13 

2 

1 

0 

.650 

.048 

1.786 

2.319 

.006 

6 

(d) 

13 

1 

1 

0 

-.593 

-.165 

2.092 

2.587 

.010 

7 

(d-) 

13 

-1 

0 

3 

-.135 

-.452 

6.724 

6.886 

-.490 

8 

(d) 

13 

1 

3 

0 

.413 

-.211 

.408 

.935 

-.484 

9 

(g) 

14 

21 

1 

1 

1.751 

.556 

2.232 

3.398 

-.682 

10 

(u) 

13 

2 

2 

0 

.263 

-.405 

.751 

1.258 

-.507 

11 

(u-) 

13 

-2 

0 

1 

-.056 

-.437 

.408 

.935 

-.460 

12 

(d) 

13 

1 

1 

0 

-.034 

-.126 

.826 

1.328 

-.149 

13 

(d-) 

13 

-1 

0 

3 

.027 

.378 

.408 

.935 

-.359 

14 

(d-) 

13 

-1 

0 

3 

-.287 

.276 

4.047 

4.364 

-.376 

15 

(d) 

13 

1 

1 

0 

.017 

.556 

.408 

.935 

-.536 

16 

(u-) 

13 

-2 

0 

3 

.396 

.020 

7.817 

7.916 

-.409 

17 

(u) 

13 

2 

3 

0 

.076 

-.052 

.473 

.996 

-.114 

18 

(g) 

13 

21 

1 

1 

-.050 

.419 

.539 

1.058 

-.400 

19 

(g) 

13 

21 

1 

2 

.197 

-.140 

1.146 

1.630 

-.263 

20 

(d) 

13 

1 

2 

0 

-.585 

.022 

.751 

1.258 

-.574 

21 

(d-) 

13 

-1 

0 

1 

-.054 

.096 

.408 

.935 

-.087 

22 

(g) 

13 

21 

1 

2 

-.223 

-.221 

.903 

1.401 

-.322 

23 

(u-) 

13 

-2 

0 

2 

-.473 

-.482 

.408 

.935 

-.683 

24 

(u) 

13 

2 

1 

0 

-.260 

.767 

.408 

.935 

-.786 

25 

(u) 

13 

2 

1 

0 

.096 

.052 

2.199 

2.622 

-.112 

26 

(u-) 

13 

-2 

0 

2 

.178 

.298 

.408 

.935 

-.336 

27 

(g) 

13 

21 

1 

3 

.286 

.066 

1.231 

1.710 

-.301 

28 

(g) 

13 

21 

1 

2 

.434 

-.103 

1.886 

2.328 

-.463 

29 

(u-) 

13 

-2 

0 

3 

.131 

-.005 

.408 

.935 

-.145 

30 

(u) 

13 

2 

1 

0 

-.138 

-.218 

.408 

.935 

-.271 

31 

(s) 

13 

3 

2 

0 

-.510 

.046 

2.092 

2.522 

-.537 

32 

(s~) 

13 

-3 

0 

2 

-.105 

-.467 

.408 

.935 

-.535 

33 

(g) 

13 

21 

1 

3 

-.238 

.440 

2.645 

3.043 

-.477 

34 

(u) 

13 

2 

2 

0 

-.026 

.071 

8.787 

8.831 

-.062 

35 

(u-) 

13 

-2 

0 

2 

.274 

-.108 

.408 

.935 

-.313 

36 

(g) 

13 

21 

1 

2 

-.266 

-.039 

47.203 

45.031 

-.309 

37 

(u-) 

13 

-2 

0 

1 

-.189 

.247 

.826 

1.328 

-.288 

38 

(u) 

13 

2 

3 

0 

-.118 

.433 

.408 

.935 

-.426 
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39 (u~) 

13 

-2 

0 

3 

.409 

.237 

2.762 

3.153 

-.473 

40 (u) 

13 

2 

3 

0 

-.384 

-.308 

.408 

.935 

-.497 

41 (u) 

13 

2 

1 

0 

-.001 

.148 

.408 

.935 

-.128 

42 (u~) 

13 

-2 

0 

2 

.040 

-.094 

.408 

.935 

-.126 

43 (g) 

13 

21 

1 

3 

.451 

.038 

1.688 

2.141 

-.463 

44 (d~) 

13 

-1 

0 

3 

-.236 

.582 

1.593 

2.051 

-.605 

45 (d) 

13 

1 

2 

0 

-.026 

-.061 

.408 

.935 

-.084 

46 (g) 

13 

21 

3 

1 

.233 

-.474 

163.962 

155.058 

-.792 

47 (u~) 

13 

-2 

0 

2 

-.126 

.025 

-9.757 

8.645 

.006 

48 (d~) 

13 

-1 

0 

2 

.526 

-.150 

-16.659 

15.157 

.010 

49 (u~) 

13 

-2 

0 

3 

-.121 

-.114 

-4.941 

4.106 

-.174 

50 (u) 

13 

2 

3 

0 

-.252 

-.044 

-1.576 

.935 

-.249 

51 (g) 

13 

21 

3 

2 

-.110 

.229 

-12.733 

11.449 

-.234 

52 (d~) 

13 

-1 

0 

3 

.051 

-.270 

-3.475 

2.724 

-.298 

53 (d) 

13 

1 

1 

0 

.175 

.560 

-1.576 

.935 

-.571 

54 (d~) 

13 

-1 

0 

2 

-.357 

-.255 

-2.231 

1.552 

-.442 

55 (d) 

13 

1 

1 

0 

-.142 

-.121 

-1.576 

.935 

-.193 

56 (d~) 

13 

-1 

0 

2 

.527 

-.078 

-3.586 

2.828 

-.547 

57 (d) 

13 

1 

2 

0 

.288 

-.539 

-1.576 

.935 

-.636 

58 (d) 

13 

1 

1 

0 

.259 

-.102 

-1.994 

1.328 

-.297 

59 (d) 

14 

1 

2 

0 

-1.532 

-.540 

-1.937 

2.087 

-.108 

60 (d~) 

13 

-1 

0 

2 

.004 

-.395 

-1.576 

.935 

-.416 

61 (d) 

13 

1 

2 

0 

.139 

-.028 

-2.314 

1.630 

-.158 

62 (d~) 

13 

-1 

0 

3 

.020 

-.413 

-1.576 

.935 

-.435 

63 (u~) 

13 

-2 

0 

1 

.071 

-.222 

-4.675 

3.855 

-.257 

64 (u) 

13 

2 

3 

0 

.287 

-.066 

-1.576 

.935 

-.311 

65 (d~) 

13 

-1 

0 

3 

.115 

-.103 

-2.071 

1.401 

-.178 

66 (d) 

13 

1 

2 

0 

.173 

.196 

-1.576 

.935 

-.255 

67 (u) 

13 

2 

3 

0 

-.800 

-.123 

-1.994 

1.328 

-.802 

68 (u~) 

13 

-2 

0 

3 

-.434 

-.218 

-1.576 

.935 

-.486 

69 (g) 

13 

21 

1 

2 

.092 

-.434 

-13.900 

12.548 

-.469 

70 (d~) 

13 

-1 

0 

3 

-.200 

-.122 

-2.231 

1.552 

-.237 

71 (d) 

13 

1 

3 

0 

-.153 

-.233 

-1.576 

.935 

-.291 

72 (g) 

13 

21 

2 

1 

-.558 

.600 

-16.397 

14.901 

-.798 

73 (g) 

14 

21 

1 

1 

-.646 

-.141 

4.334 

4.734 

-.682 

74 (d) 

14 

1 

2 

0 

-.014 

.151 

-1.915 

1.286 

-.108 

75 !g! 

21 

21 

0 

0 

2.954 

3.873 

-.435 

5.206 

-.682 

76 !d! 

21 

1 

0 

0 

-.354 

-.234 

-2.970 

2.639 

-.108 

77 !g! 

21 

21 

2 

1 

2.181 

5.310 

-3.320 

6.631 

.000 

78 !d! 

21 

1 

1 

0 

.521 

-1.530 

-.381 

1.660 

.010 

79 (g) 

14 

21 

1 

1 

2.954 

3.873 

-.435 

5.206 

-.682 

80 (d) 

14 

1 

2 

0 

-.354 

-.234 

-2.970 

2.639 

-.108 

81 (g) 

14 

21 

2 

1 

2.160 

5.284 

-3.324 

6.486 

.000 

82 (d) 

14 

1 

1 

0 

.500 

-1.556 

-.386 

1.624 

.010 

83 (g) 

14 

21 

0 

0 

2.052 

2.858 

1.999 

1.045 

-3.909 

84 (d) 

14 

1 

0 

0 

.000 

.000 

-8.008 

8.008 

.000 

85 (d) 

14 

1 

2 

1 

.000 

.000 

8.008 

8.008 

.000 

86 (CMshower) 

11 

94 

0 

0 

2.661 

3.728 

-3.710 

8.110 

5.571 

87 (d) 

14 

1 

2 

0 

1.128 

.073 

-1.364 

3.504 

3.023 

88 (cluster) 

11 

91 

0 

0 

1.116 

1.615 

-2.687 

3.412 

.750 

89 (cluster) 

11 

91 

0 

0 

-.450 

-.605 

-.228 

1.213 

.922 

90 (cluster) 

11 

91 

0 

0 

-1.528 

-.662 

1.381 

2.181 

.280 

91 (cluster) 

11 

91 

0 

0 

1.252 

.689 

-1.272 

3.136 

2.484 

92 (cluster) 

11 

91 

0 

0 

.110 

-.650 

5.219 

5.796 

2.433 

93 piO 

6 

111 

0 

0 

.660 

.523 

-.880 

1.171 

.135 

94 piO 

6 

111 

0 

0 

.461 

1.090 

-2.094 

2.377 

.135 

95 rho- 

6 

-213 

0 

0 

-.334 

-.562 

-.332 

1.068 

.769 

96 piO 

6 

111 

0 

0 

-.110 

-.046 

-.182 

.186 

.135 

97 pi- 

6 

-211 

0 

0 

-.678 

-.312 

.548 

1.066 

.140 

98 pi+ 

6 

211 

0 

0 

-.845 

-.352 

.547 

1.203 

.140 

99 a 2- 

6 

-215 

0 

0 

.161 

.747 

-.904 

1.772 

1.318 
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100 

rho+ 

6 

213 

0 

0 

1 

.097 

- 

.060 


,654 

1, 

,490 


.769 

101 

rho+ 

6 

213 

0 

0 

1 

.139 

- 

.488 

2, 

,544 

3, 

,178 


.769 

102 

piO 

6 

111 

0 

0 

-1 

.023 

- 

.165 

2, 

,389 

2, 

,849 


.135 

103 

(cluster) 

11 

91 

0 

0 

- 

.054 

- 

.150 


,360 


,686 


.562 

104 

piO 

6 

111 

0 

0 


.146 

- 

.046 


,232 


,425 


.135 

105 

piO 

6 

111 

0 

0 

- 

.194 

- 

.107 


,157 


,261 


.135 

106 

(cluster) 

11 

91 

0 

0 


.729 


.554 


,798 

1, 

,246 


.280 

107 

piO 

6 

111 

0 

0 


.491 


.351 


,622 


,810 


.135 

108 

piO 

6 

111 

0 

0 


.244 


.200 


,462 


,467 


.135 

109 

(Beam-REM) 

11 

92 

0 

0 

- 

.315 


.737 

274, 

,428 

281, 

,871 

70 

.164 

110 

(Targ-REM) 

11 

93 

0 

0 

- 

.907 

-1 

.587 

-276, 

,508 

247, 

,044 

-116 

.181 

111 

(CMF) 

15 

100 

0 

0 

-1 

.222 

- 

.850 

-2, 

,080 

528, 

,915 

528 

.909 

112 

(cluster) 

11 

91 

0 

0 

- 

.315 


.737 

274, 

,428 

281, 

,871 

70 

.164 

113 

(cluster) 

11 

91 

0 

0 

- 

.907 

-1 

.587 

-276, 

,508 

247, 

,044 

-116 

.181 

114 

(Beam-REM) 

11 

92 

0 

0 


.411 


.648 

127, 

,485 

127, 

,543 

3 

.762 

115 

(a_l-) 

17 

-20213 

0 

0 


.493 


.228 

19, 

,200 

19, 

,250 

1 

.275 

116 

(Deltat) 

17 

2214 

0 

0 

- 

.081 


.420 

108, 

,286 

108, 

,293 

1 

.232 

117 

piO 

7 

111 

0 

0 


.490 


.124 

4, 

,473 

4, 

,831 


.135 

118 

(rho-) 

17 

-213 

0 

0 


.005 


.102 

14, 

,584 

14, 

,605 


.769 

119 

piO 

7 

111 

0 

0 

- 

.160 


.015 

35, 

,475 

37, 

,042 


.135 

120 

P+ 

7 

2212 

0 

0 


.085 


.402 

72, 

,525 

75, 

,579 


.938 

121 

piO 

7 

111 

0 

0 

- 

.122 

- 

.268 

8, 

,519 

9, 

,014 


.135 

122 

pi- 

7 

-211 

0 

0 


.133 


.367 

5, 

,779 

6, 

,174 


.140 

123 

(Targ-REM) 

11 

93 

0 

0 

-1 

.252 

- 

.510 

-63, 

,766 

63, 

,825 

2 

.384 

124 

(Delta~-) 

17 

-2214 

0 

0 

- 

.779 

- 

.303 

-21, 

,201 

21, 

,253 

1 

.232 

125 

(eta) 

17 

221 

0 

0 

- 

.473 

- 

.208 

-42, 

,565 

42, 

,572 


.549 

126 

piO 

7 

111 

0 

0 

- 

.072 

- 

.254 

-3, 

,205 

3, 

,199 


.135 

127 

P~- 

7 

-2212 

0 

0 

- 

.701 

- 

.052 

-18, 

,281 

18, 

,903 


.938 

128 

piO 

7 

111 

0 

0 

- 

.179 


.069 

-12, 

,277 

12, 

,622 


.135 

129 

piO 

7 

111 

0 

0 

- 

.060 

- 

.060 

-12, 

,901 

13, 

,269 


.135 

130 

piO 

7 

111 

0 

0 

- 

.225 

- 

.220 

-17, 

,816 

18, 

,383 


.135 

131 

(Beam-REM) 

11 

92 

0 

0 

- 

.538 


.372 

74, 

,393 

74, 

,429 

2 

.212 

132 

(f_2) 

17 

225 

0 

0 

- 

.685 


.038 

30, 

,734 

30, 

,768 

1 

.273 

133 

piO 

7 

111 

0 

0 


.150 


.332 

43, 

,516 

45, 

,406 


.135 

134 

pi- 

7 

-211 

0 

0 

- 

.479 


.511 

9, 

,266 

9, 

,813 


.140 

135 

pit 

7 

211 

0 

0 

- 

.201 

- 

.476 

21, 

,182 

22, 

,184 


.140 

136 

(Beam-REM) 

11 

92 

0 

0 

- 

.016 


.852 

11, 

,680 

11, 

,879 

1 

.995 

137 

(K*_l+) 

17 

20323 

0 

0 

- 

.298 


.808 

8, 

,322 

8, 

,473 

1 

.340 

138 

(K~0) 

17 

-311 

0 

0 


.282 


.044 

3, 

,358 

3, 

,406 


.498 

139 

pit 

7 

211 

0 

0 

- 

.201 


.595 

2, 

,569 

2, 

,899 


.140 

140 

(K*0) 

17 

313 

0 

0 

- 

.094 


.211 

5, 

,610 

5, 

,686 


.892 

141 

K_L0 

7 

130 

0 

0 


.284 


.043 

3, 

,215 

3, 

,542 


.498 

142 

pi- 

7 

-211 

0 

0 

- 

.263 

- 

.066 

2, 

,426 

2, 

,691 


.140 

143 

K+ 

7 

321 

0 

0 


.174 


.274 

2, 

,899 

3, 

,222 


.494 

144 

(Beam-REM) 

11 

92 

0 

0 

1 

.284 

- 

.646 

18, 

,839 

19, 

,323 

4 

.053 

145 

(f_2) 

17 

225 

0 

0 


.927 

-1 

.937 

15, 

,847 

16, 

,042 

1 

.273 

146 

pi- 

7 

-211 

0 

0 


.360 

1 

.290 

2, 

,849 

3, 

,412 


.140 

147 

pi- 

7 

-211 

0 

0 


.541 

-1 

.029 

4, 

,275 

4, 

,752 


.140 

148 

pit 

7 

211 

0 

0 


.391 

- 

.911 

11, 

,286 

11, 

,931 


.140 

149 

(Beam-REM) 

11 

92 

0 

0 

- 

.328 

- 

.501 

4, 

,561 

4, 

,836 

1 

.491 

150 

(eta) 

17 

221 

0 

0 


.101 

- 

.808 

2, 

,154 

2, 

,367 


.549 

151 

pit 

7 

211 

0 

0 

- 

.426 


.306 

2, 

,265 

2, 

,568 


.140 

152 

piO 

7 

111 

0 

0 


.078 

- 

.224 


,225 


,457 


.135 

153 

pi- 

7 

-211 

0 

0 

- 

.073 

- 

.350 


,820 

1, 

,078 


.140 

154 

pit 

7 

211 

0 

0 


.105 

- 

.238 


,680 


,909 


.140 

155 

(Beam-REM) 

11 

92 

0 

0 

- 

.145 


.122 

2, 

,373 

2, 

,824 

1 

.518 

156 

(KO) 

17 

311 

0 

0 


.108 


.164 

1, 

,292 

1, 

,399 


.498 

157 

(K*-) 

17 

-323 

0 

0 

- 

.254 

- 

.043 

1, 

,081 

1, 

,425 


.892 

158 

(K_S0) 

17 

310 

0 

0 


.108 


.164 

1, 

,292 

1, 

,399 


.498 

159 

pi- 

7 

-211 

0 

0 

- 

.328 


.128 


,432 


,718 


.140 

160 

(K~0) 

17 

-311 

0 

0 


.077 

- 

.172 


,506 


,734 


.498 

161 

pi- 

7 

-211 

0 

0 


.018 

- 

.010 

-, 

,064 


,161 


.140 
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162 pi+ 

7 

211 

0 

0 

.096 

.171 

1.071 

1.287 

.140 

163 K_LO 

7 

130 

0 

0 

.080 

-.174 

.363 

.764 

.498 

164 (Beam-REM) 

11 

92 

0 

0 

-.335 

-.537 

2.570 

3.545 

2.358 

165 (rho+) 

17 

213 

0 

0 

-.239 

-.514 

.111 

.962 

.769 

166 (omega) 

17 

223 

0 

0 

-.097 

-.023 

2.460 

2.583 

.783 

167 piO 

7 

111 

0 

0 

.245 

-.072 

-.118 

.287 

.135 

168 pi+ 

7 

211 

0 

0 

-.478 

-.445 

-.057 

.702 

.140 

169 piO 

7 

111 

0 

0 

.120 

-.032 

.138 

.334 

.135 

170 pi- 

7 

-211 

0 

0 

-.147 

-.207 

.861 

1.087 

.140 

171 pi+ 

7 

211 

0 

0 

-.061 

.212 

1.032 

1.252 

.140 

172 (Beam-REM) 

11 

92 

0 

0 

-.022 

1.276 

2.371 

3.452 

2.159 

173 piO 

7 

111 

0 

0 

-.936 

.564 

.746 

1.475 

.135 

174 (rho-) 

17 

-213 

0 

0 

.917 

.711 

1.482 

2.033 

.769 

175 piO 

7 

111 

0 

0 

.150 

-.108 

-.042 

.247 

.135 

176 pi- 

7 

-211 

0 

0 

.773 

.816 

1.239 

1.857 

.140 

177 (Beam-REM) 

11 

92 

0 

0 

.028 

-.987 

1.765 

2.579 

1.601 

178 (omega) 

17 

223 

0 

0 

.159 

-.572 

.855 

1.303 

.783 

179 (omega) 

17 

223 

0 

0 

-.131 

-.414 

.910 

1.277 

.783 

180 piO 

7 

111 

0 

0 

.142 

.023 

.159 

.360 

.135 

181 pi- 

7 

-211 

0 

0 

-.177 

-.323 

.175 

.527 

.140 

182 pi+ 

7 

211 

0 

0 

.203 

-.276 

.092 

.436 

.140 

183 piO 

7 

111 

0 

0 

-.257 

-.051 

-.044 

.313 

.135 

184 pi- 

7 

-211 

0 

0 

-.004 

-.213 

.396 

.619 

.140 

185 pi+ 

7 

211 

0 

0 

.139 

-.155 

.129 

.368 

.140 

186 (Beam-REM) 

11 

92 

0 

0 

1.076 

.165 

1.776 

2.647 

1.633 

187 pi+ 

7 

211 

0 

0 

.429 

-.147 

-.255 

.504 

.140 

188 (rhoO) 

17 

113 

0 

0 

.650 

.311 

1.888 

2.162 

.769 

189 pi- 

7 

-211 

0 

0 

.471 

-.021 

.292 

.680 

.140 

190 pi+ 

7 

211 

0 

0 

.185 

.329 

1.310 

1.568 

.140 

191 (Beam-REM) 

11 

92 

0 

0 

.627 

-.168 

.968 

1.814 

1.390 

192 (rhoO) 

17 

113 

0 

0 

.136 

-.396 

.264 

.915 

.769 

193 pi- 

7 

-211 

0 

0 

.494 

.227 

.561 

.935 

.140 

194 pi- 

7 

-211 

0 

0 

.340 

-.491 

.007 

.654 

.140 

195 pi+ 

7 

211 

0 

0 

-.198 

.092 

-.029 

.286 

.140 

196 (Beam-REM) 

11 

92 

0 

0 

-1.186 

-.195 

1.298 

2.721 

2.067 

197 piO 

7 

111 

0 

0 

-.368 

.146 

-.593 

.641 

.135 

198 (rhoO) 

17 

113 

0 

0 

-.815 

-.343 

1.748 

2.105 

.769 

199 pi- 

7 

-211 

0 

0 

-.417 

-.533 

1.092 

1.472 

.140 

200 pi+ 

7 

211 

0 

0 

-.392 

.188 

.370 

.717 

.140 

201 (Beam-REM) 

11 

92 

0 

0 

-.184 

-.529 

.592 

1.588 

1.362 

202 piO 

7 

111 

0 

0 

.011 

.329 

-.331 

.403 

.135 

203 pi+ 

7 

211 

0 

0 

-.188 

-.861 

.637 

1.232 

.140 

204 (Beam-REM) 

11 

92 

0 

0 

.879 

.300 

.513 

2.125 

1.841 

205 (rho-) 

17 

-213 

0 

0 

.937 

.336 

.136 

1.265 

.769 

206 (rho+) 

17 

213 

0 

0 

-.058 

-.036 

.377 

.859 

.769 

207 piO 

7 

111 

0 

0 

.523 

.541 

-.130 

.794 

.135 

208 pi- 

7 

-211 

0 

0 

.420 

-.208 

-.021 

.521 

.140 

209 piO 

7 

111 

0 

0 

.103 

-.135 

-.307 

.270 

.135 

210 pi+ 

7 

211 

0 

0 

-.155 

.096 

.399 

.613 

.140 

211 (Beam-REM) 

11 

92 

0 

0 

.116 

-1.151 

.357 

2.408 

2.081 

212 (a_20) 

17 

115 

0 

0 

-.476 

-.521 

.333 

1.532 

1.318 

213 pi- 

7 

-211 

0 

0 

.595 

-.632 

-.119 

.912 

.140 

214 piO 

7 

111 

0 

0 

-.037 

.290 

.182 

.459 

.135 

215 (eta) 

17 

221 

0 

0 

-.435 

-.812 

.008 

1.073 

.549 

216 piO 

7 

111 

0 

0 

-.308 

-.509 

-.196 

.637 

.135 

217 piO 

7 

111 

0 

0 

-.048 

-.163 

-.031 

.244 

.135 

218 piO 

7 

111 

0 

0 

-.070 

-.145 

-.194 

.216 

.135 

219 (Targ-REM) 

11 

93 

0 

0 

-.558 

-.387 

-.074 

1.754 

1.615 

220 pit 

7 

211 

0 

0 

-.079 

.095 

.024 

.252 

.140 

221 (a_2-) 

17 

-215 

0 

0 

-.476 

-.483 

-.241 

1.502 

1.318 

222 piO 

7 

111 

0 

0 

.011 

-.471 

.009 

.532 

.135 

223 (rho-) 

17 

-213 

0 

0 

-.484 

-.014 

-.393 

.990 

.769 

224 piO 

7 

111 

0 

0 

-.074 

-.305 

-.459 

.465 

.135 
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225 pi- 

7 

-211 

0 

0 

-.405 

.289 

-.220 

.546 

.140 

226 (Targ-REM) 

11 

93 

0 

0 

-.584 

.922 

-.910 

2.344 

1.863 

227 (omega) 

17 

223 

0 

0 

-.237 

-.080 

-.504 

.964 

.783 

228 (omega) 

17 

223 

0 

0 

-.347 

1.002 

-.406 

1.379 

.783 

229 piO 

7 

111 

0 

0 

-.037 

.207 

-.270 

.282 

.135 

230 pi- 

7 

-211 

0 

0 

-.286 

-.188 

-.360 

.430 

.140 

231 pi+ 

7 

211 

0 

0 

.095 

-.102 

-.302 

.252 

.140 

232 piO 

7 

111 

0 

0 

-.286 

.288 

-.157 

.431 

.135 

233 pi- 

7 

-211 

0 

0 

-.047 

.065 

-.114 

.165 

.140 

234 pi+ 

7 

211 

0 

0 

-.006 

.645 

-.563 

.814 

.140 

235 (Targ-REM) 

11 

93 

0 

0 

-.472 

-1.489 

-1.060 

2.438 

1.543 

236 (rhoO) 

17 

113 

0 

0 

-.075 

-1.404 

-1.183 

1.992 

.769 

237 pi+ 

7 

211 

0 

0 

-.393 

-.087 

-.020 

.446 

.140 

238 pi- 

7 

-211 

0 

0 

-.109 

-.365 

-.137 

.405 

.140 

239 pi+ 

7 

211 

0 

0 

.039 

-1.041 

-1.332 

1.650 

.140 

240 (Targ-REM) 

11 

93 

0 

0 

.263 

.697 

-1.440 

2.287 

1.613 

241 (f_2) 

17 

225 

0 

0 

.163 

.602 

-.823 

1.639 

1.273 

242 pi- 

7 

-211 

0 

0 

.103 

.094 

-.760 

.674 

.140 

243 pi- 

7 

-211 

0 

0 

.404 

-.097 

-.916 

.923 

.140 

244 pi+ 

7 

211 

0 

0 

-.235 

.696 

-.193 

.782 

.140 

245 (Targ-REM) 

11 

93 

0 

0 

-.514 

-.122 

-1.450 

2.061 

1.366 

246 (omega) 

17 

223 

0 

0 

-.621 

.269 

-1.015 

1.450 

.783 

247 pi+ 

7 

211 

0 

0 

.110 

-.393 

-.578 

.636 

.140 

248 piO 

7 

111 

0 

0 

-.181 

.225 

-.261 

.342 

.135 

249 pi- 

7 

-211 

0 

0 

-.157 

-.169 

-.463 

.419 

.140 

250 pi+ 

7 

211 

0 

0 

-.274 

.209 

-.720 

.716 

.140 

251 (Targ-REM) 

11 

93 

0 

0 

.521 

1.161 

-2.393 

3.205 

1.711 

252 piO 

7 

111 

0 

0 

.194 

.629 

-.845 

1.011 

.135 

253 (a_20) 

17 

115 

0 

0 

.329 

.530 

-1.691 

2.233 

1.318 

254 pi+ 

7 

211 

0 

0 

.155 

.570 

-.573 

.774 

.140 

255 (rho-) 

17 

-213 

0 

0 

.177 

-.042 

-1.261 

1.488 

.769 

256 piO 

7 

111 

0 

0 

.214 

-.025 

-1.454 

1.388 

.135 

257 pi- 

7 

-211 

0 

0 

-.031 

-.019 

-.093 

.153 

.140 

258 (Targ-REM) 

11 

93 

0 

0 

.406 

-.313 

-2.421 

3.135 

1.925 

259 (f_l) 

17 

20223 

0 

0 

-.212 

-.339 

-1.807 

2.252 

1.283 

260 pi- 

7 

-211 

0 

0 

.621 

.025 

-.756 

.918 

.140 

261 (eta) 

17 

221 

0 

0 

.066 

-.030 

-.883 

1.042 

.549 

262 pi- 

7 

-211 

0 

0 

.004 

-.207 

-.974 

.902 

.140 

263 pi+ 

7 

211 

0 

0 

-.277 

-.105 

-.237 

.343 

.140 

264 piO 

7 

111 

0 

0 

.143 

-.029 

-.607 

.524 

.135 

265 piO 

7 

111 

0 

0 

.032 

-.032 

-.434 

.324 

.135 

266 piO 

7 

111 

0 

0 

-.100 

.027 

-.270 

.214 

.135 

267 (Targ-REM) 

11 

93 

0 

0 

.900 

-.332 

-6.866 

7.556 

3.006 

268 (Delta++) 

17 

2224 

0 

0 

-.308 

-.535 

-3.638 

3.890 

1.232 

269 (Delta'-) 

17 

-2214 

0 

0 

1.209 

.203 

-3.227 

3.665 

1.232 

270 pit 

7 

211 

0 

0 

.096 

-.090 

-1.275 

1.194 

.140 

271 p+ 

7 

2212 

0 

0 

-.399 

-.448 

-2.649 

2.852 

.938 

272 piO 

7 

111 

0 

0 

.571 

.009 

-1.308 

1.355 

.135 

273 p~- 

7 

-2212 

0 

0 

.644 

.191 

-2.206 

2.457 

.938 

274 (Targ-REM) 

11 

93 

0 

0 

-.423 

2.411 

-7.871 

8.635 

2.572 

275 (K*0) 

17 

313 

0 

0 

.285 

1.878 

-6.966 

7.276 

.892 

276 K- 

7 

-321 

0 

0 

-.705 

.532 

-1.047 

1.414 

.494 

277 pi- 

7 

-211 

0 

0 

.385 

.806 

-3.537 

3.652 

.140 

278 K+ 

7 

321 

0 

0 

-.095 

1.069 

-3.716 

3.914 

.494 

279 (Targ-REM) 

11 

93 

0 

0 

.454 

.090 

-6.177 

6.678 

2.496 

280 (a_2+) 

17 

215 

0 

0 

.653 

.059 

-6.253 

6.424 

1.318 

281 piO 

7 

111 

0 

0 

-.196 

.030 

-.067 

.254 

.135 

282 piO 

7 

111 

0 

0 

.231 

.367 

-1.287 

1.280 

.135 

283 (rho+) 

17 

213 

0 

0 

.425 

-.310 

-5.108 

5.192 

.769 

284 piO 

7 

111 

0 

0 

.205 

-.443 

-4.489 

4.550 

.135 

285 pit 

7 

211 

0 

0 

.227 

.131 

-.905 

.850 

.140 

286 (Targ-REM) 

11 

93 

0 

0 

.279 

-1.845 

-8.534 

9.228 

2.973 

287 (eta) 

17 

221 

0 

0 

-.507 

-1.723 

-5.714 

6.014 

.549 
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288 

(a_2-) 

17 

-215 

0 

0 


.786 

- 

.121 

-2 

.821 

3. 

.214 

1 

.318 

289 

piO 

7 

111 

0 

0 

- 

.050 

- 

.557 

-1 

.864 

1. 

.887 


.135 

290 

pi- 

7 

-211 

0 

0 

- 

.101 

- 

.501 

-2 

.002 

2. 

.010 


.140 

291 

pi+ 

7 

211 

0 

0 

- 

.347 

- 

.669 

-2 

.277 

2. 

.358 


.140 

292 

piO 

7 

111 

0 

0 

- 

.054 

- 

.298 

-1 

.183 

1. 

.135 


.135 

293 

(rho-) 

17 

-213 

0 

0 


.843 


.175 

-1 

.780 

2. 

.122 


.769 

294 

piO 

7 

111 

0 

0 

- 

.048 

- 

.113 

- 

.299 


.240 


.135 

295 

pi- 

7 

-211 

0 

0 


.897 


.286 

-1 

.767 

1. 

.956 


.140 

296 

(Targ-REM) 

11 

93 

0 

0 


.537 

- 

.194 

-4 

.665 

4. 

.841 

1 

.163 

297 

piO 

7 

111 

0 

0 


.379 


.133 

-2 

.137 

2. 

.119 


.135 

298 

(rhoO) 

17 

113 

0 

0 


.161 

- 

.328 

-2 

.671 

2. 

.803 


.769 

299 

pi- 

7 

-211 

0 

0 


.140 


.210 

-1 

.211 

1. 

.150 


.140 

300 

pi+ 

7 

211 

0 

0 


.027 

- 

.541 

-1 

.746 

1. 

.765 


.140 

301 

(Targ-REM) 

11 

93 

0 

0 

-1 

.516 


.775 

-9 

.672 

10. 

.040 

2 

.086 

302 

K+ 

7 

321 

0 

0 

- 

.822 


.830 

-2 

.936 

3. 

.191 


.494 

303 

(K*-) 

17 

-323 

0 

0 

- 

.691 

- 

.056 

-6 

.879 

6. 

.972 


.892 

304 

pi- 

7 

-211 

0 

0 

- 

.067 


.030 

-3 

.431 

3. 

.424 


.140 

305 

(K~0) 

17 

-311 

0 

0 

- 

.622 

- 

.087 

-3 

.591 

3. 

.680 


.498 

306 

(K_S0) 

17 

310 

0 

0 

- 

.622 

- 

.087 

-3 

.591 

3. 

.680 


.498 

307 

piO 

7 

111 

0 

0 

- 

.373 

- 

.224 

-2 

.763 

2. 

.766 


.135 

308 

piO 

7 

111 

0 

0 

- 

.243 


.134 

-1 

.114 

1. 

.061 


.135 

309 

(Targ-REM) 

11 

93 

0 

0 

- 

.684 

- 

.224 

-7 

.561 

7. 

.791 

1 

.736 

310 

(K*+) 

17 

323 

0 

0 

- 

.060 

- 

.267 

-5 

.535 

5. 

.613 


.892 

311 

(K~0) 

17 

-311 

0 

0 

- 

.624 


.043 

-2 

.026 

2. 

.178 


.498 

312 

piO 

7 

111 

0 

0 


.214 

- 

.170 

-3 

.101 

3. 

.092 


.135 

313 

K+ 

7 

321 

0 

0 

- 

.268 

- 

.100 

-2 

.720 

2. 

.745 


.494 

314 

(K_SO) 

17 

310 

0 

0 

- 

.624 


.043 

-2 

.026 

2. 

.178 


.498 

315 

piO 

7 

111 

0 

0 

- 

.503 


.096 

-1 

.122 

1. 

.159 


.135 

316 

piO 

7 

111 

0 

0 

- 

.115 

- 

.055 

-1 

.190 

1. 

.106 


.135 

317 

(Targ-REM) 

11 

93 

0 

0 

1 

.053 

-1 

.107 

-7 

.787 

8. 

.016 

1 

.137 

318 

(rhoO) 

17 

113 

0 

0 


.693 

- 

.993 

-6 

.816 

6. 

.965 


.769 

319 

pi- 

7 

-211 

0 

0 


.363 

- 

.115 

-1 

.114 

1. 

.093 


.140 

320 

pi- 

7 

-211 

0 

0 


.194 

- 

.534 

-1 

.738 

1. 

.767 


.140 

321 

pi+ 

7 

211 

0 

0 


.504 

- 

.462 

-5 

.364 

5. 

.477 


.140 

322 

(Targ-REM) 

11 

93 

0 

0 


.585 


.227 

-23 

.490 

23. 

.597 

2 

.149 

323 

K+ 

7 

321 

0 

0 


.796 


.530 

-7 

.564 

7. 

.799 


.494 

324 

(K*-) 

17 

-323 

0 

0 

- 

.208 

- 

.305 

-16 

.069 

16. 

.098 


.892 

325 

pi- 

7 

-211 

0 

0 

- 

.071 


.040 

-9 

.208 

9. 

.429 


.140 

326 

(K~0) 

17 

-311 

0 

0 

- 

.135 

- 

.346 

-7 

.003 

7. 

.031 


.498 

327 

K_LO 

7 

130 

0 

0 

- 

.132 

- 

.348 

-7 

.146 

7. 

.312 


.498 

328 

(Targ-REM) 

11 

93 

0 

0 

-1 

.689 


.202 

-60 

.528 

60. 

.729 

4 

.635 

329 

(Delta++) 

17 

2224 

0 

0 

-1 

.569 

-1 

.491 

-21 

.311 

21. 

.456 

1 

.232 

330 

(Delta~—) 

17 

-2224 

0 

0 

- 

.120 

1 

.693 

-39 

.217 

39. 

.273 

1 

.232 

331 

pi+ 

7 

211 

0 

0 

- 

.301 

- 

.495 

-4 

.010 

4. 

.070 


.140 

332 

P+ 

7 

2212 

0 

0 

-1 

.262 

- 

.999 

-17 

.587 

18. 

.244 


.938 

333 

pi- 

7 

-211 

0 

0 

- 

.226 


.567 

-12 

.176 

12. 

.531 


.140 

334 

P~- 

7 

-2212 

0 

0 


.112 

1 

.123 

-27 

.327 

28. 

.312 


.938 

335 

(Targ-REM) 

11 

93 

0 

0 

- 

.196 


.168 

-36 

.959 

37. 

.045 

2 

.502 

336 

(rho+) 

17 

213 

0 

0 

- 

.263 


.404 

-5 

.474 

5. 

.549 


.769 

337 

(omega) 

17 

223 

0 

0 


.067 

- 

.237 

-31 

.485 

31. 

.496 


.783 

338 

piO 

7 

111 

0 

0 

- 

.342 


.522 

-5 

.018 

5. 

.113 


.135 

339 

pi+ 

7 

211 

0 

0 


.085 

- 

.120 

- 

.742 


.657 


.140 

340 

piO 

7 

111 

0 

0 


.309 


.029 

-10 

.513 

10. 

.790 


.135 



sum: 

.00 




.000 


.000 


.000 

546. 

.000 

546 

.000 
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E Summary list of subroutines 


In the following list, the notation is S = subroutine, F = function B = block data. 

PART 1: Main steering rontines and general utilities 
S VNICLOC to measure and return the CPU time used 
S VNIXRIN to initialize the overall run and simulation procedure 
S VNIXRUN to steer the simulation of a single collision event 
S VNIXFIN to finish up the overall run and simulation procedure 
S VNICBIN to interface commonblocks of VNI and JESTSET/PYTHIA 
S VNILOGO to write program header 

S VNILIST to list event record, particle, or parameter values 
S VNIEDIT to perform global manipulations on the event record. 

F KVNI to provide integer-valued event information 
F PVNI to provide real-valued event information 
F RLU to provide a random number generator 
S RLUGET to save the state of the random number generator 
S RLUSET to set the state of the random number generator 
S VNIERRM to write error messages and abort faulty run 
S VNIFILE to open and close output files 
S VNIDUMP to dump(read) all commonblocks to(from) file 
B VNIDATA to give default values to particle and decay data 

PART 2: Initialization routines 

general as well as proeess-dependent ones 

S VNICSIN to initialize collision system and kinematics 
S VNIABXS to give parametrized cross-sections for beam/target 
S VNITMIN to initialize time slizes and form of time increment 
S VNI2JIN to set up an initial 2-jet in l~^l~ ^/Zq ^ qq 
S VNI2JFL to select flavour for initially produced qq 2-jet 
S VNI4JIN to set up an initial 4-jet in l~^l~ W~^W~ qq q'q' 

S VNI4JWW to initialize W~^ and W~ decay vertices and the 4 jets 
S VNILHIN to set up initial system for l~ + h ^ l~ -|-jet(s) -|- X 
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S VNILAIN to set up initial system for Z + A ^ I + jet(s) + X 
S VNIXDIS to generate the variables of DIS photon-hadron (nucleus) 
F VNIGAQX to give elementary photon-quark cross-sections 
S VNIGAQQ to give qq fluctuation of photon from VDM 
S VNIHHIN to set up initial system for h + h' ^ jets -|- X 
S VNIHAIN to set up initial system for h + A ^ jets -|- X 
S VNIABIN to set up initial system for A + A' ^ jets -|- X 
S VNIBTIN to set up initial state of hadronic or nuclear beams 
S VNINPDT to tabulate nuclear neutron/proton distribution 
S VNIRPIN to initialize single hadrons, or nucleons in nucleus 
S VNIQGDT to tabulate distributions in x and flavor of partons 
F VNIEFFQ to give parametrization of eff. parton resolution Q{x) 

F VNIQGSH to give simple parametrization of parton shadowing 
S VNIVQKF to extract the valence quark content of a hadron 
S VNIQGXF to initialize x and flavor of partons in parent hadron 
S VNIQGCS to arrange a set of partons as color singlet hadron 
S VNIQGIN to distribute initialized partons inside parent hadron 
S VNIPRKT to give primordial kX to initial partons 
S VNIXTIN to find spatial extension of initial particle system 

PART 3: Particle structure function package 

portable, self-contained program unit with interface to PDFLIB 

S VNISTFU calls structure functions; photon, tt, n,p, hyperon 
S VNISTGA gives photon structure function 
S VNISTAG gives anomalous part of photon structure function 
S VNISTGS parametrized parton distributions in a virtual photon 
S VNISTPI gives tt structure function 

S VNISTFL gives proton structure function at small x and/or 
S VNISTPR gives proton structure function 
S VNIFCTQ CTEQ2’94 parametrization of parton distributions 
S VNIFGRV GRV’94 parametrizations of parton distributions 
F FV Function for GRV parametrizations 
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F FW Function for GRV parametrizations 
F FWS Function for GRV parametrizations 
F FHS Function for GRV parametrizations 

F VNIHFTH gives threshold factor for heavy flavour production 
F VNIDILN gives Dilogarithm function 
F VNIGAMM gives ordinary Gamma function 
S PDFSET dummy routine, to be removed when PDFLIB is linked 
S STRUCTM dummy routine, to be removed when PDFLIB is linked 

PART 4: Main routines for cascade evolution and hadronization 

S VNIEV2J to perform evolution for + l~ ^ 2jets ^ hadrons 
S VNIEV4J to perform evolution for -\-l~ 4jets ^ hadrons 
S VNIEVLH to perform evolution for l~ + h ^ jets ^ hadrons 
S VNIEVLA to perform evolution for l~ + A ^ jets ^ hadrons 
S VNIEVAB to perform evolution for A + B ^ jets ^ hadrons 
S VNICOL2 to find 2-parton collisions and their kinematics 
S VNICOLX to give cross-sections for 2^2 and 2^1 collisions 
S VNICOLI to set up a found 2-parton collision subsystem 
S VNICOLL to find outgoing flavors and type of parton collision 
S VNISBRA to generate space-like parton branchings 
S VNITBRA to generate time-like parton branchings 
S VNICHAD to organize parton cluster formation and hadronization 
S VNICFOR to form one or two cluster(s) from a pair of partons 
S VNICREC to reconnect cluster-partons according to color flow 
S VNICFRA to fragment found parton clusters into hadrons 
S VNIBHAD to do beam remnant cluster formation and hadronization 
S VNIBFRA to fragment beam/target remnant clusters into hadrons 
S VNIPADC to decay unstable particles according to particle data 
S VNIPAKF to give a hadronic particle from a flavor pair 
S VNIBSTP to simulate baryon stopping in nuclear A-|-B collisions 

PART 5: Auxiliary routines I: 

parton cascades & cluster formation 
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F VNIALEM to give the aem value 
F VNIALPS to give the ag value 

F VNIANGL to give the angle from known x and y components 
F INVCHGE to give three times the electric charge 
F INVCOMP to compress standard KF flavour code to internal KC 
S VNICOLB to boost/rotate 2-body collision subsystem along pz 
S VNIDEC2 to generate 2-body decay 
S VNIFRAM to perform boosts between different frames 
S VNIKLIM to find kinematic limits in 2-parton collsions 
S VNICCOL to connect color indices in 2-parton collisions 
S VNIJCOL to connect partons with color flow indices 
S VNILCOL to label color indices according to color flow tracing 
F VNIMASS to give the mass of a particle or parton 
F VNIMTHR to give ‘current’ or ‘constituent’ masses for partons 
S VNINAME to give the name of a particle or parton 
F VNIGORI to give the origin of a particle in space-time history 
F VNISORI to sey the origin of a particle in space-time history 
S VNIPROP to propagate a given number of particles 
S VNIPSEA to generate a hadron pair from the Dirac ‘sea’ 

S VNIQSEA to generate a g or g, or a gg-pair from the ‘sea’ 

S VNIROBO to rotate and/or boost a given number of particles 
S VNIRSCL to rescale particle momenta for numerical precision 
S VNISORT to count and sort particles 
S VNIWIDT to calculate width of particle resonances 

PART 6: Auxiliary routines II; 

cluster-hadron decay and beam/target fragmentation 

S HWHEPC to convert VNI event record to or from standard HEP 
S HWLIST to list event record in HEP format 
S HWIGIN to set parameters of HERWIG 5.7 
S HWUINC to compute constants and lookup tables 
S HWUINE to perform re-initialization of a new event 
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S HWURES to load common blocks with particle data 
S HWUIDT to translate particle identifiers 
S HWCCUT to cut up heavy cluster into two clusters 
S HWCFLA to to set up flavors for cluster decay 
S HWCHAD to hadronize prepared beam/target clusters 
S HWDHAD to decay unstable hadrons from beam/target remnants 
S HWDTWO to generate two-body decay 0 ^ 1 -|- 2 
S HWDTHR to generate three-body decay 0 ^ 1 -|- 2 -|- 3 
F HWDPWT to give matrix element squared for phase-space decay 
F HWDWWT to give matrix element squared for Va weak decay 
F HWMNBI to compute negative binomial probability 
S HWMLPS to generate cylindrical phase-space 
S HWMULT to pick charged multiplicity of beam/target diffraction 
F HWUPCM to compute c.m. momentum for masses in 2-body decay 
S HWULOB to transform momenta from rest frame into lab frame 
S HWULOF to transform momenta from lab frame into rest frame 
S HWUMAS to put mass in 5th component of momentum vector 
S HWUROB to rotate vectors by inverse of rotation matrix R 
S HWUROT to perform rotation with rotation matrix R 
S HWUSOR to sort n-tupel A(N) into ascending order 
F HWUSQR to give square root with sign retention 
S HWUSTA to make a particle type stable 
S HWVDIF to give vector difference 
F HWVDOT to give vector dot product 
S HWVFQU to give vector equality 
S HWVSCA to give product of vector times scalar 
S HWVSUM to give vector sum 
S HWVZRO to give zero vector 
F HWRGFN to provide random number generator 
S HWRAZM to randomly rotate 2-vector p± 

F HWRFXP random number from 1-dim. exponential distribution 
F HWRFXT random number from 2-dim. exponential distribution 
F HWRGAU random number from 2-dim Gaussian distribution 
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F HWRINT random integer in finite interval 
F HWRLOG random logical 

S HWRPOW random number from power-law distribution 
F HWRUNG random number from flat-top distribution 
F HWRUNI uniform random in finite interval 
S HWUTIM dummy routine 
S HWWARN dummy routine 

B HWUDAT to provide HERWIG common blocks with data 

PART 7: Event study and data analysis routines 

S VNIANAL to steer pre-programmed event analysis 
S VNIANAl statistics on global particle properties & observables 
S VNIANA2 statistics on hadronic distributions 
S VNIANA3 statistics on pre-hadronic clusters 
S VNIANA4 statistics on partons and produced hadrons 
S VNIANA5 global thermodynamic multiparticle properties 
S VNISPHE to perform sphericity analysis 
S VNITHRU to perform thrust analysis 

S VNICLUS to perform three-dimensional jet cluster analysis 
S VNICELL to perform jet cluster analysis in (ry, (f), E±) 

S VNIJMAS to give high and low jet mass of event 
S VNIFOWO to give Fox-Wolfram jet moments 
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